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Abstract. This is a survey on eigenfunctions of the Laplacian on Riemannian manifolds 
(mainly compact and without boundary). We discuss both local results obtained by ana- 
lyzing eigenfunctions on small balls, and global results obtained by wave equation methods. 
Among the main topics are nodal sets, quantum limits, and norms of global eigenfunc- 
tions. The emphasis is on the connection between the behavior of eigenfunctions and the 
dynamics of the geodesic flow, reflecting the relation between quantum mechanics and the 
underlying classical mechanics. We also discuss the analytic continuation of eigenfunctions 
of real analytic Riemannian manifolds (Af , g) to the complexification of M and its applica- 
tions to nodal geometry. Besides eigenfunctions, we also consider quasi-modes and random 
linear combinations of eigenfunctions with close eigenvalues. Many examples are discussed. 
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Introduction 

The aim of this article is to survey both classical and recent results on the eigenfunctions 

(1) A, y,,^. = A? y.,^ 

of the (positive) Laplacian on a (mainly compact) Riemannian manifold (M, g). We con- 
centrate on the boundaryless case dM = for simplicity; when dM 7^ we impose standard 
boundary conditions. When {M,g) is compact, the spectrum is discrete and we arrange the 
eigenvalues in non-decreasing order Aq < Ai < A2 T 00. We denote by {fx^} an orthonormal 
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basis of eigenfunctions with respect to the inner product {ip\.,ip\^) = jj^^^(px^{x)ipx^{x)dVg. 
The 'topography' of an eigenfunction ideally encompasses the the shapes of pits and peaks in 
the graph of ip\, the geometry and connectivity of 'excursion sets' {x : ^\^{x) > h{\j)}, the 
norms and distribution function of (px^ , the distribution of its nodal sets and other level 
sets, the number of nodal components of different sizes, the number and distribution of its 
critical points, and the concentration, oscillation and vanishing order properties of ip\ ., often 
encapsulated by the so-called quantum limits (or microlocal defect measures), i.e. limits of 
quantum expectation values (Ayj^j , V^a^ ) as the eigenvalues tend to infinity. 

Eigenfunctions of Laplacians arise in physics as modes of periodic vibration of drums and 
membranes. They also represent stationary states of a free quantum particle on a Riemann- 
ian manifold. More generally, eigenfunctions of Schrodinger operators represent stationary 
energy states of atoms and molecules in quantum mechanics [Schj . The topography of modes 
of vibration and stationary states began with Ernst Chladni [Cl IC2j , who raised the prospect 
of 'visualizing sound' by bowing plates and observing the patterns of nodal lines (zero sets) 
of these modes. Over the last thirty years, chemists and physicists have used computers 
rather than bowed plates to visualize the energy states of atoms. Some computer graphics of 
eigenfunctions may be found in such articles as [SHMl [H] . In commemoration of the 200th 
anniversary of Chladni's diagrams, a recent volume |SS] has appeared which connects his 
work with that of contemporary physicists, chemists and mathematicians. 

In mathematics, studies of eigenfunctions tend to fall into two categories: 

• (i) analyses of ground states, i.e. ip\^ and ^jai; 

• (ii) analyses of high frequency (or semi-classical) limits of eigenfunctions, i.e. the 
limit as Xj —>■ oo. 

Behavior of ground states is very relevant to behavior of highly excited states, since an eigen- 
function ifx is always the ground state Dirichlet eigenfunction in any of its nodal domains. 
But our main emphasis in this survey is on the high frequency behavior of eigenfunctions 
rather than on the ground states. 

Studies of high frequency behavior eigenfunctions also fall into two categories: 

• Local results, which often hold for any solution of ([T]) on a (small) ball Br{x), often 
r = 0(A~^), irrespective of whether the eigenfunction extends to a global eigen- 
function on M. Doubling estimates and vanishing order estimates in terms of the 
frequency function, exponential decay bounds and nodal volume estimates often fall 
into this category. These methods often apply to large classes of functions: har- 
monic functions, poljTiomials, eigenfunctions, and more general solutions of elliptic 
equations. 

• Global results, which do use this global extension to M. The typical global assump- 
tion is that ipx is an eigenfunction of the wave group Ut = e**^. Global properties 
generally reflect the relation of the wave group and geodesic flow, particularly the 
long time behavior of waves and geodesies on the manifold. 

We aim to cover both sides of the subject. There already exists a very well written survey 
of the local aspects (the book of Q. Han and F.H. Lin [H]), so we give more details on 
the global aspects. But one of our purposes in this survey is to state the main local and 
global results so that the reader can compare approaches. As we have recently written a 
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survey on inverse spectral problems |Zlj . we concentrate on eigenfunctions and do not discuss 
eigenvalue asymptotics very much. 

The global behavior of eigenfunctions can only be fully understood by making a phase 
space analysis, where the phase space is the co-tangent bundle T*M or an energy surface 
S*M. For example, one often wishes to construct highly localized eigenfunctions or approxi- 
mate eigenfunctions (quasi-modes) of or to prove that they do not exist. Zonal spherical 
harmonics on the standard sphere (or any surface of revolution) are examples of eigenfunc- 
tions which are most highly localized at a point p (the poles). But it is more illuminating to 
observe that such modes are actually concentrating on certain Lagrangian submanifolds in 
phase space (namely, the one obtained by flowing out S*M under the geodesic flow). With 
this picture in mind, one would not expect to find modes which are highly concentrated at 
a point unless there exists a invariant Lagrangian manifold which has a large singularity 
over that point, and unless the eigenfunction concentrates in the phase space (micro local) 
sense on this Lagrangian manifold. One similarly expect modes which are extremal for lower 
l^p norms to concentrate along elliptic closed geodesies. In the case of integrable systems 
(see §12p one can prove that such expectations are correct. One further expects such special 
modes in the integrable case to be extremals for concentration. 

To obtain global phase space results relating the behavior of eigenfunctions to the behavior 
of geodesies, it is necessary to use microlocal analysis, i.e. the calculus of pseudo-differential 
and Fourier integral operators. Microlocal analysis is in part the mathematically precise 
formulation of the semi-classical limit in quantum mechanics: Pseudo-differential operators 
are 'quantizations' Op{a) of functions on the phase space T*M while Fourier integral oper- 
ators are quantizations of Hamiltonian flows (and more general canonical transformations). 
To motivate the phase space analysis, consider that one often studies the concentration and 
oscillation properties of eigenfunctions through the linear functionals px^^A) = {Aipx.,ip\.) 
on the space of zeroth order pseudo-differential operators A. The possible limits of the family 
{PAj }, known as quantum limits or microlocal defect measures, are probability measures on 
S*M which are invariant under the geodesic flow. It is difficult but important to determine 
or at least constrain the possible limits. This aim would not even be visible without a phase 
space perspective. 

We briefly review some of the basic methods and results of microlocal analysis in ^ Due 
to space limitations, we cannot go over the definitions and methods in much detail. For 
background in microlocal analysis, we refer to |EZl |DSj[ |GSj| and for systematic treatises 
to j HoI-IVj (see especially Volumes III and IV) and to [SVj (in particular for boundary 
problems). Other articles with expositions of the wave kernel of a Riemannian manifold 
are [Bel [Sul ITa2t [Z4] . Since microlocal methods and results are fundamental to the global 
theory, we use them in somewhat detailed arguments in the second half of this survey. 
Although the ideas and methods may be unfamiliar to some geometric analysts interested in 
eigenfunctions, they are quite geometric and concrete, and in a somewhat less precise form 
are common-place among physicists and chemists. While writing this survey, the author 
found it interesting to compare the range and power of the local and global methods, and 
hopes that experts in each side of the subject will find the other side of the story stimulating. 

We would like to thank H. Christianson, Q. Han, H. Hezari, D. Mangoubi, S. Nonnen- 
macher, M. Sodin, J. Toth and the referee for corrections on earlier versions of this survey. 



6 



STEVE ZELDITCH 



1. Basic Definitions and Notations 

The Laplacian of a Riemannian Riemannian (M, g) is the self-adjoint operator associ- 
ated to the Dirichlet form Q{f) = \df\'^gdVg where dVg is the volume form of {M,g) and 
I ■ Ig is the metric on one- forms. It is given in local coordinates by the expression, 

It is a natural geometric operator in the sense that commutes with all isometrics h of 
{M,g) where h acts on functions by translation f^{x) = f{h~^x). 

The eigenvalue problem ([T]) is a stationary version of the homogeneous wave equation 

(3) Oguit, x) = 0, where Ug = — + Ag. 

The eigenvalue problem ([T]) is equivalent to finding the periodic solutions e^'^'^^Aj (2;) in time. 
The reader should carefully note two conventions we use throughout this article: 

(1) We define the Laplacian to be positive in ([2]). This is opposite to the convention 
of many authors, but saves us from writing many minus signs in wave and heat 
operators; 

(2) We denote eigenvalues by in ([1]), so that A is the eigenvalue of a/A. This again 
saves many square root signs in wave operators. Physicists often denote A by k. 

1.1. Planck's constant and eigenvalue asymptotics. The eigenvalue problem for fixed 
A is a model elliptic equation. A standard fact (elliptic regularity) is that (px G C°°(M) for 
any C°° metric g and that ipx G C"^(M) (real analytic) when {M,g) is real analytic. 

However, it is more illuminating to regard A as the inverse of a semi-classical parameter 

A = h'\ 

where A is regarded as an 'operator' of order 1. Often one writes 

A-^A - 1 = /i^A - 1, 
in particular when constructing approximate eigenfunctions and normal forms. 

1.2. Spectral kernels. The individual eigenfunctions are very difficult to study directly. 
One generally approaches them through various kernel functions, i.e. Schwartz kernels of 
functions of A^. A basic one is the spectral projections kernel, 

(4) E{\,x,y)= ipj{x)ipj{y). 

j:Xj<X 

Semi-classical asymptotics is the study of the A — > 00 limit of the spectral data {ipj, \j} or of 
E{X,x,y). The (Schwartz) kernel of the wave group Ut = e**^ can be represented in terms 
of the spectral data by 

Utix,y) = U{t,x,y) := ^ e**^^(^j(a;)(^,(?/), 

j 

or equivalently as the Fourier transform J^e^^''^dE{X,x,y) of the spectral projections. (Note 
that we sometimes write Ut{x,y) and sometimes U{t,x,y) for the wave kernel, whichever 
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notation is more convenient in its context.) Hence spectral asymptotics is often studied 
through the large time behavior of the wave group. It is more or less equivalent to study the 
resolvent kernel 



for A e C lying on horizontal or logarithmic curves in C 

To obtain relations between geometry and eigenfunctions, it is necessary to give approx- 
imate formulae for these kernel functions in terms of geometric invariants. Such formulae 
originate in works of Hadamard and Riesz and are often called Hadamard-Riesz parametrices 
for manifolds without boundary. Other useful parametrices have been constructed by Lax 
and Hormander. The constructions are reviewed in ^ based on the exposition in [D.GllBe] . 
Another excellent reference is |HoI-IVj . 

We only discuss wave kernel parametrices for manifolds without boundary. Parametrices 
for the wave kernel in the boundary case are very difficult due to 'grazing rays'. For further 
discussion we refer to |HoI-IV[ ISV] . 

1.3. Geodesic flow. The geodesic flow is the Hamiltonian flow of \C,\g on the cotangent 
bundle T*M of M, equipped with its canonical symplectic form Yli dxi A d^i. By definition, 
g^{x,^) = {xt,C,t), where {xt,^t) is the terminal tangent vector at time t of the unit speed 
geodesic starting at x in the direction ^. Here and below, we often identify T*M with the 
tangent bundle TM using the metric to simplify the geometric description. The geodesic 
fiow preserves the energy surfaces {\^\g = E}, i.e. the co-sphere bundles S^M. Due to the 
homogeneity of H, the fiow on any energy surface {\C,\g = E} is equivalent to that on the 
co-sphere bundle S*M. 

This definition of the geodesic fiow makes homogeneous of degree one, and is slightly 
different from the geometer's definition as the Hamiltonian fiow of |^|^. The homogeneous 
geodesic flow is not well-deflned on the zero section of T*M, and so it is punctured out. The 
punctured cotangent bundle is denoted T*M\0. 

1.4. Closed geodesies. By a closed geodesic 7 one means a periodic orbit of (?* on S*M. 
The period is the length of 7 viewed as a curve in M (i.e. projected to M), thus 
g^"'{x,^) = (x,^) for all (x,^) G 7. The set of all periodic points {x,^), of all possible 
periods, is often called 'the set of closed geodesies'. A basic dichotomy in spectral theory is 
as follows: 

(1) Aperiodic; The Liouville measure of the closed orbits of (7* is zero; or 

(2) Positive measure of periodic orbits: the Liouville measure of the closed orbits is 



(3) Periodic : If the entire geodesic flow is periodic, g^ = id for some T > 0, the metric is 
said to be ZoU. The common Morse index of the T-periodic geodesies will be denoted 



In the real analytic case, {M,g) is automatically of type either (1) or (3), since a positive 
measure of closed geodesies implies that all geodesies are closed. In the C°° case, it is simple 
to construct examples with a positive but not full measure of closed geodesies (e.g. a pimpled 
sphere) . 




positive. 



by /3. 
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1.5. Jacobi fields and lineEir Poincare map along a closed geodesic. We recall that 
Jacobi's equation for a normal vector field Y along 7 is Y" + -R(T, Y)T = where T = 7 
and where R is the curvature tensor. We denote the space of normal Jacobi fields along a 
closed geodesic 7 by J'^. It is a symplectic vector space of dimension 2(dimM — 1) with 
respect to the Wronskian 

Mx,r) = s(x,£y)-s(£x,r). 

The linear Poincare map is the (real) linear symplectic monodromy map on (J'^,u) 
defined by PjY{t) = Y{t + L^). 

To diagonalize it, we complexify to obtain a linear complex symplectic map on the space 
JT"^ ® C of complex normal Jacobi fields. Since e Sp{J'^ ® C, a;) (the symplectic group), 
its spectrum cr(P^) is stable under inverse and complex conjugation: thus, if p e a{P^), 
then also p~^,p,p~^ e '^{P^)- The closed geodesic is non- degenerate if 

pT^i . . .p;^" = 1 =^ mi - (Vi,mi e N). 

In particular, the eigenvalues are simple and ±1 ^ cr(P^). It is called elliptic if all eigenvalues 
of Pj are of modulus one. We then denote them by {e^'^°'\j — 1, ...,n}. Thus, {ai, ...,«„}, 
together with tt, are independent over Q. The closed geodesic is called hyperbolic if all of 
the eigenvalues are real. They then come in inverse pairs and we denote them by {e^^^ ,j = 
l,...,n}. When dimM > 2, there are mixed hyperbolic and elliptic geodesies, and more 
general ones where the eigenvalues are complex and not of modulus one. For simplicity of 
exposition, we only consider elliptic and hyperbolic geodesies. 

The associated normalized eigenvectors will be denoted {Yj, Yj, j — 1, ...,n}, 

P^Yj = e'^^Yj P^Yj = e-^"^F,- u;(Yj,Yk) = 5jk 

and relative to a fixed parallel normal frame e{s) :— (ei(s), ...,e„(s)) along 7 they will be 
written in the form Yj{s) — Yl^=iyjk{s)^k{s). 

1.6. Geodesic fiow as a uniteiry operator. Since the Hamiltonian flow of \^\g preserves 
the canonical symplectic form cu — dxj A d$,j of T*M, it preserves the volume form a;"* 
where m = dimM. It also preserves the one form d\^\ and hence it preserves the Liouville 
form dfiL := on S*M. By definition, the quotient is the unique 2m — 1 form whose 
wedge product with d\^\g equals uj"^. 

We define the unitary operator Vf on L'^{S*M, d^i) 

(5) Vt{a):=aog\ 

It is sometimes called the Koopman operator associated to the geodesic flow. 

1.7. Spectrum and geodesic flow. One of the principal emphases in this survey is on 
the relations between the global dynamics of the geodesic flow and the eigenfunctions of the 
wave group. These relations have a long tradition in quantum mechanics, and are referred 
to as semi-classical analysis or study of the classical hmit. There exists a large speculative 
physics literature on relations between eigenfunctions and eigenvalues and the underlying 
classical Hamiltonian system. Very little is understood, even conjecturally, outside of the 
following model types of geodesic flow: 
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• Ergodic or weak mixing geodesic flows, reviewed in ^ |T5l More quantitative results 
hold if the geodesic flow is Anosov (e.g. rates of quantum ergodicity ^ IT5l entropy of 
quantum limits ^ }T6!) . 

• Integrable systems, reviewed in §1211141 

• KAM systems, i.e. small perturbations of integrable systems. There are few if any 
results on eigenfunctions, but there exist results on approximate eigenf unctions or 
quasi-modes [D, ICV21 |Pop| . 

• Special systems with a 'divided phase space' possessing an open set of periodic orbits 
[MUZ] . 

• Special cases where group theory is available (quotients of reductive, or solvable, or 
nilpotent groups by discrete subgroups). 

1.8. Ergodic, weak mixing and Anosov geodesic flows. Ergodicity and weak mixing 
are spectral conditions on the geodesic flow. The geodesic flow is called ergodic if the only 
invariant functions, Vtf = f with / G L'^{S*M, dfii) are the constant functions. Equivalently, 
1 is an eigenvalue of multiplicity one. It is called weak mixing if the spectrum of Vt is 
continuous on the orthogonal complement of the constant functions. Mixing systems add a 
'smoothness' assumption on the spectral measures and is also a spectral condition. 



1.8.1. Examples. 

• The most familiar examples are (M, g) of strictly negative sectional curvatures. Their 
geodesic flows are Anosov. Model examples are compact quotients of hyperbolic 
space. 

• In |BD] . many embedded surfaces in with ergodic geodesic flow are described. 
They may have any genus. There exist real analytic metrics on S"^ with ergodic 
geodesic flow. 

• Although we do not discuss manifolds with boundary in detail, there are many ex- 
amples of domains in M*" or other Riemannian manifolds with ergodic geodesic flow. 
The most famous are the Bunimovich stadium (a rectangle with semi-circular ends) 
and a Sinai billiard table (a rectangle or torus with a disc removed). References to 
the literature are given in [HZt IZZwt IGLj . where ergodicity on the quantum level is 



studied. 

We recall that a geodesic flow g^ is Anosov if, for each p G S*M, the tangent space TpS*M 



splits into g* invariant sub-bundles 

TpSlM = E-{p) © E'^ip) (BRXh{p) 

where E^ is the unstable subspace and E^ the stable subspace and where Xh is the Hamil- 
tonian vector field of the function H{x,^) = \^\g on T*M. The unstable Jacobian J"(p) at 
the point p is defined as the Jacobian of the map g~^, restricted to the unstable subspace at 
the point g^p: J^{p) = det {dg^^ui^gip^) (the unstable spaces at p and g^p are equipped with 
the induced Riemannian metric) 

1.9. Completely integrable geodesic flow. By a (classical) completely integrable system 
on T*M with dim M = n, we mean a set of n independent, C°° functions pi, . . . , pn, on T*M 
satisfying: 
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• {Pi^Pj} = for all I < i, j < n; 

• dpi A dp2 A ■ ■ ■ A dpn 7^ on an open dense subset of T*M. 
The associated moment map is defined by 

(6) V = {pu...,Pn):T*M^BcR''. 

We refer to to the set B as the 'image of the moment map,' and denote the set of regular 
values of V by Breg- When V is proper and 6 is a regular value of V, 

(7) V~\b)=A^'\b)U---UA^^^'\b), {be Breg) 

where each A^''\b) ^ T" is an n-dimensional Lagrangian torus. The number mci{b) = 
j^V~^{b) of orbits on the level set V~^{b) is constant on connected components of Breg and 
the moment map ([H]) is a fibration over each component with fiber ([7]). 
The Hamiltonians pj generate an action of MP' defined by 

(8) $t = exp tiHp^ o exp t2'B.p^ ■ ■ ■ o exp t„Hp„ . 

We denote the $t-orbits by M" ■ (x, and the isotropy group of (x, ^) by 1{x,i,)- When 
M" ■ (x, ^) is a compact Lagrangian orbit, then is a lattice of full rank in M", and is 

known as the 'period lattice', since it consists of the 'times' T G M" such that ^T\KW{h) — ^d,- 
By the Liouville-Arnold theorem, the orbits of $f are diffeomorphic to M.^ x T™ for some 
{k,m),k + m < n. In sufficiently small neighbourhoods of each component torus, 

A*^'^(fe), the Liouville-Arnold theorem also gives the existence of local action-angle variables 
\ In\6^\ ...,6lP) in terms of which the joint fiow of Ep-^,...,Ep^ is linearized. For 
convenience, we henceforth normalize the action variables jj^...,/^'' so that jj'^ = 0; j = 
1, n on the torus A^^\b). 

Some examples of integrable systems are as follows. 

1.9.1. T" = R"/Z". In the case of a fiat torus, the C,j variables are action variables and the 
moment map is the projection V : T*T" — » M". The joint eigenvalues of the action operators 
Ij = is the standard integer lattice Z" and rays are multiples of lattice points. 

1.9.2. S"^ . The classical action variables are Ji = Peix,^) = {^,-§q), known as the Clairaut 
integral, and I2 = \^\g- The moment map I = {Ii, I2) maps T*S'^ to the triangular cone 
{{x,y) ■.y>0, \x\<y}c 

1.9.3. Simple surfaces of revolution. These are metrics on S"^ for which the geodesic fiow and 
eigenfunctions are almost as simple as for the standard sphere. The key feature is that they 
are 'toric integrable'. 

To explain this notion, we assume that the action is given by rotations around the 
Xs-axis. We are not primarily interested in this 5*^ action but rather the 5*^ x M action on 
T*M\0 generated by rotations and by the geodesic fiow. The basic assumption defining 
simple surfaces of revolution is that the M x S*^ action simplifies to a (2-torus) action on 
T*S'^\0. That is, can define two global action variables /i,/2 with 27r-periodic Hamiltonian 
fiows such that = H{Ii,l2). If we write the metric as dr'^ + a{r)'^d6'^ in geodesic polar 
coordinates centered at a fixed point, which we visualize as the north pole, then a sufficient 
condition that the geodesic fiow be toric is that the distance function to the axis of revolution 
possesses precisely one local maximum. For instance, a convex surface of revolution is simple. 
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A simple surface of revolution possesses precisely one invariant closed geodesic, which we 
refer to as the equator. 

The first action function is the Clairaut integral generating the action, defined by 
P0{v) :— {v, To define the second action variable, we need to consider the moment map 

V = mg,pe) : T*S^ -^B:= {(61, 62) : I62I < a{r,)h] C M x M+ 

of the Hamiltonian M x 5'^-action defined by the geodesic flow and by rotation. The singular 
set of V is the closed conic set Z := {(ro, 6', 0,pe) : 9 G [0,27r),pe G IR}, i-e. Z is the cone 
through the equatorial geodesic (in either orientation). The map V\T*SgS'^-z is a trivial 
X S*^ bundle over the open convex cone Bo (the interior of B), so that V~^{b) is an 
invariant torus for the geodesic flow and for rotation. The second action function is the 
function hiV) given by 

I2{bi,b2) = - / \ bl--\^dr + \b2\ 

^ Jr-(b) V 

where r±{b) (with b = (61, 62) arc the extremal values of r on the annulus tt o V~^{b) (with 
IT : SgS'^ ^ S"^ the standard projection). 

These action variables are best thought of in the following classical way (Liouville, Jacobi, 
Arnold): For each b e Bo , let Hi{Fb,Z) denote the homology of the fiber V~^{b). This 
lattice bundle is trivial since B is contractible, so there exists a smoothly varying homology 
basis {71(6), 72(6)} G Hi(V^^{b),'Z) where 71 are the orbits of the action, and 72 is a 
fixed closed meridian 'Jm when b is on the center line M"*" • (1,0). Then the actions are the 
integrals of the 'action form' over this moving homology basis: 

r f 1 r+^^ I /i2 

h{b) = / ^dx^pe, h{b) = / idx=- Jbl - -^dr+\b2\. 

On the torus of meridians in S^S"^, the value of I2 equals ^ and it equals one on the equatorial 
geodesic. So extended, /i,/2 are smooth homogeneous functions of degree 1 on T*S'^, and 
generate 27r-periodic Hamilton flows. 

The pair X := {h:h) generate a global Hamiltonian torus {S^ x 5'^)-action commuting 
with the geodesic flow. The singular set of I equals Z :— {I2 — ipe}, corresponding to the 
equatorial geodesies. The map 

I : T*S^ -Z^To:= {(x, y) e M x M+ : |x| < y) 

is a trivial torus flbration. 

1.10. Quantum mechanics: Wave group and pseudo-differential operators. The 

quantization of the Hamiltonian is the square root \fK of the positive Laplacian, of 
{M,g). It generates the wave 

Ut = e''^, 

which is a group of Fourier integral operators which propagates singularities along geodesies. 
In particular, for flxed t,x, the singular support of the wave kernel e**^(t,x, •) hes on the 
distance sphere centered at x of radius t. This is the wave front of a spherical wave launched 
at X. Moreover, the wave front propagates outward along geodesies normal to the distance 
spheres. Thus, the 'singular directions' lie along the geodesic rays emanating from x. 
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Above, -\/A, resp. e**^ are defined by the spectral theorem: i.e. they have the same 
eigenfunctions (p\ as A on {M,g) and with eigenvalues e**^^ resp. Xj. 

1.10.1. Pseudo- differential operators and expected values. The relations between the wave 
group and geodesic flow are fundamental to the global analysis and will be discussed at 
greater length in ^ They originated in Dirac's Principles of Quantum Mechanics |Dirj as 
the relations between a classical Hamiltonian flow and its quantization as a unitary group 
on a Hilbert space. The theory of pseudo-differential and Fourier integral operators is the 
rigorous mathematical framework for the somewhat heuristic ideas of Dirac's book. We refer 
the reader to textbook treatments |DSj[ lEZt |GSj| and to the comprehensive treatise [HoI-IVj 



for background. 

We denote the class of pseudo-differential operators of order m on M by We fix a 
quantization a — >• Op{a) of pseudo-differential operators to functions a{x,$,) G C°°(T*M\0) 
which are symbols of order m in the sense of admitting poly-homogeneous expansions, a ~ 
Yl^o (^m-j where am-j is homogeneous of order m—j on T*M\0. It is the behavior at infinity 
and not at ^ = of a symbol which is important and one usually cuts off the homogeneous 
functions in a small ball around 0. The principal symbol (Ta of A = Op{a) is the leading 
homogeneous term = flm- 

On M™, Op{a) is given by the formula 



Op{a)fix)= / a(x,Oe^<^-^'^^/(l/)«. 



Equivalently, Op{a) is defined by its actions on exponentials, 

(9) Op(a)e'<"'«> = a(x,Oe'<"'^^ 

Since a is assumed to be a symbol, the right side is a perturbation expansion with leading 
term ao{x, 

^)ei(^.5>. When a is independent of ^ one has a multiplication operator and when 
it is independent of x one has a convolution operator. 

We also follow the notation of Dimassi-Sjoestrand |DSj| for operator classes: Given an 
open U C M", we say that a{x,^; K) G C°°{U x M") is in the symbol class S"^'''(f/ x M"), 
provided 

\d:dla{x,i-n)\<Capn-^{i + \i\f-\^\. 

We say that a G S^'^{U x M") provided there exists an asymptotic expansion: 

oo 

a(x, ^; h) ~ aj{x,^)K' , 

3=0 

valid for |^| > ^ > with aj{x,^) G S^^'^'^U x W) on this set. The associated h Kohn- 
Nirenberg quantization is given by 

Opn{a){x,y) = (271^)"" ! e'^^-y^^/'^ a{x,^;h) d^ 

Jr"- 

As is well-known, the definition can be globalized to M using a partition of unity. We denote 
this class by Opn{S"^'^){T* M) . The symbol of the composition is given by the usual formula: 
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Given a e S"''''"' and b G S"'^^''\ the composition Opn{a) o Opnib) = Opn{c) + in 
L'^{M) where locally, 

c(x,e; h) ~ r(™^+'"^) V ^ ' / (d?a) ■ (d^b). 

|a|=0 



For further details, we refer to DSj 



On a manifold one patches together such local expressions using a partition of unity. There 
is no unique way to do this, and in fact there is no unique definition of Op{a) on R'". The 
one above is the 'Kohn-Nirenberg' definition; other natural choices are the Weyl definition 
and the Friedrichs positive quantization. 

On special manifolds, such as symmetric spaces, one may define Op{a) using the Fourier 
transform of the symmetric space in place of the Euclidean Fourier transform in local coor- 
dinates. We will discuss the definition in the case of hyperbolic space. 

Pseudo-differential operators greatly enlarge the class of operators which can be used to 
test eigenf unctions for their properties. One often uses matrix elements {Aip\.,ipx^) for this 
purpose. The diagonal matrix element is viewed as the 'expected value of the observable A 
in the state in quantum mechanics. Such quadratic forms in ipx often arise in the study 
of eigenfunctions, usually with Op{a) being just multiplication by a function, or a gradient 
operator. It is useful to analyze quadratic expressions of a more general kind. 

1.11. Modes and quasi-modes. Quasi- modes are approximate eigenfunctions. It is im- 
portant to discuss quasi-modes along with modes (true eigenfunctions) for the following 
reasons: 

• Most theorems and proofs concerning modes also apply to quasi-modes. There are 
few techniques that distinguish modes from quasi-modes. 

• It is often possible to construct quasi-modes with special properties. In many appli- 
cations, they are just as useful as modes, and often have clearer geometric properties. 

There are several ways to define 'approximate eigenf unction'. The classical definition 
of Keller [K], Babich [j] (see also [El]), Laz utkin [Ll], Arnold [Xj] Ral ston [R^ [Ra2] . 
Guillemin-Weinstein |GW] . Colin de Verdiere |CV2] (see also Popov |Pop| ) and others is 
that a quasi-mode {ipk} of order zero is a sequence of L^-normalized functions satisfying 

(10) ||(A-/Xfe)V^fe|U2 = 0(1), 

for a sequence of quasi-eigenvalues ^k- By the spectral theorem it follows that there must 
exist true eigenvalues in the interval [iik — 5, yUfc + S\ for some 5 > 0. Moreover, if Ek^s denotes 
the spectral projection for the Laplacian corresponding to this interval, then 

WEk^si^k - i^kWv^ = 0{k~^). 

One can refine the definition by demanding that the remainder in (ITUl) is of order 0{n^^) 
and define a quasi-mode of order s by 

(11) ||(A-/i,)V^fc|U. = 0(/x,^). 
We refer to |CV2^ IZ9j for other modifications. 
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In the 'classical' work on quasi-modes of Babich, Lazutkin, Arnol'd, Ralston and others, 
quasi-modes are often constructed as oscillatory integrals (or Lagrangian quasi-modes) 

with special phases and amplitudes designed so that ||(— A — A^)'?/'a|| is small. An important 
example is the construction of a quasi-mode associated to a stable elliptic orbit, reviewed 
in ^ ITOl Another important quasi-mode of this type is the quasi-mode associated to the 
central rectangle of a Bunimovich stadium, or more precisely to the Lagrangian cylinder 
with boundary consisting of unit tangent vectors to the 'bouncing ball orbits' in the vertical 
direction in the rectangle. A recent study of such quasi-modes and of the possible behavior 
of actual modes is given in the articles [BZl IBZ2] of Burq-Zworski. However, the definition 
(ITU]) does not force quasi-modes to be Lagrangian and includes 'random' combinations of 
eigenfunctions with eigenvalues in a small interval around A^, which need not have such an 
oscillatory integral structure. 

An important example of such a quasi-mode is a sequence of "shrinking spectral projec- 
tions", i.e. the L^-normalized projection kernels 

^z^^^ ^ X[\,,\,+e,]{x,z) 
^ VX[A,,Aj+e,](2;,2;) 

with second point frozen at a point z E M and with width ej 0. Here, xiXj,Xj+ej]{x, z) is 
the orthogonal projection onto the sum of the eigenspaces Vx with A G [Xj, Xj + ej]. The zonal 
eigenfunctions of a surface of revolution are examples of such shrinking spectral projections 
for a sufficiently small ej, and when z is a partial focus such are generalizations of 

zonal eigenfunctions. On a general ZoU manifold, shrinking spectral projections of widths 
ej = 0{XJ^) are the direct analogues of zonal spherical harmonics, and are quasi-modes of 
order 1. 

1.12. Heuristics and intutions. There are several key intuitions to keep in mind from the 
outset: 

• Local intuition: Eigenfunctions of eigenvalue A^ resemble polynomials of degree ~ CA 
in terms of their local complexity and growth, e.g. vanishing order at zeros, volumes 
of nodal hypersurfaces, growth rates on small balls. 

• Global intuition: Eigenfunctions are stationary states of the quantization U{t) = 
expity/A of the geodesic flow. Their high-frequency limits A — >■ oo should reflect 
the dynamics of the classical geodesic flow. When the geodesic flow is integrable, 
eigenfunctions should localize on the invariant tori (or more correctly, on level sets 
of the moment map). When the geodesic flow is ergodic, eigenfunctions should be 
diffuse (i.e. not localize). 

• Modes versus quasi-modes and random waves: Most results about eigenfunctions 
apply to quasi-modes, i.e. linear combination of eigenfunctions with very close by 
eigenvalues. More precisely, when \Xj — Xk\ < In integrable cases, one can 
spectrally separate out true eigenfunctions from such 'random waves' but in general 
one cannot. In ergodic cases, eigenfunctions in many respects resemble random waves. 
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The key difficulty in relating classical to quantum mechanics (e.g. in quantum chaos) is 
that it involves a comparison between long-time dynamical properties of and Ut through 
the symbol map and similar classical limits. The classical dynamics defines the 'principal 
symbol' behavior of Ut and the 'error' UtAU^ — Op{aA ° g^) typically grows exponentially 
in time. This illustrates the ubiquitous 'exponential barrier' in the subject. The classical 
approximation is not clearly valid after the so-called 'Heisenberg time' (see ^ IT6l) . The articles 
O lANl lANK] have excellent discussions of these problems. 

1.13. Notational Index. 

• = 5'(^, [g'^] is the inverse matrix to [g,j]. 

• r = r{x,y) denotes the distance function of {M,g). 

• B{xo,r) C M denotes the geodesic ball of radius r centered at xq; 

• dVg denotes the volume density of {M,g)] 

• Q{x, y) denotes the volume density in normal coordinates at x, i.e. dVg{y) = 9(x, y)dy. 

• T*M is the cotangent bundle of M, and T*M\{0} is the puncture cotangent bundle 
where the zero section is deleted. 



• \^\g = \ T!'ij=i9'^{^)iiii ■ T*M\{0} R+ denotes the length of a (co)-vector. 



• The unit (co-) ball bundle is denoted B*M = {{x,^) : \^\g < 1}. Its boundary 
S*M = {\C,\g = 1} is the unit cosphere bundle; 

• H is the Liouville measure on S*M, i.e. the surface measure d/j, = induced by 
the Hamiltonian H = \^\g and by the symplectic volume measure dxdC, on T*M. 

• u is the linear functional on C{S*M) defined by co'(cr) = ^(^gl^^^ js*M^^I^- '^^^ same 
notation is used for the functional (state) on the algebra \E'°(M) of zeroth order 
pseudo-differential operators defined by uj{A) = uj{aA)- 

• g^ : T*M\{0} T*M\{0} denotes the geodesic flow, i.e. the Hamilton flow of l^l^. 

• 7 denotes a closed geodesic, i.e. closed orbit of g^ in S*M. Thus, 7(t) = g*{xo,^o) G 
S*M where g'"{xo,C,o) = {xq,C,o). L = is the period of the closed geodesic. By 
abuse of notation, we sometimes also use 7 to denote its projection to M, where 

is the length of 7. 

• Geodesic loops versus closed geodesies: Viewed as curves on M both satisfy a(0) = 
a{L), but the latter also satisfy a'(0) = a'{L). 

• inj(M, (yf) denote the injectivity radius. 

• A(M) denote the loopspace of M; 

• Q{M) denote the subset of closed geodesies in A(M); 

• Q[^j denote the set of closed geodesies in Q{M) whose free homotopy class is [7]; 



There are only a few Riemannian manifolds (M, g) where one has explicit formulae for 
eigenf unctions. In this section, we briefly review these examples. What they have in com- 
mon is that in each case is completely integrable, i.e. commutes with a maximal family of 
(pseudo-differential operators). The joint eigenfunctions of such quantum integrable Lapla- 
cians have very special properties reflecting the complete integrability of the geodesic flow of 




2. Explicitly solvable eigenfunctions 
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(M, g). Our intuition is that such integrable eigenfunctions should be models of extremal be- 
havior: for instance, extremal growth and concentration. Quantum integrable eigenfunctions 
will be discussed in depth in ^ IT2l Here, we only wish to go through the simplest examples. 

2.1. M". The eigenspaces of the Laplacian on M" are defined by 

where 5'(M") is the space of tempered distributions. Since eigenfunctions are C°° the tem- 
peredness only constrains the growth to be polynomial, i.e. rules out exponentially growing 
eigenfunctions such as e^'*'^'^^ 

Since the Euclidean motion group En commutes with the flat Laplacian A = A^n it 
preserves the eigenspaces. Hence the infinitesimal translations and the infinitesimal 
rotations commute with A. The joint complex valued eigenfunctions of the translations are 
the Euclidean plane waves e*^^'^^ with ^ G M". The eigenspaces £x are infinite dimensional 
but are spanned in the following sense by the e*^^'^^ with |^| = A (i.e. of frequency A): 
There is a Poisson type integral formula for eigenfunctions: for any (px ^ £\ there exists a 
distribution d/i e V'{S"'~^) such that 

(12) ^x{x) = [ e^^«'-)rf/.(0. 

We refer to }Hel] for further background. 

Let (fx E Sx. Since SO{n) acts on Sx we may decompose it into isotypic subspaces, 

OD 

N=0 

where Sx{N) is the subspace of eigenfunctions transforming by the A^th irreducible repre- 
sentation of SO{n), realized by the space of spherical harmonics of degree on 5*""^ 
(see §2.31) . As this implies, the elements of Sx{N) may be expressed as sums of of separation 
of variable eigenfunctions JN^n{Xr)ipM{u) where ipN G T-Cn{S"'~^). The radial factor is a 
temperate solution of the n dimensional spherical Bessel equation 

^'t^ + - 1)^T- + - + ^ - '^JV,n(Ar) = 0, 
dr"^ ar / 



and is given explicitly by 

Jjv,„(Ar) = (Ar)-"^ J|^n (Ar) = C^^n^ (Ar)^ / cos(Ar cos 9) sin2^+'^-2 9d9, 

^ ^ Jo 

where C^_,_n^ = „_2 _ \ — ^ and where Ju is the Bessel function of order z/. 



2 2^r(2^+i)r(i) 



/•TT 

Jn+^ (rA) = C^+n^ (Ar)^+'^ / cos(Ar cos 9) sm^^^""^ 9d9. 
^ ^ Jo 

The space Sx carries an inner product which is invariant under En+i- An ortho normal 
basis is given by 

(13) u^x) = [ e^"«'^>Xa(O^K.(0, 



n-l 
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where {Xa} is an orthonormal basis of spherical harmonics of L'^{S"~^^). The notation is 
simplest when n = 2: Then an orthonormal basis of S\ is given by {e^ := j£{Xr)e^^^}(^_^. 
Equivalently, the norm of u in Sx is the norm of the distribution dfi in L^(S'"~^). We denote 
by S^^ the subspace of elements of finite norm. 

In view of the infinite dimensionality of Sx, it is possible to construct eigenf unctions with 

(2) 

very special properties. First, we consider the orthogonal projection onto Sj^ , which is given 
by the Bessel kernel 

(14) Exf{x)= f J^{\\x-y\)f{y)dy. 

In the notation of |HoI-IVj (vol III, Chapter XVII), the spectral projection for Ajjn for the 
spectral interval [0, A^] is given by 

(15) eo(x,A2) = (27r)-" / e'<"'«>rfe 

J\i\<^ 

Hence, the spectral projection onto £x is given by 

(16) -^eo(x,A2) = (27r)-" / e*<-'«>rf^, 

J\i,\=\ 

where S is the standard surface measure. 

If we fix y we obtain a normalized Euclidean coherent state (p\ = Jn-2 (A|x — y\). By the 

standard asymptotics of Bessel's function, ~ {r{-,y)\)~^^'^ as r{-,y)\ oo. Thus, (f^^ 
peaks at X = y and decays at a rate r{x, y)'^^"^ away from the peak point. Also, in the high 
frequency limit it decays like A"^/^ for fixed x ^ y. 



2.2. Flat tori. A fiat torus is a compact quotients M"/L where L is lattice such as Z". The 
Laplacian A of the fiat metric again commutes with the n vector fields In the compact 

case, an orthonormal basis of joint eigenfunctions (px is provided by the exponentials e*^'^'^^ 
where A G L*, the dual lattice. The corresponding Laplace eigenvalue is |Ap. Usually (for 
instance when studying nodal sets), we prefer the real orthonormal basis sin(A, x), cos(A, x). 
A key feature of these eigenfunctions is that they are linear combinations of a finite number 
(two) functions of the form a(a;)e*'^'-^^/'^ with h = |A|~^. Such functions are known as WKB 
modes or Lagrangian states. The heuristic scaling fH6|) is exactly true, with 



i.e. with dT^° = e^''^'^°'> 5 j^. The 'phase function' S{x) = (x, -j^) is the generating function 
of a Lagrangian submanifold A = {(x, ^ = j^)}, which is a Lagrangian torus in T^(R"/L). 
The eigenfunctions have norm one, |e*^^''^^| = 1. We also note that 



{Op{a)ex,ex) = / a{x,X)dx= / ad/Ji- 
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2.3. Standard Sphere. Eigenfunctions of the Laplacian Asn on the standard sphere 
are restrictions of harmonic homogeneous polynomials on M""*"^. 

Let Ajgn+i = —(4-^ + ■■■ + — ) denote the Euclidean Laplacian. In polar coordi- 

nates {r,u) on M"+^ we have Ai-n+i = - (^|^ + + ^'^s^- A polynomial P(x) = 

P{xi, . . . , Xn+i) on MJ^~^^ is called: 

• homogeneous of degree k if P[rx) = r^P[x). We denote the space of such polynomials 
hj Vk- A basis is given by the monomials 

= x"^ ■ ■ ■ , |q;| = ai + ■ ■ ■ + a^+i = k. 

• Harmonic if A^n+iP[x) = 0. We denote the space of harmonic homogeneous polyno- 
mials of degree k by Hk- 

Suppose that P{x) is a homogeneous harmonic polynomial of degree k on M'^^^. Then, 
= Ak^+iP = -{-g, + ^fKPM + ^AcjnPM 

=^ A5nP(cj) = {k{k - 1) + nk)P{uj). 

Thus, if we restrict P{x) to the unit sphere S*" we obtain an eigenfunction of eigenvalue 
k{n + k-l). Let Hk C ^^(S") denote the space of spherical harmonics of degree k. Then: 

• L'^{S"') = '^fc- "^^^ orthogonal. 

• Sp{Asn) = {XI = k{n + k-l)}. 

• dim?-^fc is given by 

+ k — l\ /n + A; — 3^ 



k J \ k-2 

The Laplacian A^n is quantum integrable. For simplicity, we restrict to S"^. Then the 
group SO (2) C 5*0(3) of rotations around the X3-axis commutes with the Laplacian. We 
denote its infinitesimal generator by L3 = r^. The standard basis of spherical harmonics is 
given by the joint eigenfunctions (|m| < k) 

As.Yj^ = kik + 1)YJ^; 

_d_Yk _ jyiY^ 
idd m '"'-'m- 

Two basic spherical harmonics are: 

• The highest weight spherical harmonic Y^. As a homogeneous polynomial it is given 
up to a normalizing constant by {xi + 1x2)^ in with coordinates (xi, 0:2, x^)- It is 
a 'Gaussian beam' along the equator {xs = 0}, and is also a quasi-mode associated 
to this stable elliptic orbit. These general notions will be discussed in ^ |TOl 

• The zonal spherical harmonic Yq . It may be expressed in terms of the orthogonal 
projection 11^ : L'^{S'^) — > Hk- 

We now explain the last statement: For any n, the kernel Ilk{x, y) of 11^ is defined by 

Ukf{x)= [ Uk{x,y)f{y)dS{y), 
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where dS is the standard surface measure. If {Yj^} is an orthonormal basis of Tik then 



m=l 



Thus for each y, Ilk{x, y) G Tik- We can L? normahze this function by dividing by the square 
root of 

I|nfc(-,i/)||i2 = / llk{x,y)Ylk{y,x)dS{x) = Ylk{,y,y)■ 
J s-^ 

We note that Ilkiy^y) = Ck since it is rotationally invariant and 0{n + 1) acts transi- 
tively on S"'. Its integral is dimTCk, hence, Ilk{y,y) = voi\s") dimT^fc. Hence the normalized 
projection kernel with 'peak' at yo is 

^ Ukix,yo)^/VdiS^) 
° y/dimHk 

Here, we put yo equal to the north pole (0, ■ ■ ■ , 1). The resulting function is called a zonal 
spherical harmonic since it is invariant under the group 0{n + 1) of rotations fixing pq. 

One can rotate Yq{x) to Y^lg ■ x) with g G 0{n + 1) to place the 'pole' or 'peak point' at 
any point in 5*^. 



2.4. Surface of revolution. By a surface of revolution is meant a surface, necessarily 
M = S*^ or M = M^/Z^, whose metric g is invariant under an action by isometrics. In 
this case, one can separate variables to analyze eigenfunctions, i.e. the joint eigenfunctions 
of Ag and ^ have the form e^"'^ipn,j{r). But this is not necessarily the best way to analyze 
eigenfunctions. In §12.2.2[ we will discuss how to express A in terms of action operators 
and how to obtain results on eigenfunctions that would be awkward if one used separation 
of variables. In the case of simple surfaces of revolution, the existence of a 'global Birkhoff 
normal form' for gives much more control than separating variables. 

2.5. H". Much of the discussion of §2.1l has an analogue on hyperbolic space. For the sake of 
simplicity we assume n = 2, so that the group of isometrics is SL{2, M). Then the analogue 
of Euclidean plane waves are the horocyclic (or hyperbolic) plane waves e'-*'^"'"^^^^'^^ where 
{z, b) is the function on x B (with B = the ideal boundary of H^) equal to the signed 
distance from to the horocycle passing through z and b. Equivalently, 

e^'''^ = 7^S = ^D(^,&), 
\z — op 

where Puiz, b) is the Poisson kernel of the unit disc. (We caution again that e^^'''^ is written 
g2(2,f)> yL2\ ). Helgason defines a non- Euclidean Fourier transform by means of these 

plane waves, and there exist analogues of the objects and results of §2.1[ Further one can 
define a covariant calculus of pseudo-differential operators by using the non- Euclidean Fourier 
transform by replacing ([9]) by 



(17) 



Op(a)e(^^+i)<^'^> = a(z,6, A)e(^^+i)<^'^\ 
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The most interesting aspect of lies in its quotients by discrete groups T C S'L(2,M). 
Eigenfunctions on the quotient H^\r are the same as automorphic eigenfunctions on sat- 
isfying ip{'~fz) = ipiz). Helgason has introduced a generahzed Poisson formula for eigenfunc- 
tions of exponential growth in the sense that there exists C > such that |v5(2)| < Ce'^'^°^^'^^ 
for all z. We assume Aip = X^ip and following a traditional notation put = ^ + r"^. Then 
in ([H], Theorems 4.3 and 4.29, it is proved that there exists a distribution Ti^^^^^ G V{B) 
such that 

(18) v^r{z) = [ e^'^+'^^^'^'^T„,^Jdb), 

Jb 

for all 2; G D. The distribution is unique if | + 7^ 0, —1, —2, ■ ■ • . 

The distribution Tir^^p.^ is called the boundary value of ipir and may be obtained from ipir 
in several explicit ways. One is to expand the eigenfunction into the "Fourier series" , 

(19) fir{z) = y^^ttn^rA^)^ 

in the disc model in terms of the generalized spherical functions defined by ([H], Theorem 
4.16) 

(20) e^^+^r){z,b) ^ J2 <l>r,n{^)b'\ b G B. 

Then (cf. H], p. 113) 

(21) T,r,^Adb)=Y,('nb''\db\. 

A second way is that, at least when Re (ir) > 0, the boundary value is given by the limit 
([H], Theorem 4.27) 

d{0,z)—>oo 

where c is the Harish-Chandra c-function and d{0, z) is the hyperbolic distance. 

In particular, eigenfunctions on the quotient H^\r are generalized Poisson transforms of 
boundary values. We fix an orthonormal basis {(firj} of eigenfunctions on H^\r and denote 
their boundary values by Tj^^, . As observed in |Z2] , when ipir^ is a F-invariant eigenfunction, 
the boundary values Tirj (db) have the following invariance property: 

ip„, M = {z) =^ e(5+-.)(7-'7f>>T,,^ ^d-fb) = e(5+-.)(-.^)Ti,^. (db) 

(22) 

This follows from the uniqueness of the Helgason representation and by the identity {'jz, 'jb) = 
{z,b) + (70,76). Conversely, any distribution with such a F- invariance property defines an 
eigenfunction on H^\F. It follows that the study of eigenfunctions in this setting is equiv- 
alent to the study of such boundary values. It is proved in |GUj that the boundary values 
are derivatives of C2 functions in the co-compact case. Similar results are proved in [MSJ 
with graphs of the boundary values. The boundary values will be discussed again in §15.81 
in relation to quantum chaos and in ^ IT71 in relation to analytic continuation. 
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2.6. The Euclidean unit disc D. Although our emphasis is on manifolds without bound- 
ary, we review the standard orthonormal basis of eigenfunctions for the disc. 

The standard orthonormal basis of real valued Neumann eigefunctions is given in polar 
coordinates by ipjn^n{r,0) = Cm,nSmm9Jmij'm,n'r'), (I'esp. Cm,n cos m9Jm{j'm,n'^)) where „ 
is the nth critical point of the Bessel function and where Cm,n is the normalizing constant. 
The A-eigenvalue is = {j'mn)'^- The parameter m is referred to as the angular momen- 
tum. Dirichlet eigenfunctions have a similar form with j'^ ^ replaced by the nth zero jm,n of 
Jm- Nodal loops correspond to zeros of the radial factor while open nodal lines correspond 
to zeros of the angular factor. 

If we fix m and let Xm,n oo we obtain a sequence of eigenfunctions of bounded angular 
momentum but high energy. At the opposite extreme are the whispering gallery modes which 
concentrate along the boundary. These are eigenfunctions of maximal angular momentum 
(with given energy), and ~ rn. As discussed in |B.B] . they are asymptotically given by 
the real and imaginary parts of e^^"''^ Aip{p~^^^ xU^y) . Here, Aip{y) := Ai{—tp + y) where Ai 
is the Airy function and {— tp} are its negative zeros. Also, s is arc-length along dD, p is a 
normalizing constant and y = 1 — r. 

2.7. An ellipse. Ellipses have several special sequences of eigenfunctions. One is a se- 
quence of eigenfunctions concentrating in a Gaussian fashion along the minor axis. Such 
eigenfunctions are known as Gaussian beams, or as bouncing ball modes, and quite some 
effort has gone into their generalizations to more general domains and to manifolds without 
boundary. Hence, we briefly review the existence of exact eigenfunction of this kind. We 
will discuss the construction of approximate eigenfunctions of such Gaussian beams in ^ ITOl 

One can separate variables in special coordinates on the ellipse and obtain exact, if rather 
esoteric, formulae for eigenfunctions of the Dirichlet or Neumann problem. We express an 
ellipse in the form + = 1, < a < 1, with foci at {x,y) = (±a, 0). We define 
elliptical coordinates {(p,p) by {x,y) = (a cos coshp, a sin sinhp). Here, < p < pmax = 
cosh~^ a"^, < (p < 2n. The lines p = const are confocal ellipses and the lines cp = const 
are confocal hyperbolae. The foci occur at if = 0,n while the origin occurs at p = 0, = f • 
The eigenvalue problem separates into a pair of Mathieu equations, 

{9^Grn,n C COS ipG^^n \^ n^fn,n 

2 
dpF,fYi,n C cosh pFjYi^n ^m,n'^n^tn 

where c is a certain parameter. The eigenfunctions have the form \E'm,n(</', p) = Cm,nFm,n{p) ■ 
GmAv^) where, Fm,n{p) = Ccmip, ^) and GmAv^) = ccmif, ^) (and their sin analogues). 
Here, ccm, Ccm are special Mathieu functions (cf. [C] (3.10)-(3.2)). The Neumann or Dirich- 
let boundary conditions determine the eigenvalue parameters knC. 

There exists a special sequence of eigenfunctions which are like Gaussian beams along the 
minor axis. For this special sequence, G = Gm,n are asymptotic to ground state Hermite 
functions. More precisely, 

(24) Gm,n{^; Xm,n) = c^,n(A^,.)e-"-" '^(1 + 0{X;^]J), 
while 

(25) F^,„(p; Xm,n) = e*^'"-/oVcosh2-+i'i-a+(p; A^,„) + e-^-"/oVcosh2-+w-a_(p; A^,„) 
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where a±(p; A ni,n) ~ ^^j=o ^i,jip^^rn,n '^^^ determined by the Dirichlet or Neumann boundary 
conditions. Moreover, from the L^-normahzation condition jj \'^m,n{p, f )P'^P = 1 it follows 

that Cm,n{^m,n) ~ ^rn,n- 

From and (l25l) . the Gaussian beams are roughly asymptotic to superpositions of 
g±ifcsg-Am,ny |B.Bj ). where s denotes arc- length along the bouncing ball orbit and y 

denotes the Fermi normal coordinate. It follows that outside a tube of any given radius 
e > 0, the Gaussian beam decays on the order 0(e^'^™'"'^ ). 

Before leaving this example, we should point out two further interesting sequences of 
eigenfunctions. As in the case of the disc, there exists a 'whispering gallery' sequence which 
concentrates on the boundary of the domain. As with Gaussian beams, there exists a gen- 
eralization of this sequence to any convex smooth domain in the form of 'quasi-modes'. A 
second interesting sequence becomes highly enhanced at the two foci. 

3. Local behavior of eigenfunctions 

By the local behavior of eigenfunctions, we often mean methods and results which pertain 
to all solutions of Am = X'^u on a ball B{xo,R), not just to solutions which extend to 
global eigenfunctions on (M, g) . More generally, we consider results which are obtained from 
a small ball analysis and which use covering arguments to draw global conclusions from 
local arguments. But we emphasize that some aspects of eigenfunctions discussed in this 
section are truly global, for instance the estimate of the frequency function in terms of the 
eigenvalue. Thus, we do not aim to segregate local from global results in this section, but 
do try to indicate when a result assumes that the eigenfunction is global. 

Local properties include: 

• Frequency function estimates; we also include Carleman estimates, although they 
often require integration of global eigenfunctions over all of M; 

• Vanishing order estimates at points; 

• Doubling estimates; 

• Bernstein estimates; 

• Lower bounds on masses in small balls. 

• Local structure of nodal sets. 

For background and references on general elliptic estimates, see |GiTrt [HL] . 

3.1. Eigenfuntions and harmonic functions on a cone. One can easily convert eigen- 
functions !f\ to harmonic functions in a space of one higher dimension. Two closely related 
options are: 

• Form the cone R+ x M and consider the metric g = dr"^ + r'^g. Let ip\ = r^u where 

« = 1 (^^4X + {n-iy -{n- 1)) . 

Let A be the Laplacian on the cone. Then, 

Av?A = 0. 

• Form M+ x M and consider c'^Va- Then {d'f + A){e^^*ipx) = 0. 
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This approach was first taken in [GaLl ILin] and used further in |Kut ICM] , among other 
places. It is a useful approach when the frequency function is employed, since the latter 
has its best properties for harmonic functions. In |CMj it allowed for the use of Harnack 
inequalities on balls of M+ x M where ^p\{x) is positive. 

3.2. Frequency function. The frequency function N{a,r) of a function m is a local 
measure of its 'dee; polynomial like function in Br{a). More precisely, it controls the 

local growth rate of u. In the case of harmonic functions, it is given by 

/ N / N rD(a,r) 

where 

H(a,r) = / u^da, D(a,r) = / iV-uPrfx. 



'dBr{a) J Br (a) 

A well-written detailed treatment of the frequency function and its applications can be found 
in ^ IKu] , following the original treatments in |GaLt IGaL2t ILin] . 

To motivate the frequency function, let us calculate it in the special case of a global har- 
monic function u on M". Then u may be decomposed into a sum of homogeneous components 
of non-negative integral order: 



UN, UN{ruj) = ttNr^^Niuj). 



u = 

N=0 

Then at is a spherical harmonic of degree on S"""^, and we choose a^v so that its norm 

^n— 1 



is equal to one. From the fact that the LfN are orthogonal on S*" , one easily calculates that 

d 

2 d log r 



(27) N(r) = ^^=°^^l"^l " = ^ log V 







2^N=0 







1 2 2N 



N=0 



I.e. 



A H{r)\ D(r) 
(28) -rA^og-^)=2- 



dr \ r"^i / H(r) 

Here, one also uses that |VmP = jg^ ud^u, where is the unit normal. From the 
formula it is immediate that if u is homogeneous of degree N, i.e. has only one component 
ttNT^^N, then N{r) is the constant function N{r) = N equal to its degree. 

We note that N{r) is analogous to the average energy of a statistical ensemble with 
partition function Z{r) = Yl^=o kivP^^- We recall a partition function for the canonical 
ensemble at temperature /? is given by Z{P) = J e~^^duj{E) where uj{E) is the density of 
states. The average energy is then {E) = — ^logZ(/5). Monotonicity of the average energy 

follows from the fact its derivative {{E— {E)Y) = — ^ log Z(P) is the variance of the energy. 

Frequency functions may also be defined for eigenf unctions. At least two variations have 
been studied: (i) where the eigenf unctions are converted into harmonic functions on the 
cone M"^ x M as in §3.1j (ii) where a frequency function adapted to eigenfunctions on M is 
defined. 

We first consider method (i) in the case of an eigenfunction cpx on 5". The associated 
harmonic function on the cone is precisely the homogeneous harmonic polynomial on M"+^ 
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which restricted to S*" gives (p\. By the previous calculation, A^(0,r) = where A = 
A^(A^ + n — 1). We note that on the cone M""*"^, the ball of radius r has the form [0, r] x S*", 
i.e the frequency function is global on S"". On a general manifold, the analogous global 
calculation is cleanest if we define N{r) with -8^(0) everywhere replaced by [0,r] x M. If we 
'harmonize' an eigenfunction Lfx as in §3.H we obtain Nir) = a. 

The second method is to define a frequency function on balls of M itself. The generalization 
to eigenf unctions ([T]) is as follows (see |GaLl IGaL2l IKu] ) . Fix a point a E M and choose 
geodesic normal coordinates centered at a so that a = 0. Put 

CJijXiXj 



^(a;) - I 12 ' 
Ixr 



and put 



Dia,r) := {d'^^^ + AV^) dV, resp. ff(a,r) := £ 



(29) D{a,r):= I [ g^^ + X^, ] dV, lesp. H{a,r) := f^^pi 

'dB, 

By the divergence theorem, one has 



(30) D{a,r)= / ipx ^ 
Define the frequency function of (p\ by 

(31) N{a,r) := 



Hir) 



As in the case of harmonic functions, the main properties of the frequency function of an 
eigenfunction are a certain monotonicity in r in small balls of radius O(^) (see Theorem 13. II) 
and the fact that N{a, r) is commensurate with N{h, r) when a and b are close. 

Simple examples show that, despite its name, the frequency function measures local growth 
but not frequency of oscillations of eigenfunctions, and therefore is not necessarily comparable 
to A. For instance, the frequency function of sin rax in a ball or radius < is bounded. 

An example considered in |DF] are the global eigenfunctions e*^ sin ty on of A-eigenvalue 
— t"^ and frequency function of size s. One could let s, t — >• oo with — bounded 
and obtain a high frequency function but a low eigenvalue. However, as discussed in [L] (p. 
291), if one forms the harmonic function from a global eigenfunction as in §3.H then one has 
^(0,2) < CX where C depends only on the metric. This is a global estimate since a ball 
centered at in the cone will cover all of M. An application is given in Theorem 13.21 

Let us work out the frequency function for a global eigenfunction on M", parallel to the 
discussion above for harmonic functions. We use the notation of §2.11 Since the tempered 
solution is unique up to constant multiples, we have Sx{N) ~ T-Cn{S"^^), as with harmonic 
functions. In the notation §2.1l of may therefore write 



<fx{ruj) = ^ a7vJAr,„(Ar)v9Ar(u;), 

Ar=0 

where as before | I^^atI |/,2((5n-i) = 1. Then, 

CO oo 

rD{0,r) = Xr^J2M''-JN,nJNA^r), H{0,r) = r"-^ ^ M^JN^nfiXr 



N=0 N=0 
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and 

Thus, the frequency function is again analogous to an average energy, but with respect to 
a partition function Z{r) = ^^=q |aArp( J^Ar,n)^(Ar) which no longer has the canonical form 
(i.e. it involves Bessel functions rather than powers). Unlike r", Bessel functions oscillate 
on the scale j and hence the frequency function is generally not monotone non-decreasing in 
r; but it is monotone for r in a small interval [0,ro(A)] with ro(A) = | for some a > 0. The 
following weaker monotonicity result is however good enough for applications to vanishing 
order and doubling estimates. 

Theorem 3.1. Theorem 2.3 of [GaL] (see also |GaL2l [Lml [H] and [Ki] (Th. 2.3, 2.4);. 
There exists C > such that e^'^{N{r) + A^ + 1) is a non- decreasing function of r in some 
interval [0, 'ro(A)] . 

Another basic fact is that the frequency of ip\ in B^{a) is comparable to its frequency in 
Bjiih) if a, h are close and r, R are close. More precisely, there exists Nq{R) << 1 such that 
if iV(0, 1) < Nq{R), then v^a does not vanish in Br, while if A^(0, 1) > Nq{R), then 

(32) iV(p,i(l-i?))<CiV(0,l), \JpeBR. 

3.3. Doubling estimate, vanishing order estimate and lower bound estimate. Dou- 
bling estimates and vanishing order estimates give quantitative versions of unique continua- 
tion theorems. Given a partial differential operator L and a solution u of Lu = 0, the unique 
continuation problems asks whether u is uniquely determined in a ball B by its values in a 
smaller set E C B7 That is, if m = in must u = in B7 In the limiting case where 
K = {xo} is a point, if u vanishes to high enough order k at xq must m = 0? A related 
question is the 'proximity to zero' of u in the sense of Nevanlinna theory, i.e. how small can 
sup^g^ \u{x) \ be in a ball B (see §3.41) ? These questions are answered by frequency function 
estimates and by Carleman estimates. 

An early doubling inequalities was proved by Bernstein for polynomials of one variable: 

(33) max |pAr(x)| < R^ max |pAr(x)| 

-R<x<R -l<x<l 

for any polynomial of degree N . A generalization known as Remez's inequality allows one 
to compare the growth on an interval to growth on any measurable set in the interior. Let 
cr G R+ and denote by Vn{<7) to be those polynomials pn of degree N such that UpatH^; < 1 
for some subset E C [— 1, 1 + cr]. Then 

sup < ||T7v||[-l,l+<x]- 

Here, \\f\\E = sup^; |/(a;)| and Tjv is Tchebychev's polynomial. 
In the case of harmonic functions on one may set 

- H{a,r) 
H{a, r) = — . 

The monotonicity of N{a, r) immediately implies the doubling formula 

(34) H{a, 2R) = H{a, R) exp (^j^"^ ^^("'^) ^^^ < 4^Ki-«)if(a, R). 
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Integrating in R and using ( l32i) gives 



^ ^ 1^0/(53 (a)) (a) - VolB.ia) Js.ia) 

As mentioned in §3.21 the frequency function of a global eigenfunction may be estimated 
in terms of the eigenvalue. 

Theorem 3.2. |DFl ILinj and [H] (Lemma 6.1.1) Let (fx be a global eigenfunction of a C°° 
(M, g) there exists C = C{M, g) and tq such that for < r < tq, 

Further, 

(36) max |o9A(a;)l < ( — ; ) max |o9A(a;)|, (0 < r' < r). 

B{p,t) \r' / x(iB{p,r') 

The doubling estimates imply the vanishing order estimates. Let a E M and suppose that 
u{a) = 0. By the vanishing order z/(m, a) of m at a is meant the largest positive integer such 
that D°'u{a) = for all |a| < u. The vanishing order of an eigenfunction at each zero is 
of course finite since eigenfunctions cannot vanish to infinite order without being identically 
zero. The following estimate is a quantitative version of this fact. 

Theorem 3.3. (see |DFj .- |Linj Proposition 1.2 and Corollary 1.4; and [H] Theorem 2.1.8.) 
Suppose that M is compact and of dimension n. Then there exist constants C{n),C2{n) 
depending only on the dimension such that the the vanishing order z/(m, a) of u at a E 
M satisfies a) < C{n) A^(0, 1) + C2{n) for all a e Bi/^^O). In the case of a global 
eigenfunction, v^ipx.a) < C{M,g)X. 

In the case of harmonic functions, one may write u = Pi, + ip^, where P^, is a homogeneous 
harmonic polynomial of degree u and where ipi, vanishes to order z/ + 1 at a. We note that 
highest weight spherical harmonics Cn{xi + 1x2)^ on S*^ are examples which vanish at the 
maximal order of vanishing at the poles xi = X2 = 0, 0:3 = ±1. 

3.4. Semi-classical Lacunas. For the purposes of this article, we define a 'semi-classical 
lacuna' to be an open subset U <Z M for which there exist a sequence {^x,^ } of L^-normalized 
eigenfunctions (M, g) and constants C, a > so that 

The sup-norm could be used in place of the norm. Another descriptive term is 'exponential 
trough'. 

Lacunae often arise as 'classically forbidden regions' of quantum mechanical systems, 
and in fact we do not know of any other examples. Consider for instance a semi-classical 
Schrodinger operators /i^A -|- V for which the classical energy level -|- V{x) = E projects 
over a compact subset Ke = {x : V{x) < E} C M. Then eigenfunctions of h'^A + V with 
eigenvalues Ej{h) E [E — 0{h),E + 0(h)] decay exponentially outside Ke at a rate given 
by 0{e~h^^^'^'^^) where d{x, Ke) is the distance from x to Ke in the Agmon metric, i.e. the 
metric {E — V{x))dx'^. 
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Classically forbidden regions also occur on compact Riemannian manifolds with integrable 
geodesic flow. For instance, on the round sphere, we consider sequences {Yj^} of joint 
eigenfunctions for which m/N ^ E. As will be discussed below, this sequence concentrates 
on the Lagrangian torus A^; C S*S'^ for which (|||, ^) = E. This torus projects to an 

5*^ invariant annulus Ke on 5*^. This annulus is the 'classically allowed region' and its 
complement is the classically forbidden region. It is not hard to show that |y^(a;)| < 
^-NdA(x,KE) ^]-^jg example, where cIa is a suitable Agmon distance. We will discuss this 
and related examples in more detail in §121 

To our knowledge, there are no converse results characterizing lacunae in terms of clas- 
sically allowed or forbidden regions. For instance, can lacunae occur if the geodesic flow of 
(M, g) is ergodic? Theorem 115.11 shows that lacunae cannot occur in the full density sequence 
of 'ergodic eigenfunctions', but might occur in a possible sparse subsequence. 

A more refined notion is that of microlocal lacunae, i.e. open subsets U C S*M such that 
{x{x, D)(px^ ,(p\ ) < CrXJ^ for all > and all xi^^D) G ^^(M) which are microsup- 
ported in U. In the real analytic case, one can insist on exponential decay, but in the general 
C°° case one could always add a smoothing operator to x(a;, D) to ruin exponential decay. 

One could also define lacunae depending on A. Namely, one could term a sublevel set of 
the form {x : Iv^aI^;)! < Ce~^^} a lacunae if its volume Vol{x : Iv^aI^;)! < Ce'"^^} is bounded 
below by some constant e > for some C, A. Eigenfunctions with a semi-classical lacuna 
clearly have this property, but there could exist other examples. It would be interesting to 
know if ergodic eigenfunctions can have lacunae of this form, or even if their sublevel sets 
can be larger than a j tube around the nodal set. 

The doubling estimates and Carleman estimates give quantitative lower bounds on the 
exponential decay rate of eigenfunctions in balls as the eigenvalue tends to infinity and show 
that the rate is never faster than in the definition above of semi-classical lacunae: 

Corollary 3.4. |DFj Suppose that M is compact and that ip\ is a global eigenf unction, 
Aifx = X^ifx- Then 

max \^px{x)\ > C'e"^^. 

As an illustration, Gaussian beams such as highest weight spherical harmonics decay at 
a rate e~'"'^^^^'"'^ away from a stable elliptic orbit 7. Hence if the closure of an open set is 
disjoint from 7, one has a uniform exponential decay rate which saturate the lower bounds. 
To our knowledge, it is unknown whether semi-classical lacunae can occur in more general 
situations than 'classically forbidden regions'. It even seems to be unknown whether semi- 
classical lacunae can occur on (M, g) with classically chaotic (i.e. highly ergodic) geodesic 
flows. 



3.5. Three ball inequalities and propagation of smallness. These inequalities gener- 
alize Hadamard's three circle theorem and Nevanlinna's two constants theorem. We only 
briefly mention the results and refer to [Ml IKu] for discussion of the work of Nadirashvili, 
Korevaar-Meyers and others and for further background. 

Let f2 be a domain in R" {n > 2), let G Q he a. non-empty subdomain, and let E G Q 
be a non-empty compact set. Define ||m|U = sup^g^ Then there exists a constant 

a = a{E,QQ,Q) G (0,1] such that \\u\\e < ll'^llnoll'^llfr" ^'^^ complex-valued harmonic 
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functions u on Q. If \\u\\q < C and if ||«||r2o << 1; the smallness of u on the subdomain VLq 
propagates to any compact subset E <zVt. 

In particular, one has a three spheres inequality (see (!29l) for notation) 

Theorem 3.5. (see |Kuj ) Lei < ri < r2 < ra. T/ien, 



-^^(^"2) < C*! J if(ri)°o+/3o if(r3)°o+/3o, 
where 

Q;o = log— , /5o = Clog— . 
r2 ri 

3.6. Bernstein inequalities. There are a variety of types of inequalities known as Bernstein 
inequalities. The original inequalities were proved for polynomials of one variable; a survey 
is given in |RY] . Further results in this direction (for analytic functions) are in |Brull [Bru2j . 
We will discuss these further in ^ IT71 

Among the classic Bernstein inequalities are the gradient estimates on polynomials of 
degree N: 

(1) |p'^(x)| < A^(l — a;^)~^''^||pAf||[-i,i] for x G [—1,1] (the Bernstein-Markov inequality); 
it implies that \p'^{x)\dx < 7rA^||p7v||[-i,i]. 

(2) \\p'n{x)\\[-i,i] < A^^||pjv||[-i,i] for x G [-1,1] (Markov's inequality); 

(3) \p'j,ix)\ < iV(l -a;2)-V2 (^\\p^\\2_^^^^ _ ^^^^y^ fo, ^ ^ 

The following result of Donnelly-Fefferman generalizes Bernstein's gradient inequality to 
eigenfunctions: 

Theorem 3.6. |DF3j Local eigenfunctions of a Riemannian manifold satisfy: 
(1) Bernstein estimate: 

(37) (! \V^,\Hv]'\ — ( ! W,\'dV^"' 



' B{p,t) / ^ \JB(p,r) 

(2) Bernstein estimate: There exists K > so that 

(38) max |V<y5A(a;)| < max |(y9A(a;)|- 

x&B{p,r) r xeB{p,r) 

(3) Dong's improved bound: 



IT7 i^c^Vx 

max Vv?A < max 09a 

Br{p) r Brip) 



for r < CaA-i/^ 



3.7. Carleman inequalities. Carleman inequalities are weighted integral inequalities which 
are an alternative to frequency function estimates in giving quantitative unique continua- 
tion results. We only indicate some results in this section and refer to the original articles 
|ArllDF|, lTa] and more expository articles |Ta2[ [H nXllEZj for further discussion. As discussed 
in [JLj . the idea of Carleman estimates is to use weights if) which are largest on the set from 
which one wants uniqueness (or other features) to propagate. On the other hand, there is 
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a constraint on ip in order that the Carleman estimate be true: it needs to be convex in a 
suitable sense. 

We first follow [E2l[lj. Let 

A{x, d) = ^ ajk{x)djdk + ^ ajdj + a{x). 

j,k j 

One searches for weights cp = e^"^ so that the weighted Carleman estimate holds: 



(39) 



[ e^^^ {t'\u\^ + r|Mp) dx<cJ [ e^^^\Au\^ + / e^^^ {t'\u\^ + t\Vu\^) dx) 
Jn \Jn Jan J 



for all T > Ci and all u G H(2){^). 

Theorem 3.7. [I] If ip is pseudo-convex with respect to A on Vt, then there exist constants 
Ci(A),C2 so that for \ > C2 and r > Ci{X), the Carleman estimate holds. 

Here, a function ip is called pseudo-convex with respect to A on f2 if 

dA — 
/l(x,O = 0; =0, for x G fi, ^ G M"\{0}, 



imply that 



sr^^ ^ dAdA ^ f.^ dA.dA ^ a^A \ ^ ^ ^ 



and if 

A{x,Vi){x)) ^0, Vx gH. 

(1) If A is elliptic then any ip G C^(fi) with Vip 7^ in f2 is pseudo-convex with respect 
to A on Vt. 

We are particularly interested in the Helmholtz operator A = A — A^, and wish to estimate 
the dependence of the constants on A. In this case one can put f{x) = r{x, b)"^. 

Theorem 3.8. [I], Theorem 3.1: Let (p{x) = \x - 6p. Then 

f e-'^^{T^\u\' + T\u\^)dx<c( [ e^^^\{A-e)u\^dx+ [ e^^"^ {t'\u\^ + T\Vu\'dS) 
Jn \Jn Jan 

for T > C and all real k and u G H(2) (^) ■ 

There is a microfocal interpretation of the convexity condition. Let P(x, D) be a pseudo- 
differential operator, and let h G [0, Hq) be a small parameter, and put 



(40) P^ = e^Pe-^. 
Then, 

P^ = P{x,D + ^V^). 

We denote the principal symbol of by p^. The Carleman weight condition (in the 
Hormander formulation) is 

(41) p^(x,O=0 =^ {Rep^,lmp^}(a;,0 > 0, 
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and in the boundary case, 

(42) dM^) ^ 0. 

||P^/||2 = ||Op(Rep^)/|p + ||Op(Imp^)/|p + (/, [Op(Rep^), Op(Imp^)]/). 
One has, 

[Op{Kepip),Op{lnip^p)] » modulo terms depending on Op{Rep^) , Op(lmp^) , 

and this is the condition {Rep^,lmp^}{x,^) > 0. 

In |DFl IDF2] (see also |JL] ). Donnelly- Fefferman use singular weight functions to prove 
vanishing order estimates and doubling estimates (Theorem 13. 4p . Their Carleman inequali- 
ties involve a function f which is comparable to the distance r. Let u G Cc{B{p, h)) and let 
r(x) = d{x,p). Put f(x) = Jq^^^ e~^^^ ds = r{x) + 0{r{x)^)) for some constant z/ >> 0. 

Theorem 3.9. [Xr] |DF3j (Lemma A) There exist constants C,B,5 such that, if a > (7(1 + 
A) and if u vanishes in a ball B{p,6) C B{p,h), then 

(43) / f-2°|(A-A')M|W>Sa' / f-2"-VrfK 

3.8. Geometric comparision inequalities. We now state a theorem due to P. Kroger 
and Bacry-Qian which compares eigenfunctions of (M, g) to eigenf unctions of model oper- 
ators L^ fn on an interval (a, 6). We follow the exposition in |BQ| . The model operators 
are essentially the radial parts of the Laplacian in geodesic polar coordinates on spaces of 
constant curvature . To motivate them, we note that the radial part equals 

• l^ + ^f on [0,oo] onM™; 

• J3 + ('^~1) cot on [0, vr] on the standard sphere 5*"^. A better comparison interval 
is to shift the pole to the center of the interval to obtain ^ — (m — 1) tanr^ on 

r 7r 7r] . 

[ 2 ' 2J' 

• ^ + ('^ ~ 1) cothr^ on [0, 00] on hyperbolic space TV^. 
We put 

Ji U /2 U Is = (0, 00) U (-00, 00) U (-00, 0). 
The model operators are defined as follows: 
(1) If > 0,m > 1, 



Lr.m 



on 1 , n A/ 1 |. 

I- 2 y m— 1' 2 y m— IJ 

(2) If i? < 0,m > 1 then 



(f)(x) = v"{x) — ^/ R{m — 1) tan v'{x) 

R TT I R 1 



^r,m(f)(a;) = t^"(x) + ^/^^Eirrl^^ coi\i (^y^^xj t;'(a;), on Ji U Jg 
-^^r-,m(^^)(x) = v"{x) + ^J-R{m - 1) tanh (^-^^xj t;'(x), on /2, 
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(3) If = 0,m > 1 then, 

Lo,mv{x) = v"ix) + ^ on Ji U Is 

Lo^mv{x) = v"{x) on I2, 

Theorem 3.10. Let {M,g) be a compact Riemannian manifold of dimension m with Ricci 
curvature hounded below by R. Let ipx be an eigenfunction of eigenvalue , let (a, 6) C M &e 
a finite interval, and let v\ be a solution of 

Lji^rnV = -X^v, on{a,b) 

satisfying: 

• v'{a) = v'{h) = 0. 

• v' ^ on (a, b). 

• [mmipx,ma.xipx] C [min ^>;^ , max w;^] . 
Then 

|VK-io<^,)l<i. 

In other words, consider an interval (a, b) so that is the lowest eigenvalue of the Neumann 
problem for Lpt^m on (a, b) and let vx be the corresponding ground state eigenfunction. Then 
if the range of (fx is contained in the range of Vx then at any points x G M and y G (a, b) 
such that ipx{x) = vx{y), \Vipx{x)\ < \Vvx{y)\- 

3.9. Symmetry of positive and negative sets. Let M be a smooth compact manifold 
and let (fx be a real nonconstant eigenfunction of the Laplacian on M. Let (f1^, resp. (f^ 
denote the positive and the negative part of (fx- Since /^_^ (fxdVg = 0, it is obvious that 
J^j (fxdVg = Jj^^ (f^dVg. This represents a symmetry in L^ between the positive and negative 
parts. On the other hand, eigenfunctions are not necessarily symmetric in this sense for 
higher norms. The article |JN2j of Jakobson-Nadriashvili discusses the extent of such 
symmetries and proves: 

Theorem 3.11. |JN2] Then for p ^ 1 there exists a positive constant C, depending only on 
p and M such that 1/C ^ IIv^a ||lp/||v^a lU'' ^ ^■ 

3.10. Alexandroff-Bakelman-Pucci-Cabre inequality. The following is Cabre's improve- 
ment of the Alexandroff-Bakelman-Pucci estimate. In the following, | |u| | denotes the norm 
in Q. 

Theorem 3.12. (see [Cab\ . Theorem I.4) There exists a constant C = C{M,g) independent 
of fl so that, for any subdomain fl C M, 

^ limsup \u\ + 2 1 |Am| |. 

Hence for any smooth function vanishing on d^, \\u\\i^oo(q^ < C\n\2\\Au\\. 

Corollary 3.13. If (fx satisfies Aux = A^m and(px\dn = 0, then \ \(px\\L°°(n) 

An interesting application of this inequality to extrema of eigenfunctions on nodal domains 
is given in [PSj . 
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3.11. Bers scaling near zeros. By Theorem 13.31 we know that at each xq G A/'^;^, the 
order of vanishing k of (fx is finite. 

Proposition 3.14. [Bersj IHW2] Assume that (px vanishes to order k at xq. Let (px{x) = 
V'A:°(^) + V^fc+i + ' ■ ■ denote the C°° Taylor expansion ofipx into homogeneous terms in normal 
coordinates x centered at xq. Then ipl°{x) is a Euclidean harmonic homogeneous polynomial 
of degree k . 

To prove this, one substitutes the homogeneous expansion into the equation Aipx = X'^'fx 
and rescales x Xx, i.e. one apphes the dilation operator 

II 

(44) Dl^(px{u) = (p{xo + j). 

The rescaled eigenfunction is an eigenfunction of the locally rescaled Laplacian 

A- := X-'DI^A,{DIT' = E ^ + • • • 

j=i j 

in Riemannian normal coordinates u at Xq but now with eigenvalue 1, 

A.Pa°)"V(^o + !) = AV(^o + f ) 

(45) 

=^ A^V(a^o + f)=^(xo + f). 

Since (p{xo + j) is, modulo lower order terms, an eigenfunction of a standard fiat Laplacian 
on M", it behaves near a zero as a sum of homogeneous Euclidean harmonic polynomials. 

The Bers scaling is used by S.Y. Cheng (see also earlier results of Hartman-Wintner 
|HWl IChll ICh2] ) to prove that at a singular point of (px in dimension two, the nodal line 
branches in k curves at xq with equal angles between the curves. For further applications, 
see |Bes] . Unfortunately, no analogue of the Bers scaling seems to exist at critical points 

3.12. Heuristic scaling at non-zero points. One may ask what becomes of the scaling of 
eigenfunctions on the length scale of 'Planck's constant' h = Xj^ around a point where the 
eigenfunction does not vanish. Since p{xq + j) is, modulo lower order terms, an eigenfunction 
of a standard flat Laplacian on M", it is reasonable to think that asymptotically (px behaves 
like a sum of plane waves, i.e. that there exists a distribution dTx on the unit sphere in 
momentum space such that 

(46) ¥^a(xo + y)~ / e^<«'">rfT,^°(0- 

This is not a rigorous definition of dTx since ~ has not been precisely defined. 

However, examples suggest that there do exist rigorous plane wave approximations. For 
instance, if we rescale the Helgason representation (fT8|) for an eigenfunction on a hyperbolic 
quotient, we obtain 

^A(.o + ^) = /^e(^^-^^)<-+i'^)rfr,. 

We have {zq + j,b) = {zq, 6) + | • 6 + 0(A^^) where u ■ b denotes the Euchdean inner product 
of u with the element of represented hy b E B = S^. Hence, 6^'^^^^^^^°''''^ dTx is a plausible 
candidate for the distribution. In this heuristic discussion, we neglect the remainder estimate 
which requires an estimate of the A dependence of dTx- 
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4. Nodal sets on C°° Riemannian manifolds 

In this section we review results on the nodal, critical and singular sets of eigenfunctions 
in the smooth case. In ^IT71 we review the much stronger results in the real analytic case. 
The nodal set of an eigenfunction ^px is the zero set 

(47) Z^^ = {xeM: ^x{x) = 0}. 

The global structure of the nodal set is determined by integrals fdHJ^'^ of continuous 

functions (or characteristic functions of nice sets) over Z^p^. This seems very difficult, so 
we study first the local structure of the set, e.g. its Hausdorf dimension and local Hausdorf 
measure. In ||T7] we present some global results on the nodal set when the geodesic fiow is 
ergo die. 

The following theorem, due to Briining after an observation of R. Courant (see [Br])), is 
used to obtain lower bounds on volumes of nodal sets: 

Theorem 4.1. If {M^,g) is a C°° compact Riemannian manifold without boundary, then 
there exists C{M,g) > so that, in each ball of radius > j there exists a point where ipx 
vanishes. 

Proof. Fix xq^t and consider B{xQ,r). If ipx has no zeros in B{xo,r), then B{xo,r) C Dj-x 
must be contained in the interior of a nodal domain Dj-x of (px- Now = Xl{Dj.x) where 
Xl{Dj.x) is the smallest Dirichlet eigenvalue for the nodal domain. By domain monotonicity 
of the lowest Dirichlet eigenvalue (i.e. Ai(fi) decreases as fl increases), < \1{Dj.x) < 
\l{B{xo,r)). To complete the proof we show that Xl{B(xo,r)) < ^ where C depends only 
on the metric. This is proved by comparing Xl{B{xo,r)) for the metric g with the lowest 
Dirichlet Eigenvalue Xl{B{xo,cr); go) for the Euclidean ball B{xo,cr; go) centered at Xq of 
radius or with Euclidean metric equal to g with coefficients frozen at Xq; c is chosen so that 
B{xo,cr; go) C B{xo,r,g). Again by domain monotonicity, Xl{B{xo,r, g)) < Xl{B{xo,cr; g)) 

for c < 1. By comparing Rayleigh quotients ^'Ji^v^'' one easily sees that Xl{B{xo, cr; g)) < 

CXl{B{xo, CT] go)) for some C depending only on the metric. But by explicit calculation with 
Bessel functions, Xl{B(xo, or; g^)) < ^. Thus, A^ < 

□ 

For background we refer to [Chj . A nice variation on the proof is given in |CMj . where 
eigenfunctions are converted to harmonic functions as §3.11 and Harnack's inequality on 
positive harmonic functions is used to prove Theorem 14. 1[ Another use of positivity is given 
in [Hj to prove Xi{B{xo,r)) > A: Let Ur denote the ground state Dirichlet eigenfunction for 
B{xo,r). Then > on the interior of B{xo,r). If B{xo,r) C Dj-x then also ipx > in 
B{xQ,r). Hence the ratio ^ is smooth and non-negative, vanishes only on dB{xQ,r), and 
must have its maximum at a point y in the interior of B{xo, r). At this point (recalling that 
our A is minus the sum of squares), 

so at y, 

^ f Ur\ (fX'^Ur - UrAipx (X'i(B(xo,r))-X'^)ipxUr 



34 



STEVE ZELDITCH 



Since > 0, this is possible only if Xi(B(xo,r)) > A. 

We recall that the nodal set of an eigenfunction Lpx is its zero set. When zero is a regular 
value of (px the nodal set is a smooth hypersurface. This is a generic property of eigenfunc- 
tions mi . It is pointed out in [Baej that eigenf unctions can always be locally represented in 
the form 

in suitable coordinates {xi,x') near p, where <fx vanishes to order k at p, where Uj{x') vanishes 
to order k — j a.t x' = 0, and where v{x) 7^ in a ball around p. It follows that the nodal 
set is always countably n — 1 rectifiable when dim M = n. 

For a general C°°{M,g) of dimesion n, S.T. Yau |Yll IY2] has conjectured that 

(48) CM,gX < Vr-\Z^^) < CM,g\. 

Here, Ti^ is the k dimensional Haussdorf measure. The conjecture was proved in [DF] in the 
real analytic case, which will be discussed in detail in ^ IT7t see Theorem 117.11 In the C°° 
case, the lower bound was proved in dimension 2 by J. Briining [Br] and S.T. Yau. Because 
of this result, and lower bounds in some other calculable cases, Yau conjectured the same 
lower bound in all dimensions. It remains an open problem in dimensions > 3. 
Regarding upper bounds, in dimension two one has 

Theorem 4.2. [DFl, [Dong1 Suppose {M,g) is C°° and that dimM = 2. Then, 

Dong's proof was based on the following integral formula: 
Theorem 4.3. |Dong| Let q = |Vup + ^m^. Then 



1 f {A + \^)\u\ 



K'-'iZu) = - / ^ ^ '-^dV. 

2 Jm Q 

In higher dimensions, the best estimate to date is the following: 
Theorem 4.4. [HS] For any {M,g) of dimension n one has 

4.1. Courant and Pleijel bounds on nodal domains. Another well-known result is 
Courant's nodal domain theorem: 

Theorem 4.5. The number rik of nodal domains of the k satisfies Uk < k. 

Here, a nodal domain is a component of M\Z^^^ . This estimate is not sharp. Pleijel [Pj 
proved: 

Theorem 4.6. For any plane domain with Dirichlet boundary conditions, limsupfc^^o^^ ^ < 
~ 0.691..., where ji is the first zero of the Jq Bess el function. 

ii 

He conjectured that the same result should be true for a free membrane, i.e. for Neumann 
boundary conditions. This was recently proved in the real analytic case by I. Polterovich 
[Poj . using a result |TZ3] counting the number of nodal lines which touch the boundary. 
Another recent result is the proof that there exists an asymptotic mean number of nodal 
domains of random spherical harmonics |NSj : see §18. 6[ 
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4.2. Critical and singular sets of eigenfunctions on Riemannian manifolds. 

Critical points are points where V'^\{x) = 0. We denote the critical set by S = Vv9^^(0). 
Singular points are critical nodal points, i.e. critical points lying on the nodal hypersurface, 
Sq = S n Af. The reader is warned that the term 'critical point' is sometimes used for 
singular points. In general, the singular set is simpler to study than the critical set. 

There exist simple examples such as surfaces of revolution where the critical point sets of 
eigenfunctions are of codimension one. For instance, the rotationally invariant eigenfunctions 
on a surface of revolution have invariant critical point sets, and not all consist of fixed 
points. Thus, zonal spherical harmonics have critical point sets consisting of unions of 
lattitude circles. 

For generic metrics, all of the eigenfunctions are Morse functions and consequently their 
critical point sets are discrete [U]. However, there are no known bounds on the number of 
critical points. In the analogy of eigenfunctions of eigenvalue to polynomials of degree 
A, a very naive application of Bezout's theorem suggests that the number of critical points 
should be bounded above by CA*^, since the critical point equation is a system of n equations 
of degree A — 1 in n unknowns. Even for real analytic metrics, no rigorous results in this 
direction are known. Any such bounds would have to refiect the non-degeneracy of the 
critical points; a simple count could be unstable if a small perturbation of the metric had 
a sequence of eigenfunctions with codimension one critical point sets. On the other hand, 
no lower bound on the number of critical points exists: it is proved in [JN] that there exists 
a Riemannian surface possessing a sequence of eigenfunctions with a fixed finite number of 
critical points, answering a question of S.T. Yau |Y3j . which asks if the number of critical 
point for the eigenfunction grows when the eigenvalue grows (Yau believes that the answer 
is positive for most metrics, and that it is interesting to understand it and to give some 
estimate of the growth.) In [Y3j . Yau also proved the existence of at least one nontrivial 
critical point for surfaces that can not be guaranteed by standard Morse theory. Up to now, 
this is the only result of proving existence of nontrivial critical points. 

Singular sets have been studied for harmonic functions, eigenfunctions, and more general 
solutions of elliptic equations. The singular set is somewhat simpler than the critical point 
set. In the real analytic case, the singular set is known to be at most of codimension 2 and 
n"'~'^{T.) < oo (see ^e^, 3.4.8). This was extended to metrics in jHHL] . One has: 

Theorem 4.7. [HNOOj IHHLj Let dimM = n. There exists a constant C^^ depending only 
on the frequency function of ip in -B/j(0) such that the singular set So satisfies: 

7Y"~^(So n Bji/2) < Cmh- 

A key point in the proof is that the complex critical set of a homogeneous harmonic 
polynomial in dimension n is of dimension <n — 2. 
Lin has conjectured: 

Conjecture 4.8. //dimM = n, the singular set S satisfies: 



Theorem 4.9. Dong] Suppose that dimM = 2 and let m{Xj) denote the number of singular 
points of (fx counted with multiplicity (i.e. m{Xj) equals the sum over all singular points of 



36 STEVE ZELDITCH 

the order of vanishing minus one). Then 

m{\^) < ^ (\'^Vol{M,g) -2 [ mm{K,0)dVg 



M 



where K is the Gaussian curvature. 



5. The wave kernel of a compact Riemannian manifold 

Global properties of eigenfunctions often arise from the fact that global eigenfunction are 
eigenfunction of the wave group Ut = e**^. We begin by reviewing some basic theory of the 
wave equation on a compact Riemannian manifold. 

The wave group of a Riemannian manifold is the unitary group Ut = e**^ is defined by 
the spectral theorem, 

(49) U{t,x,y) = Y,e''^^Vj{x)^M- 

j 

Closely related but simpler wave kernels are the even part of the wave kernel, cos ty/A which 
solves the initial value problem 



a 



-A)u = 



Similar, the odd part of the wave kernel, is the operator solving 



(51) 



u\t=o = fMt=o = 9 

To employ wave kernels in spectral geometry it is indispensible to have approximations 
defined in terms of geometric data. They were first constructed by Hadamard and Riesz and 
are usually referred to as Hadamard(-Riesz) parametrices. There are alternative parametrices 
due to Lax and Hormander. 

We begin with the small-time Hadamard Riesz parametrices for cos t^, sii^. These 
kernels only involve A and their kernels can be constructed in the form 



r°° °° n-l ■ 

(52) / e*^('''-*') Vvrj(x,i/)^."^"'d^ mod C 



j=0 

where Wj are the Hadamard- Riesz coefficients determined inductively by the transport equa- 
tions 

e' Ti/ , awo _ n 
2e^o + — - U 

(53) 

'^^r{x,y){{^ + g)iy.+i + ^} = AyW,. 

Here, r = r{x,y) is the geodesic distance and is a regularization of 6*^ at 6* = 0; t"" 
is the distribution defined by = Re{t + zO)"" (see [Be]). We recall that {t + iO)"" = 
e~''^2p^ Jp°° e**^x"~^(ix and also that e^^^x"'~^dx has precisely the same singularity at 

t = as the sum EZo'^^'^'^^^Kk + f )""'• 
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The solutions are given by: 

Wo{x,y) = Q-^x,y) 

(54) 

Wj+i{x,y) = 6 2(x,y) s''Q{x,Xs)2A2Wj{x,Xs)ds 

where Xg is the geodesic from x to y parametrized proportionately to arc-length and where 
A2 operates in the second variable. 

Performing the integrals, one finds that 

(55) costVAix,y) ^ CM ^(-l)^t/;,(x, y) ] mod 

j=0 ^ [J 2 ^) 

where Co is a universal constant and where Wj = Co(^~^''^4:~^Wj{x,y),. Similarly 

(56) -—^(x,!/) ^ Cosgnit) ^(-l)%,(x,t/) ^ . >- mod C^ 

Here, ~ means that the difference of the two sides is a C°° function on M x M, or more 
precisely, that if one truncates the sum after a number of terms, the difference lies in 
C^{M X M). The formulae are only valid for times t < inj{M,g) and for this reason are 
called small-time parametrices. When the metric is real analytic, the series for the ampltiude 
converges for t sufficiently small and {x,y) sufficiently near the diagonal |Bej . 

To obtain truly global results on eigenfunctions, one actually needs large time paramet- 
rices. These are very complicated and in many respects are not understood. The simplest 
way to obtain one is to use the group property of U{t) = U(t/N)^ to determine the wave 
kernel for all times from the wave kernel at a small time. It shows that for fixed (x, t) the 
kernel U(t){x,y) is singular along the distance sphere St{x) of radius t centered at x, with 
singularities propagating along geodesies. 

By a similar but more complicated calculation (using the action of \/A on the oscillatory 
integral for '^"^y^ , one has 

/•CO °o 

(57) f/(t, x,y)= / e^^(^'(^'^)-*') J] W,{x, y)e'^-'d9 {t < inj(M, g)) 

-^0 k=0 

where f/o(a;,y) = Q^^{x,y) is the volume 1/2-density, where the higher coefficients are 
determined by transport equations, and where again 6^ is regularized at 0. 
An alternative parametrix has the form 

(58) U{t,x,y)= [ e^*l«l««e-^<«'""P«"'(^)M(t,x,i/,Of^^ 

where \C,\g^ is the metric norm function at x, and where A{t,x,y,^) is a polyhomogeneous 
amplitude of order 0. The expression exp~^{y) is again only defined in a sufficiently small 
neighborhood of the diagonal. For background, we refer to |HoI-IVt ID.G] . 
The existence of these parametrices is sufficient to prove: 
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Theorem 5.1. U{t,x,y) E P'(R x M x M) is a a Fourier intetgral operator of order — | 
associated to the canonical relation, 

r = {(t,r,x,e,l/,r/) : r + 1^1, = 0, = {y,v)} CT*{RxMx M). 

This means that U {t, x, y) can be locally written as a finite sum of oscillatory integrals 
/irm e^'^^^'^''^'^^ a{t, X, y , ^)d$, whose phases locally parametrize T in the following sense: De- 
fine the critical set along the fibers by 

= {{t,x,y,0 eRx M X M X R^'^ : d^if = 0} 

and define the immersion i^ : ^ C T* (R x M x M) by 

i^^t, a;, y, i) = (t, dt^, x, d^ip, y, -dyip). 

For instance, when 



one has 



Vit^ X, y, = -^l^lgy + ('^' exp^ (x)) 

= {{t, X, y, : expy(t^) = x}, 
so that X, y are linked by a geodesic segment of length t and 

^^ = {it,-\^U,x,^,g\x,0)}- 

5.0.1. Manifolds without conjugate points. The Riemannian manifolds with the simplest 
wave groups are those without conjugate points. A clear exposition is given in |Bej . We 
recall that a Riemannian manifold {M,g) is without conjugate points if the exponential 
maps expj, : T^M —>■ M have no singular points. In this case, the universal Riemannian 
cover M, the total space of covering map tt : {M,g) — > {M,g), is diffeomorphic to R" 
(n = dimM) and the exponential maps expx '■ T^M — ^ M are diffeomorphisms for all x. 
The distance function f is globally well defined and real analytic away from the diagonal 
A^oiM X M. 

On a manifold without conjugate points, the Hadamard and Lax-Hormander parametrices 
for the wave kernel ?7(t, x, ?/) on M x M are well-defined for all t. The wave kernel on the 
quotient M can be expressed in the form 

(59) U{t,x,y)=Y,U{t,x,^y) 

where F is the deck transformation group of M — > M and where we implicitly identify x, y 
with one of their lifts to M. The series converges when one takes for U the costv^A or ^is^^: 
by finite propagation speed, there are only a finite number of terms for each t namely terms 
where f(x, 7?/) < t. The unitary kernel e**^(t, x,j/) does not have finite propagation speed, 
but outside the light cone the kernel is smooth and the same finite number of terms of the 
sum determine the singularity completely. 
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6. Methods for global analysis 

The global analysis is based on properties of eigenfunctions which derive from the fact that 
they extend globally to M and hence solve the wave equation eigenvalue problem e^^^ip\ = 
e^^^ipx. The wave operator is very useful because it propagates singularities along geodesies. 
Studying singularities of various expressions involving eigenfunctions and eigenvalues gives 
further information that the local methods do not provide. 

To take full advantage of the wave group it is important to use pseudo-differential opera- 
tors, which transform well under conjugation by the wave group. We briefly review some of 
the key properties of pseudo-differential operators in this section. For a detailed treatment, 
we refer to [Hol-lVl lEZl [D5j| [Ugjj 



We denote by ^''^(M) the space of pseudo-differential operators of order m on M. The 
principal symbol of A is denoted ga- 

6.1. Egorov's theorem. Egorov's theorem for the wave group concerns the conjugations 

(60) at{A) := UtAU; , A G ^'"(M). 

Such a conjugation defines the quantum evolution of observables in the Heisenberg picture. 
Egorov's theorem is the following: 

Theorem 6.1. at defines an order-preserving automorphism of \E'*(M), i.e. at{A) G 
^"^(M) if Ae ^™(M), and that 

(61) (yu,Aut{x,0 = ^a($*(x,0) := Vt{aA), (x,0 e T*M\0, 
where Vt is the unitary operator 

6.2. Sharp Garding inequality. The Garding inequality addresses the question: To what 
extent is the quantization Op{a) of a positive symbol a positive pseudo-differential operator? 
Can one give lower bounds for expressions {Op{a)(pj, ipj) in terms of a? 

The answer depends on the precise definition of Op{a). There is a somewhat complicated 
quantization Op^{a) due to Friedrichs in the case of such that Op^{a) > if a > 0. But 
for a general definition of Op{a) one can only expect the following sharp Garding inequality: 

Theorem 6.2. For any f e C°°(M), we have 

(Op(a)/,/)>(infa)||/|p-q|A-i/Vir, 
where \ \f\\ is the LF'-norm of f . 

This immediately implies the 

Corollary 6.3. If a>0, then 

{Op{a)^x,Vx) > -C\-~\ 

6.3. Operator norm and symbol norm. Another natural question comparing the prop- 
erties of the classical observable a to its quantization Op{a) concerns traces and norms. For 
instance. 

Theorem 6.4. For any A G ^1^°, ||cr^||L°c = infi^ ||y4-|-fr|| where the infimum is taken over 
the set of compact operators K. 
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6.4. Quantum Limits (Microlocal defect measures). One of the principal problems 
in the global analysis of eigenfunctions is to determine the possible weak limits of the Wigner 
measures, or microlocal defect measures, which arise from matrix elements. 
To define them, we regard diagonal matrix elements as linear functionals on 

(62) pk{A) = {A^k.^k)- 

We observe that pk{I) = 1, that pk{A) > if A => and that 

(63) pk{UtAU;)= pu{A). 

Indeed, if A > then A = B*B for some i? G and we can move B* to the right side. 
Similarly ( 163|) is proved by moving Ut to the right side and using the fact that the eigenvalues 
of Ut are of modulus one. In quantum statistical mechanics, these properties are summarized 
by saying that pj is an invariant state on the algebra or more precisely, on its closure in 
the operator norm. An invariant state is the analogue in quantum statistical mechanics of 
an invariant probability measure. 

We denote by A^/ the convex set of invariant probability measures for the geodesic flow. 
Further, we say that a measure is time-reversal invariant if it is invariant under the anti- 
symplectic involution (x, ^) — (x, — ^) on T*M. We denote the time-reversal invariant 
elements of M.i by 

Proposition 6.5. Any weak limit of the sequence {pk} on is a time-reversal invariant, 
(yf* invariant probability measure on S*M, i.e. is an element of . 

Proof. For any compact operator K, {Kipj, ipj) 0. Hence, any limit of {A(pk, ^k) is equally 
a limit of {{A + K)ipki Vk)- By the norm estimate, the limit is bounded by infi^- HA + i^'H (the 
infimum taken over compact operators). Hence any weak limit is bounded by a constant 
times IIu^IIloo and is therefore continuous on C{S*M). It is a positive functional since each 
pj is and hence any limit is a probability measure. By Egorov's theorem and the invariance 
of the pfc, any limit of Pk{A) is a limit of Pk{Op{aA o $*)) and hence the limit measure 
is invariant. It is also time-reversal since the eigenfunctions are real-valued, i.e. complex 
conjugation invariant. □ 

Problem 6.6. Determine the set Q of 'quantum limits', i.e. weak* limit points of the 
sequence {$fc} of distributions on the classical phase space S*M, defined by 

ad^k ■■= {Op{a)ipx^,(fxJ 
where a G C'^{S*M). 

The set Q is independent of the definition of Op. The simplest examples are the expo- 
nentials on a fiat torus M™/Z™. By definition of pseudodifferential operator, ^e^'^*^'^''^'* = 
a(x, /s)e^'^*^'^'^'^ where a(x, k) is the complete symbol. Thus, 

(64) (^^^2ni{k,x)^^2ni{k,x)^ ^ f a{x,k)dx^f aA{x,^)dx. 

A subsequence e^'"'^'^^^'^'^ of eigenfunctions has a weak limit if and only if tends to 
a limit vector in the unit sphere in M". In this case, the associated weak* limit is 
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f^„^j^„aA{x,^o)dx, i.e. the delta- function on the invariant torus T^^ C S*M for g^, de- 
fined by the constant momentum condition ^ = ^o- The eigenfunctions are said to locahze 
on this invariant torus. Given ^q, we can always define a sequence kj so that — and 
thus, every invariant torus measure arises as a quantum limit. 

In general, there are many possible limit measures. The most important are: 

(1) Normalized Liouville measure. In fact, the functional uj of fl72]) is also a state on \E'° 
for the reason explained above. A subsequence {^Pj^} of eigenfunctions is considered 
diffuse if p^^, — * u. 

(2) A periodic orbit measure fi^ defined by n-yiA) = f CAds where is the length 
of 7. A sequence of eigenfunctions for which p^^ /i^ obviously concentrates (or 
strongly 'scars') on the closed geodesic. 

(3) A finite sum of periodic orbit measures. 

(4) A delta-function along an invariant Lagrangian manifold A C S*M. The associated 
eigenfunctions are viewed as localizing along A. 

(5) A more general measure which is singular with respect to dfi. 

All of these possibilities arise as (M, g) varies among Riemannian manifolds. Indeed, the 
standard sphere provides an extreme example: 

Theorem 6.7. [JZ] For the standard round sphere S*", Q = Mj. 

An important case is when p^^ w, i.e. the limit measure is Liouville measure. The 
corresponding eigenfunctions become uniformly distributed on the energy surface S*M and 
are sometimes called 'diffuse'. By testing against multplication operators, one gets 

' \<pkAx)rdVol ■ 



Vol{M) Je ' Vol{M) 

for any measurable set E whose boundary has measure zero. In the interpretation of 
\{pkj{x)\'^dVol as the probability density of finding a particle of energy at x, this says 
that the sequence of probabilities tends to uniform measure. However, pk^ — > u; is much 
stronger since it says that the eigenfunctions become uniformly distributed on S*M and not 
just on the configuration space M. For instance, on the fiat torus ]R"/Z", the standard ex- 
ponentials e^'^*^*^'^^ satisfy [e^'^^^'^'^^p = 1, and are thus uniformly distributed in configuration 
space. On the other hand, as seen above, in phase space they localize on invariant Lagrange 
tori in S*M. 

The fiat torus is a model of a completely integrable system, on both the classical and 
quantum levels. On the other hand, if the geodesic flow is ergodic one would expect the 
eigenfunctions to be diffuse in phase space. The statement that the all eigenfunctions are 
diffuse, i.e. Q = {u}, is known as quantum unique ergodicity. It will be discussed in ^ IT5l 

Off-diagonal matrix elements 

(65) pjk{A) = {A^x,,'^x^) 

are also important as transition amplitudes between states. They no longer define states 
since Pjk{,I) = 0, are no longer positive, and are no longer invariant. Indeed, pjkiJJtAU^) = 
Q^K^j-^k) p.^i^j^'^^ so they are eigenvectors of the automorphism at of fl60|) . A sequence of such 
matrix elements cannot have a weak limit unless the spectral gap \j — Xk tends to a limit 
r G M. In this case, by the same discussion as above, any weak limit of the functionals 
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Pjk will be a time-reversal invariant eigenmeasure of the geodesic flow which transforms by 
e*"^* under the action of g*^. Examples of such eigenmeasures are orbital Fourier coefficients 
r~ lo^ ^~^^^'^A{9'^{x,^))dt along a periodic orbit. Here r G We denote by Qr such 

eigenmeasures of the geodesic flow. Problem 1 has the following extension to off-diagonal 
elements: 

Problem 6.8. Determine the set Qr of 'quantum limits', i.e. weak* limit points of the 
sequence {pkj} on the classical phase space S*M . 

As will be discussed in §15. 3[ the asymptotics of off-diagonal elements depends on the 
weak mixing properties of the geodesic flow and not just its ergodicity. 

7. Singularities pre-trace formulae 

One of the principal methods for relating eigenfunctions and geodesic flow are the local 
Weyl laws. One form of the local Weyl law is to find the asymptotics and remainder for the 
pointwise sums, 

N{X-x)= \v,^{x)\\ 

j:\j<X 

Another is to find, for A G the asymptotics and remainder for 

In both cases, the asymptotics are determined by a Fourier Tauberian method, by studying 
the singularities of the dual trace, 

S{t,x)=J2e'''^^xM"^ 
j 

resp. 

j 

7.1. Duistermaat-Guillemin short time pre-trace formula. By a pre-trace formula 
one means a 'geometric' singularities formula for the wave kernels 

U{t,x,x) = Y^'''^xA^)' 
j 

on the diagonal. The actual trace formula is the integral over M of the pre-trace formula. 
The trace formula clearly gives only spectral information, while the pre-trace formula gives 
information on eigenfunctions. 

Proposition 7.1. |D.G] Let {M,g) be a C°° compact Riemannian manifold of dimension 
n. Then there exists a sequence u)i,uj2, ■ ■ ■ of real valued smooth densities on M such that, 
for every p G with supp p contained in a sufficiently small neighborhood ofO and p = 1 

in a small neighborhood of 0, 

oo 

j k=0 
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as A — > OO (and rapidly decaying as X ^ —oo) with 

u!o{x) — Vol{S*M), uji — — cun', i^k — for odd k. 
Proof. One uses a short-time paramctrix, 

U{t, X, y) = (27r)-" / e^'^(*'^'^''')a(i, x, y, rj)drj 

where a is a classical symbol of order m and where ip{t,x,y,ri) = ip{x,y,ri) — t\ri\, with 
ip{x, y,r]) — if {x — y,r)) — 0. Hence, 

p*dN{X,x) = / e'^*p{t)U{t,x,x)dt 

= (27r)-" / / e'^'p{t)e-''^"^a{t,x,x,ri) dri dt . 

One now changes variables rj Xrj, puts the d$, integral into polar coordinates ^ = rcu, |a;| = 1 
and carries out the dtdr integral by the method of stationary phase. 

□ 

The same kind of argument applies to Na_{X): 

PROPOSITION 7.2. For A e ^'"(M), let Na{X) = ^x.Kxi^fj^'Pj)- Then for any p G S{R) 
with p e C^(M); supppfl Lsp{M,g) — {0} and with p = 1 in some interval around 0, we 
have: 

oo 

p * dNA{X) ~ Yl (A ^ +oo) 

k=0 

where: 



n = dimM, cto = / o'Adp, ak = oukdp 

Js*M Js*M 

where ujk is determined from the k-jet of the complete symbol a of A. 

Proof. The only new step is to apply A to the parametrix for Ut- Applying a ipDO to 
ae^'^ produces an expression ae^"^ with the same phase and only a change in the amplitude. 
Hence, 

AU{t, X, y) = (27r)-" f a{t, x, y, r;)e'^(*'^*'')dr; 
where a is a classical symbol of order m. Now proceed as before. 

7.2. Long time pre-trace formulae. We now discuss pre-trace formulae for U {t, x, x) for 
long times t. The main difficulty is that there is no better parametrix for U{t, x, x) than to 
use the group formula U{t) = U{t/N)^ to reduce to times < inj{M, g). But this will require 
N xm integral signs and amplitudes that are difficult to control. Moreover, the singularities 
become very difficult to control. 

One can get a good idea of the difficulties by studying the simplest case, that of manifolds 
without conjugate points. This case is a rather straightforward generalization of the Selberg 
pre-trace formula for compact hyperbolic manifolds. The key feature is that there exists a 
global in time parametrix on the universal cover, so it is not necessary to assume that the 
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suppport of p is small. As in the classical trace formula, one organizes the elements 7 e F 
(the deck transformation group) into conjugacy classes 7. One then has 

U{t,x,x) = U^jt, X, x), 
7 

where 

Uxf{t,x,y)— U{t,x,a''^jax) — U{t,ax,^ax), 

since a is an isometry. Here, is the stabilizer in F of 7. 

In this case, it is best to use the Hadamard parametrix. Since we are interested in the 
long time singularities, we can use the phase function r — t instead of — i^. Either one 
parameterizes the graph of the geodesic flow away from r — 0. One then has, 

p*dN{X,x) = E-yer / e'^'p{t) U{t,x,'yx)dt 

Jr 

= (2vr)-"E-yEaer/r, / / e'''p{ty'^'^^''''^'''^-'^^a{ax,^ax,e) dO dt 

Jr Jr+ 

Jr Jr+ 

To determine the asymptotics as A cxd one again applies stationary phase. The phase 
function is critical when 9 = l,t = r{ax,^ax). This corresponds to the times t when there 
exists a geodesic loop at x G M (i.e. at ax e M). The geodesic loops are in one-one 
correspondence with conjugacy classes in F and hence form a countable set. However, the 
growth rate of this set as t 00 is often exponentially large. Thus, if supp p C [— T", T], then 
there are often e^^ terms in the sum over 7. This makes it difficult to control the remainder 
terms. 



7.3. Safarov trace formula. Safarov added some precision to the Duistermaat-Guillemin 
singularities pre-trace formula. For fixed x, Given x e M, we let denote the set of loop 
directions at x: 

(66) C^^i^e SIM : 3T : exp, = x). 

We let Tj. : S*M ^ U {00} denote the return time function to x, 

mi{t > : exp^ — x}, if ^ G jCx', 

+00, if no such t exists. 



We then define the first return map by = : S*M. We also define T^''\0 to be 

the time of kth return for directions which loop back at least k times. 

We then consider the positive partially unitary operator (the Pcrron-Frobenius operator) 



: L'{C,, \du;\) ^ L'{S:, \du;\), UJ{0 



0, e ^ U 
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Here, Jx is the Jacobian of the map i.e. = Jx{,i)\di\. We have: 

keiUx = {f ^L\Sl): supp /n = 0}; Im[/, = {/GL2(^:):supp/c£4. 

We further define 

Let pT be the dilated test function satisfying pr(T) = p(^)- The pre-trace formula then 
has the form, 

(67) pT*dN{\x) = ao{x)X'-' + y-' / /)(^y^)t/.(A)'=Me| + ot,.(A"-1). 

"^^^ fc=i 

A key point in the proof is that the phase of the oscillatory integral for the left side only 
has a stationary phase point at ^ e C^- That reduces the integral to one over modulo an 
error ot(A'^~^). For the details, we refer to Proposition 4.1.16 of [SV| . 

8. WEYL law AND LOCAL WEYL LAW 

The classical Weyl law asymptotically counts the number of eigenvalues less than A, 

(68) N{\) = #{j : A, < A} = -^Vol{M,g)\- + Oir-'). 

Here, \Bn\ is the Euclidean volume of the unit ball and Vol{M,g) is the volume of M with 
respect to the metric g. Equivalently, 

(69) T..B. = M|iLi^ + 0(A»-). 

where Vol is the symplectic volume measure relative to the natural symplectic form Yl]=i dxjA 

d^j on T*M. Thus, the dimension of the space where H = ^/A is < A is asymptotically the 
volume where its symbol \^\g < A. 

The Weyl law with remainder is proved by using the integrated version of Proposition 
17.11 together with the Fourier Tauberian estimate A^(A) — p * N{X) = 0(A"^^), valid when 
p*dN{X) = 0(A"^^). The latter is visibly true from Proposition 17. 1[ See the Appendix ^ IT9| 
Theorem 119.11 

Proposition [7]T] assumes that p is supported in a small interval around t = 0. Thus it only 
takes into account the generic singularity at t = 0. An improved, two-term Weyl law has 
been proved which takes into account the singularities of Tr cos ty/A for larger values of t. 
The singular t ^ are the lengths of the closed geodesies 7 of (7*. The size of the remainder 
reflects the measure of closed geodesies. 

(1) In the aperiodic case, Ivrii's two term Weyl law states 

N{X) = #{j : A, < A} = Vol{M,g) + o(A™-^) 

where m = dimM and where Cm is a universal constant. 

(2) In the periodic case, the spectrum of a/A is a union of eigenvalue clusters Cn of the 
form 

C^ = {(^)(A^+^)+/iiv., 1 = 1. ..d^} 
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with /iATj = 0(A^ ^). The number of eigenvalues in Cn is a polynomial of degree 
m — 1. 

We refer to [D.Gt IHoI-IVt ISV] [ZT] for background and further discussion. 
The integrated Weyl law is relevent to spectral theory, while the pointwise local Weyl law 
is relevant to eigenfunctions. By the same Fourier Tauberian theorem, one has ( jPTGj ) . 



(70) Y. = (^'^"1^" + 



A,<A 



where -R(A, x) = 0(A"~^) uniformly in x. In this case one obtains a o(A"~^) remainder if the 
set of geodesic loops at x has measure zero. Such refinements will be discussed in §14. 2[ 
The other local Weyl law concerns the traces TrAE{X) where A G \1/™(M). It asserts that 

(71) Yl ^^■) = WVn{l ^^dxdA A" + 0(A""i). 

A,<A ^^^^ ^-^B'^I ) 

When the periodic geodesies form a set of measure zero in S*M, one could average over 
the shorter interval [A, A + 1]. Combining the Weyl and local Weyl law, we find the surface 
average of a a is a limit of traces: 



(72) ^ 
Here, /i is the Liouville measure on S*M. 



9. Local and global estimates of eigenfunctions 

One of the applications of pointwise and local Weyl laws is to obtain bounds on norms 
of L^-normalized eigenfunctions, when p is sufficiently large. Estimates of the L°° norms can 
be obtained from the local Weyl law (112.41) . We use the notation /(A) = 0{g{X)) if there 
exists a constant C such that /(A) < Cg{X), and /(A) = Q{g{X)) if there exists a constant 
C such that /(Aj) > Cg{Xj) for some sequence Xj oo (in other words, the negation of 
/ = o{g{X)). 

Proposition 9.1. Let {M,g) be a compact m- dimensional C°° Riemannian manifold. Then 

m — 1 

IIv^aJU^^ = 0{Xj 2 ). 

Proof. Since the jump in the the left hand side of (170]) at A is J2j x^=x the jump 

in the right hand side is the jump of R{X, x), we have 

_ _ m — 1 

(73) Y \^xM' = OiXf-') =^ ||^,J|^. = 0(A/ ). 

j:Aj=A 

□ 
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We note that this estimate is stronger by the factor A^^ than the estimate obtained by 
Sobolev estimates of eigenfunctions: If dimM = m, then 



A \Li- 



Thus, when p = 00, the Sobolev estimate gives ||v2a||l°° < A™ ||v9a||l2. 

For general L^-norms, the following bounds hold on any compact Riemannian manifold 



Theorem 9.2. Sog 



(74) = 0(A^(^')), 2 < p ^ cx). 



where 

(75) 5{p) 



2(n+l) ^ ^ 



n—1 

These estimates are sharp for some {M,g). For instance, on the unit sphere 5*" C M""^^, 
the zonal spherical harmonics achieve the maximal L°° norm. Moreover, the highest weight 
spherical harmonics on S"^ saturate the bounds for small norms: using Gaussian integrals, 
one easily finds that \\{x + iy)''\\L2{s^) ~ k~^^'^. For instance, 

/jg3(x^ + y'^Y e^^'^^^y^^'"'^'' dxdydz = ||(a; + iy)''| 1^2(52) r'^^e'^'^r'^dr, 

^ ll(- + iy)'llV, = iJ|i~*-"^. 

Thus, the normalized highest weight vector has the form k^^^(x + iyY- It achieves its 
L°° norm at (1,0,0) where it has size k^^^. More importantly, it is an extremal for for 
2 < p < 6. For instance, 

/^3(x' + y')'^e-(^-'+2^'+^') = ||(x + «2/)'=||i6(52)/o°°r6V^V2c/r, 

^ \\(^^A,.\k\\Q _ r(6fc+l) _ , 7 -1/2 

Hence, the norm of k^/^[x + iyY equals 

Since A^ ~ A; and 5(6) = | in dimension 2, we see that it is an extremal. 

It is natural to ask if extremals (in order of magnitude) for the norms necessarily 
resemble these examples. This motivates the problems: 

(1) What is the structure of eigenfunctions which are maximal for the W norms? 

n — 1 

(2) Determine the {M,g) for which L°°{X,g) = 

(3) At the other extreme, determine the structure of eigenfunctions which minimize the 

norms, and determine the (M, g) with a uniformly bounded orthonormal basis of 
eigenfunctions. Irrational fiat tori are examples. 
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9.1. Sketch of proof of the Sogge estimate. We now sketch the proof of the Sogge 
estimates, following the appraoch in |KTZ J. 

The estimates are proved by interpolation from three estimates: p = oo, "^^^^^ (for n ^ 
1). The estimate is of course trivial. The estimate already followed from the un- 
integrated local Weyl law with remainder estimated by the Fourier Tauberian method. We 

2(ra+l) 

will nevertheless sketch an alternative proof which has more in common with the L "-i 
estimate. 



9.1.1. The estimate. 

Proof. The first observation is that the estimates are local in phase space T*M. That is, the 
estimate holds for (fx if it holds for X~^D)(px for any microfocal cutoff X^^D) = 
Opx-i{x), where x is a smooth function on T*M\0 which is supported in a small conic subset 
(i.e. the cone through a small open subset of S*M.) 

The semi-classical symbol of A~^A — 1 equals p{x, ^) = |^|^ — 1. We may choose coordinates 
in the support of x so that <9gi(|'C|g ~ 1) 7^ 0. Locally write p = p(x,^',^i), where C,' = 
{^2, ■ ■ ■ , Cm)- By the implicit function theorem, there exists a local function a{x, such that 
,^1 = a{x, C,') when p = 0, and there exists a symbol q{x,C,) > 0, on the support of x so that 

(76) p{x,0 = (l{x^0{Ci~(^{x^O)- 

Quantizing the right factor and using the support properties of the symbol, we get 
(77) 

(g(x, X-'D){X-^D,,) - a{x, X-^D')) (x(x, X-'D)^x) = {X'^A - l)x{x, X~'D)ipx 

= [A-2A,x(a;,A-iZ})]^A = 0(A-i). 

Since g(x, X~^D) is elliptic we may invert it to obtain, 

(78) /a(xi, x') := {X-^D.,, - a(x, X~^D'))x{x, X-'D)ipx = 0{X-'). 
Let E{t) be the unitary evolution operator which solves 

(79) (A-^ A - a{t, x', X-^D'))x{t, x' , X'^D)u = f 

with E{0) = Id. The estimate then follows from the Sobolev bound in dimension m — 1 
as long as 

(80) \\x{x,X~'D)ipx{xir)\\LHR^-^) = 0{l). 

This follows by solving the ODE (179!) in the Xi variable. If the right side were equal to 
zero, it would be solved by a unitary evolution operator E{t){x{x, D)ipx)fx{0, x'). One then 
solves the inhomogeneous equation by the Duhamel formula, 

{X{x, X-'D)^x){xi, x') = E{x^){x{x, X-'D)^x){0, x') + iX T E{t ~ s)0{X-')ds, 

Jo 

proving fl80l) and completing the proof of the Theorem. 

□ 
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2(n+l) 

9.1.2. The L "-i estimate. The estimate to be proved is that 

n-l 

(81) WfxW 2(n + l) < CA2{" + 1). 

L 

We start again with the symbol factorization (ffGll and the estimate (fTfll . which implies 

(82) [ \\f{xi,-)\\mR^-i)dxi < C||/||i2(Mn) = 0{X-'). 

Jr 

We now use the evolution operator E{t) of (!79|) . The equation is time dependent, so the 
solution operator solves an operator ODE of the form 

(83) X-^DtE{t,s) + A{t)E{t,s) = 0, E{t,t) = Id. 

Then there exists a parametrix for E{t) of the type (11401) . More generally, suppose that 
F(t, s) solves the equation (119. ip with initial condition F{s, s) = G{s){x, X~^D), a smoothing 
operator. Then for small times t G (0, to) there exists an oscillatory integral parametrix 

F{t, s)u{x) = X-'' [ e'^('^(*'^'"''?)-<J''^>)6(t, X, r/; X)dydri + R{t, s)u{x), 



where R{t, x) is a smoothing operator. 

Let U{t, s) = ip{t)F{t, s)x{x, X~^D) where x ^ C'^(]R). Then the main point is to prove. 



(84) snp\^JJ\Uit,s)f\\l,^^^„^dtj <CA-^||/|U2(Mn). 

2{n4-i) 2,(n-\-l) 

For the apphcation to the L "-2 estimate one takes p = q = , and we assume this 

from now on. The idea of the proof is to use the oscillatory integral parametrix for F and 
U{t, s) and the stationary phase method to prove, 

(85) \\U{t,r)U{s,r)*f\\Lo. < CA^^-^^/^^^-i ^ _ ^|)-(n-i)/2^ 

Then (!84l) follows by an abstract Strichartz estimate (see |KTZ] for references). 
One can choose ip and x so that 

,-1 



X{x, X-^D)^x{xi,x') = iX I U{xi, s)fx{s, x')ds + 0{X~°^). 

Then, 







n — 1 

Mx,X-^D)^x{xi,x')\\ < XX'^(-+^) j^\\h{x,-)\\L\m^-^)ds + 0{X-'^) 

T, n — l 

(86) 

n-l 

< CA2("+i). 

9.2. Generic non-sharpness of Sogge estimates. As the proof indicates, the Sogge 
bounds are '(micro-) local results', and do not take the global dynamics of the geodesic 
flow into account. The dynamics of the geodesic flow has a strong impact on the growth of 
1/ norms of eigenfunctions. In ^SoZj . it is shown that the Sogge U' estimates for large p are 
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very rarely sharp. To state the result, we need some notation. Let Vx := {ip : Aipx = Xfx} 
denote the A-eigenspace for A G Sp{A) and define 

(87) L^iX,g)= sup Ml^, £°°(A,^) = inf ( sup 

V^v,^ ONB{ipj]<E_Vx j=l,...,dimV';, 

l¥'ll^2=l 

1 

Thus, L'^i^g) = 0(A— ) for any {M,g). 

We define the measure \Cx\ of loops at x (see ( 166!) for the definition) as the surface measure 
on S*M induced by the metric gx on T*M. For instance, the poles xn,xs of a surface of 
revolution {S'^,g) satisfy l^^,! = 27i. 

^ n — 1 

Theorem 9.3. If L°°{X,g) = Q{X~), then there exists a point xq such that \Cxo\ > 0. In 
particular, if g is analytic, then all geodesies leaving xq must return to xq at the same time. 

The proof is based on a study of the remainder term in the local Weyl law. The main step 

is: 

Lemma 9.4. Let R{X,x) denote the remainder for the local Weyl law at x. Then 

(88) i?(A,x) = o,(A"-^) zf \Cx\=0. 

Additionally, if \Cx\ = then, given e > 0, there exists a neighborhood J\f of x and A < oo, 
both depending on e so that 

(89) |i?(A, y)\ ^ 1/ G AT, A ^ A. 

Thus, there are topological conditions on M which are necessary for M to possess a real 
analytic metric such that some sequence of eigenfunctions has the maximally growing L°° 
norms. In dimension two, only the sphere possesses such a metric. Moreover, the maximal 
growth rate of exhibited by zonal spherical harmonics can never occur for a metric on 
S"^ with ergodic geodesic fiow. 

Theorem 19.31 is certainly not the end of the story. In all known cases where 'maximal 
eigenfunction growth' occurs, all geodesies leaving the 'pole' Xq return to it at the same time 
T, not just a set of positive measure. Further, the 'first return map' : S'^^M S*^M is 
the identity map. One would hope to prove that the second property holds at least on a set 
of positive measure. 

When the geodesic fiow is 'chaotic' (i.e. highly mixing), it is expected that the eigen- 
functions resemble Gaussian random functions. The random wave model (see ^ IT8|) then 
predicts that eigenfunctions of Riemannian manifolds with chaotic geodesic fiow should have 
the bounds Hv^aHlp = 0(1) for p < oo and that ||v2a||l°° < Cv^log A. But no rigorous PDE 
methods to date have done better than 0{^^). There also exist counterexamples to the 
logarithmic estimate on special arithmetic hyperbolic quotients (see |RSt HSl lDon2j ). As in 
the case of zonal spherical harmonics on surfaces of revolution (but using Hecke operators 
in place of the rotations), there exist special eigenfunctions which take large values at 'fixed 
points' of the Hecke action or theta-corresopondence. In general, the exponential growth of 
the geodesic fiow is a huge obstacle to improving the estimate beyond the logarithm. Further 
discussion of L°°-norms, as well as zeros, will be given at the end of ^ITSlfor ergodic systems. 
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10. Gaussian beams and quasi-modes associated to stable closed geodesics 

We noted above that highest weight spherical harmonics maximized U' norms for 2 < p < 
'^^n-i^ , and we noted that they had the shape of an oscillating bump along the equator with 
Gaussian decay in the transverse direction. Such eigenfunctions are only known in a few 
cases, all of them quantum completely integrable: convex surfaces of revolution, ellipsoids 
and ellipses. What they have in common is the existence of a stable elliptic closed geodesic 
along which to construct a Gaussian beam. 

The construction of Gaussian beams is possible on any compact Riemannian manifold 
with a stable elliptic closed geodesic. However, in general it only produces an approximate 
eigenfunction or quasi-mode. We review them in this section since they are a good introduc- 
tion to quantum Birkhoff normal forms and local models. In the following section we will 
consider quantum integrable systems where the construction produces true eigenfunctions. 
References for Gaussian beams and more general quasi-modes are |B.B[ IRaj IRa2[ ITZ21 IW] 
among many other places. We follow [Z8] and |B.B] Chapter 9. 

Quasi-modes along stable elliptic orbits are approximate eigenfunctions (sometimes exact) 
which have the form of Gaussian beams along geodesics. A Gaussian beam is a simple 
oscillatory function e*'^'^ along the geodesic times a transverse Gaussian or (more generally) 
Hermite function. Hence, the quasi-modes are very localized along closed geodesics and in a 
sense are the best localized approximate eigenfunctions. 

Let (M, g) be a Riemannian manifold of dimension m = n+1 with a stable elliptic closed 
geodesic 7 of length L. We assume 7 is an embedded (non self-intersecting) curve with an 
orientable normal bundle, and use Fermi normal coordinates [s, y) along 7. That is, we fix 
an origin 7(0) and let s denote arc-length on 7. The exponential map exp : A^^ — > T^ij) 
from the e-ball in the normal bundle along 7 to a tubular neighborhood of radius e is a 
diffeomorphism and any choice of linear coordinates on endows T^{'y) with Fermi normal 
coordinates. 

We denote the eigenvalues of the linear Poincare map by {e*"-' }. For q G we put 

1 "1 
rkg = -{2Txk + ^{qj + -)aj) 

i=i 

In this section, Planck's constant takes the form, 

h = r^q- 

10.1. Local model. The model space is the normal bundle A^^ along 7, which may be 
identified with M" x S^. On the quantum level, the model Hilbert space is isomorphic to 
n = H'^iSl) ® L2(R"), where H'^{Sl) is the Hardy space of the circle of length L = L^. An 
orthonormal basis of L^(M") of joint eigenfunctions of Harmonic oscillators 

i, = i,{y,Dy) :=^-{Dl+y]) 

is provided by the Hermite functions Dq,q E N". Here, Dq is the Gaussian -Do(y) = e"^!?'!^, 
while Dq := C^Af \..A;9"Do(g e N'^), with Cg = (27r)""/2(g!)-i/2^ gl = q^\...qj and where 



Aj ■.= yj + iDy^ A* = yj - iDy^ 
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are the creation/annihiliation operators. An orthonormal basis of H'^{S\) ® L^(M") of joint 
eigenfunctions of ^ and the Ij is then furnished by 

iplq{s, y) ■= efc(s) ® Dg{y), e^is) := e'-^"''. 

This model needs to be adapted to 7 in the sense that x 5^ needs to be converted to 
the normal bundle A^-y. Without explaining this in detail, it may be illustrated by the fact 
that the adapted model transverse ground state Hermite function is 

Uo{s,u) = (detr(s))-i/2g*i(r(.K«>^ 

where Y = (yjk) is the matrix whose columns are a basis of vertical Jacobi fields along 7 
and where r(s) := ^Y~^. Higher transverse Hermite functions are obtained by applying 
standard creation/annihilation operators 



(adapted to 7): 

U, = Kf..Al-Uo. 

Roughly speaking, the metaplectic representation maps the model eigenfunctions on x 
to the normal bundle N^. Introducing the symplectic matrix 

ImY{s)* ImY{s)* 
ReY{sy ReY{sf 

with s G 5*1,, we introduce the unitary metaplectic operator (depending on the parameter s) 
fj,{a) := fj,{as)ds on Jg^ L'^{W)ds. Then fi{a) conjugates standard creation/annihilation 
operators to the adapted ones, and thus explains why the full quasi-mode is obtained by 
applying adapted creation/annihilation operators. We refer to |Z8 t IB.Bj for further details. 

10.2. WKB ansatz for a Gaussian beam. The most direct way to construct the Gaussian 
beam is to try to construct an approximate solution u = e'^^^U (s, yi, . . . , yn) of (A — A^)u = 0. 
Assuming that U is localized in a A"^/^ neighborhood of 0, one obtains a parabolic transport 
equation for f/, 

QTJ a2TT 

(90) CU := 2zX— + ^ ^ - A ^ K,h{s)y,yhU = 0, 

^ i=i j,h=i 

where Kjk are components of R{T, Y)T where R is the curvature tensor and where T = 7. 
The detailed ansatz {$^^(3, ^/Vk^y)} of |B.B] is defined by 

00 

(91) ^kq{s, ^y) = e*"^-«^ rljul (s, r'^) 

with U° = Uq and where are constructed inductively so that 

(92) A ^kq{s, y/r^y, r^g) ~ Afcg$fcg(s, ,/r^y), 
where 

(93) A,, , + £iM + PaW + .... 
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The numerators Pj{q) are polynomials of degrees j + 1- Thus one simultaneously constructs 
approximate eigenfunctions and eigenvalues around a stable elliptic geodesic. For each q and 
large k, the right side solves the eigenvalue problem to high approximation with a sequence 
of functions concentrated near 7. 

Instead of scaling the WKB quasi-mode, one can scale the Laplacian. More precisely 
(see [Z8]), one rescales the Laplacian around 7 in T^i^j) and conjugates by e*^'?*'^ using the 
operators (with h = rj^g), 

T,(/(s, u)) := h-^/'fis, h-'^u), M.ifis, u)) := e^^fi-s, y). 

We define the scaling of an operator A as Ah := T^M^AThMh- Then a simple computation 
gives 

n n 
ij=l i=l 

where the subscript [h] indicates to dilate the coefficients of the operator in the form, 
fh{s,u) := f{s,h2u). One then has, 

00 

m=0 

where £2 = L~^, £3/2 = and where 



r\ ^ It II 



j=i ij=i 
One then plugs into the ansatz ( 1911) to get 



/2 



=0 / i=0 j=0 

and solves recursively for Uq^"^ . A key point is that the metaplectic operator /i intertwines 
Dg to the operator 

n n 

C := fi{a)*Dsfi{a) = - -(^ + X! Kij{s)u^Uj) 



of am. 



2 

j=i ij=i 



10.3. Quantum Birkhoff normal form: intertwining to the model. A second ap- 
proach, which perhaps originates in [W] and which has been developed in many articles 
including |^ [Z8| \TZ\ ITZ2j . is the method of interwining operators. In this approach, the 
construction of quasi-modes follows from the dual construction of an interwining operator 
to a normal form of A around 7. The intertwining operator maps the model eigenfunctions 
to the detailed quasi-mode. 
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where ^kq solve ( 192|) . Rather than constructing the quasi- mode, one instead constructs 
so that it conjugates to a function of the model operators on the model space. We are 
assuming that 7 is elliptic, so the model operators are harmonic oscillators. 

This function is the quantum Birkhoff normal form. The normal form is implicit in (!92|) . 
i.e. the normal form is the same as the quasi-eigenvalues (p3|l interpreted on the operator 
level. The precise statement is: 

Theorem 10.1. There exists a microlocally unitary Fourier integral operator defined 
near the cone R'^'j generated by 7 in T*{S}^ x M") such that 

W-'^^AW^ = Pi(£, J^i, V) + Po(£, V, V) + . . . , 



where 



ri{L,, i^i, i-ynj = H ^ 1 TTT^ 



k= 



I^l , . . . , I^n) 



with p^^ for m=-l,0,l,..., homogenous of degree k-m in the variables (J^i, J^„) and of 
weight -1. The kth remainder term lies in the space ©^=0 02(fc+2-j)^^~-' 

Here, 0„\l/^ is the space of pseudo differential operators of order r whose complete symbols 
vanish to order n at {y,ri) = (0,0). Thus, the remainder terms are 'small' in that they 
combine in some mixture a low pseudo differential order or a high vanishing order along 7. 

If one plugs in the eigenvalues G 1^^ for the harmonic oscillators (J^, the operator expansion 
in Theorem (110. ip becomes the eigenvalue expansion in (l92l) . 



11. Birkhoff normal forms around closed geodesics 

The classical Birkhoff normal form of a Hamiltonian around a closed geodesic (or an in- 
variant torus) is a local approximation around 7 of the Hamiltonian by completely integrable 
Hamiltonians modulo errors which vanish to higher order around the orbit. One constructs 
local action variables Ij (Hamiltonians with 27r-perodic flows) and approximates the Hamil- 
tonian by a function Hk{I) up to an error l'^. 

On the quantum level one has two notions of order: 'order in W and vanishing order at 
7. As in Theorem 110.11 one constructs local 'action operators' and expresses H = ^/A as a 
polyhomogeneous function of the action operators modulo small remainders. 

In the general non-degenerate case, the normal form involves a greater variety of qua- 
dratic normal forms or 'action operators' than in Theorem 110. 1[ In addition to the elliptic 
action operator Ij there can also occur the real hyperbolic action operators Ij and complex 
hyperbolic (or loxodromic) action operators f'j^'^^ ^ jch,im^ 

In the elliptic case was a direct sum of rotations, and the quantum normal form of C 
had the form 

1 " 
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where the spectrum o"(P-y) — {e^*"^}. In the general non-degenerate case the normal form 
will similarly depend on the spectral decomposition of P^. Recall that, since is symplectic, 
its eigenvalues pj come in three types: (i) pairs p, p of conjugate eigenvalues of modulus 1; (ii) 
pairs p, p~^ of inverse real eigenvalues; and (iii) 4-tuplets p, p, p~^p~^ of complex eigenvalues. 
We will often write them in the forms: (i) e^*"^, (ii)e''=^^, (iii) e^'*^^*'^-'' respectively (with 
aj, \j, pjjUj e M), although a pair of inverse real eigenvalues {— e''^'^} could be negative. 
Here, and throughout, we make the assumption that P^ is non- degenerate in the sense that 

Ut,pr ^ 1, (Vpi e a(P^), (mi, ... , m2„) G N^"). 

Each type of eigenvalue then determines a different type of quadratic action, both on the 
classical and quantum levels (cf. [Ho, Theorem 3.1],[Ar]): 



Eigenvalue type 


Classical Normal form 


Quantum normal form 


(i) Elliptic type 




P := \a{Dl + y^) 


(ii) Real hyperbolic type 


= 2\yr] 


:= X{yDy + Dyy) 


(iii) Complex hyperbolic 
(or loxodromic type) 


ph,B. ^ 2/x(yi77i + 2/2^72) 
jch,Im ^ - y277l) 


jch,Re ^ ^(y^Dy, + Dy,yi + y^Dy, + Dy,y2), 



In the case where the Poincare map P^ has p pairs of complex conjugate eigenvalues 
of modulus 1, q pairs of inverse real eigenvalues and c quadruplets of complex hyperbolic 
eigenvalues, the linearized a/A will have the form: 

j=i j=i j=i 

The full quantum Birkhoff normal form is then given by the analogue of Theorem B of [Z.l]: 



Theorem B Assuming 7 non-degenerate, there exists a microlocally elliptic Fourier integral 
operator W from the conic neighborhood 0/ R+7 in T*{N^) to the corresponding cone in 
T*S^ in T*{S^ X W) such that 



V 

j=l j=l j=l 

■n { f<^ f e fh fh fch,Re fch,Im fch,Re fch,Im\ 

/^ll-'l) ■ ■ ■ ) -'p) -'1 ) • • • ) -'q ) -'1 J-*! )---)-'c '-'c J 

where the numerators Pj{If, . . . ,Ip, Ii, . . . , /c'*'^™) are polynomials of degree j+1 in the vari- 
ables (/^, . . . ^/c'i'''^"*) and where the kth remainder term lies in the space ©*^^^02(jt+2-j)^^~"' 
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Here, On'^^ is the space of pseudo differential operators of order r whose complete symbols 
vanish to order n at {y,r)) = (0,0). Thus, the remainder terms are 'small' in that they 
combine in some mixture a low pseudodifferential order or a high vanishing order along 7. 

11.1. Local quantum Birkhoff normal forms. Readers unfamiliar with microlocal Birkhoff 
normal forms might keep in mind a comparison to the local models used in §3.81 to study 
gradient estimates and ranges of eigenfunctions. The local models were spaces of constant 
curvature, a classical kind of comparison in geometric analysis. The quantum Birkhoff nor- 
mal forms are of a different nature: they are more global, since they are local around closed 
orbits rather than points. Moreover, they are dynamical and geodesic-based rather than 
curvature-based. The local models around regular orbits will be used in §12.41 to obtain 
localization results for quantum limits and in normal forms around singular closed or- 
bits will be used to deduce estimates on quantum integrable eigenfunctions. The local 
models are enumerated and described in §11.21 Note that they model operators are often 
Schrodinger operators with quadratic potentials rather than Laplacians for special metrics. 

Throughout this section, the parameter h = plays the role of Planck's constant. 

11.2. Model eigenfunctions around closed geodesies. Model quantum completely in- 
tegrable systems are direct sums of the quadratic Hamiltonians: 

• := h{Dyy + yDy) on L^(R+) (hyperbolic Hamiltonian) , 

• J*^ := h'^Dy + y'^ on L^(]R) (elliptic harmonic oscillator Hamiltonian), 

• i^'^ := h [{yiDy^ + y2Dy^) + v^(z/i£'2 - 1/2/^5,1)] on L'^{^^) (complex hyperbolic 
Hamiltonian) , 

• / := UDq, on 5*^ (regular Hamiltonian). 

The corresponding model eigenfunctions are: 

• Uh{y]\n) = \\ogn\-^l^ [c+{fi)Y{y) |y|-i/2+^^(^)/Hc_(fl)F(-y) |y|-i/2+iA(?i)A]. |c_(n)p+ 
|c+(n)|2 = 1; \{h) e R. 

• Ue{y;n,h) = h~'^^'^exp{-y'^/h) $„(^"^/^2/); n eN. 

• Uchlr,9;ti,t2,h) = |logn|-V2^(-i+iti(?i))Meit2(/i)0. ti{h)M{K) E R. 

• Ureg{0] m, K) = 6*™"^; m G Z. 

Here, Y{x) denotes the Heaviside function, $„(?/) the n-th Hermite polynomial and {r,6) 
polar variables in the {yi,y2) complex hyperbolic plane. 

The important part of a model eigenfunctions is its microlocalization to a neighborhood 
of X = ^ = 0, so we put: 

^(x; h) := Opni x{x)x{y)x{0 ) " u{y; K), 

where e > and x ^ C^{[—e, e]). In the hyperbolic, complex hyperbolic, elliptic and regular 
cases, we write iphiy'i K),'il}ch{y] h),4'e{y] and ipregiy,^) respectively. A straightforward 
computation [T2] shows that when ti{h),t2{h),nh,mh = 0{h) the model quasimodes are 
L^-normalized; that is 

(94) ||Op,(x(x)x(y)x(0)^(i/;^)IU^~i 

a.s h 0. Note that, although the model eigenfunctions above are not in general smooth 
functions, the microlocalizations are C°° and supported near the origin. 
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12. Quantum integrable Laplacians 

We now go into more detail on (globally) integrable Laplacians. As discussed in ^ 
quantum integrable Laplacians are the only examples where one has explicit formulae for 
eigenfunctions. They are also examples where one can explicitly determine the norms 
of the eigenfunctions. Moreover, there are geometric properties of the geodesic flow which 
account for the behavior of the norms. 

Quantum integrable Laplacians are fundamental for two reasons besides the relatively 
explicit computability of their eigenfunctions. The first reason, which is rigorous and now 
well-understood, is that quantum integrable Laplacians provide local models for all Lapla- 
cians around closed geodesies. The sense in which they are local models is made precise by 
quantum Birkhoff normal forms. The second reason why they are fundamental is heuristic: 
to date, all extremals for norms occur in integrable systems and are explained by the 
geometry of the geodesic flow. It is possible that eigenfunction sequences which achieve 
maximal growth rates must resemble sequences of eigenfunctions of integrable systems. 
At least, one hopes that they must do so locally, e.g. in the vicinity of a closed geodesic. 

In this section, we define quantum complete integrability and give some examples. In 
succeeding sections, we survey results on the norms and concentration properties of the joint 
eigenfunctions. 

12.1. Quantum integrability and ladders of eigenfunctions. Definition: We say 
that the operators Pj G \I'^(M); j = l,...,n, generate a semiclassical quantum completely 
integrable system if 

[P„P,] = 0; Vl<^,j<n, 

and the respective semiclassical principal symbols pi,...,pn generate a classical integrable 
system with dpi A dp2 A ■ ■ ■ A dpn 7^ on a dense open set fl C T*M\0. 

Semiclassical limits are taken along ladders in the joint spectrum. For fixed b = {bi, 62, &n) ^ 
M", we define a ladder of joint eigenvalues of the original homogeneous problem Pi = 
yA,P2,...,P„by: 

(95) Sfe := {(Aife, ...,A„fe) e 5'pec(Pi,...,P„); Vj = l,..,n, lim = bj}, 

fc— +00 Afc 

where |Afc| := ^/Xlk + - + >^lk- 

We define the joint eigenspace corresponding to Hb as follows: For b E B := V{T* M\0), 
where V is the classical moment map, define 

(96) Vb := {v?/.; W^P^^Wl^ = Iwith/x G SJ. 

12.2. Geometric examples. 

12.2.1. The round sphere. On the quantum level we define A = ^ and J2 = y^A + so 

that /2^ = ^ + The joint spectrum of Ji, 12 equals {(m, + |) : \m\ < k}. 

Let us compare classical and quantum values. The image of the classical moment map 
{pe{x,^), is the triangular region of points (x, y) satisfying < y,y > 0. It is easy to see 
that the inverse image of any point in the interior is invariant under the a;3-axis rotations and 
under the geodesic flow; this is a Lagrangian torus. The boundary of the image corresponds 
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to singular points of the moment map, namely the equatorial geodesic, traversed in either of 
its two orientations. 

The joint eigenfunction F^, with joint eigenvalue (m, k + ^), corresponds to the Lagrangian 
torus with pe{x,^) = m and |^| = k. If we rescale back to S*S'^ we obtain the Lagrangian 
torus Pq = m/k. The Lagrangian corresponding to zonal spherical harmonics is the meridian 
torus pe = 0, i.e. to longitudinal great circles which depart from the north pole, converge at 
the south pole and then return to the north pole. The highest weight spherical harmonics 
correspond to the boundary points of the triangular image of the moment map, hence to the 
equatorial great circle. Indeed, is the restriction of the harmonic polynomial (xi + ix2)^ 
(up to normalization). This polynomial is independent of 0:3 and is a holomorphic function 
of xi,X2 so it is certainly a harmonic homogeneous polynomial. Also, by its form it clearly 
is largest on the unit circle in the (xi, X2) plane and tends to zero as (xi, X2) — > (0, 0). Hence 
on the sphere it defines a spherical harmonic which is large on the equator and tends to zero 
at the poles. 

12.2.2. Simple surfaces of revolution |CV3] . The classical action variables can be quantized 
to produce quantum action variables Ij, which are first order pseudodifferential operators 
with the property that e^'^*^^ = Cjld for some constant Cj of modulus one. From the fact 
that e^'^*^-' = Cjld for a quantum torus action, it follows that the joint spectrum of the 
quantum moment map 

Sp{i) c n r + {/i} 

is the set of integral lattice points, translated by /i, in the closed convex conic subset F C M^. 
The vector fi = {^1,^2) can be identified with the Maslov indices of the homology basis of 
the invariant torii. In the case of -^A^, /i = (0, 1/2) |C V3] . 

The operator /i = 4^ is the standard infinitesimal generator of the 5*^ action, but the 

second action operator I2 is a rather unexpected operator that does not appear in the classical 
'separation of variables' approach to eigenfunction. We may define spectral projectors Hm 
projecting onto the eigenspace of I2 of eigenvalue A^. In the case of the round sphere it is the 
projection onto the eigenspace Hat of A of spherical harmonics of degree N , but in general 
it does not project onto eigenspaces of the Laplacian and is a rather novel feature of the 
pseudo-differential approach. 

One may express y^A^ as a function H{Ii,l2) where H{^i,^2) is a polyhomogeneous 
function on F. Its principal symbol is the function \^\g = H{I) expressing the metric norm 
function in terms of the action variables. These expressions are known as the Birkhoff normal 
forms of the classical, resp. quantum, Hamiltonian. In generic situations, the action variables 
are only locally defined and there exists only a local Birkhoff normal form. The interesting 
feature of simple surfaces of revolution is that the Birkhoff normal forms are global. 

One can then characterize the the packet of eigenvalues of -\/A which corresponds to the 
eigenfunctions with eigenvalue for I2: the eigenvalues form the set {H{m, N) : \m\ < N}. 

12.2.3. Hyperbolic cylinders with boundary |CVP] . A hyperbolic cylinder with boundary is 
a metric tube X around the unique closed geodesic 7 of a hyperbolic cylinder H/Z where 
Z acts by the cyclic group generated by a hyperbolic element. Its axis then projects to the 
closed geodesic 7. 
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As with any surface of revolution, the Laplacian commutes with the generator ^ of the 
action. Here, one puts either Dirichlet or Neumann boundary conditions on the boundary of 
the tube X. The novel feature is that the closed geodesic is hyperbolic and lies on a singular 
level of the moment map and the normal form of around 7 is given by /'^ := h{Dyy+yDy). 
The (not very close) analogue of the Gaussian beam or highest weight spherical harmonic 
is the sequence of eigenfunctions where the weight, i.e. -eigenvalue, is asymptotically 
the same as the \/A eigenvalue. Equivalently, if one separates variables to obtain Sturm- 
Liouville operators Lm for each weight, then the sequence is obtained by taking the smallest 
Lm-eigenvalue for each weight. 

12.3. Localization of integrable eigenfunctions. In this section, we consider the con- 
struction of highly localized or concentrated eigenfunctions (and quasi-modes). 

12.3.1. Toric integrable systems. Let A G ^°(M) denote any zeroth order pseudodifferential 
operator and dfi\ denote Lebesgue measure on the Lagrangian torus Tx. In the toric case 
we have the following localization theorem: 

Proposition 12.1. |Zlj For any ladder {kX + z/ : = 0, 1, 2, . . . } of joint eigenvalues, we 
have: 

{Aipkx,ipkx) = (TAdfix + 0{k~^). 

We thus have: 

Corollary 12.2. For any invariant torus Tx C S*M , there exists a ladder {({>kx-, k = 
0, 1, 2, . . . } of eigenfunctions localizing on Tx- 

We illustrate the localization result on S"^. The proof is the same in all toric quantum 
integrable cases, but is harder for general integrable systems for reasons explained in the 
next section. 

The image of T*S'^ — under the moment map fi{x,^) = (^^(x,^), |.^|) is a vertical trian- 
gular wedge. It is a cone, reflecting that fi{x,r^) = rfi{x,C,) is homogeneous. We can break 
the homogeneity by taking a base for the cone with |^| = 1, i.e. by considering points [x, 1). 
This corresponds to looking at pe : 3*3"^ ^ M. 

Thus, we consider pairs {rrij, kj) in the joint spectrum of De, A = ^^/A'^^T/2 — 1/2 whose 
projection to the base of the cone has a limit (c, 1). 

Theorem 12.3. Suppose thatruj/kj c. Then 

{Op{a)Y^,Yj^) ^ / aodx, 

where is the principal symbol of Op{a) . 

Thus, the eigenfunctions in this ray localize on the invariant torus ^(c). 
We define t/(ti, ^2) = e*^*^^''^*^"^'' and note that it is a unitary representation of the 2-torus 
T2 on L2(52). Further 



(Op(a)r^,r^) = {U{t,MrOp{a)U{t,M)YLY^). 
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Indeed, the eigenvalues cancel out. Average this formula over T2. We note that 

{A) := [ U{tiMyOp{a)U{h,t2)dtidt2 

commutes with both Dg and A. Indeed, the commutator with A gives ^ under the integral 
sign, and the integral of this derivative equals zero. 

But Dg, A have a simple joint spectrum: the dimension of the joint eigenspace equals one. 
Hence, any operator which commutes with them is a function of them. Thus, 

{A) = F{Dg,A). 

The function F must be homogeneous of degree zero. Also, the right side is a \E'DO. Its 
symbol is 

(ao) : / ao{^'"''{x,0)dtidt2. 

It follows first that 

(Op(a)y^,F^) = ((Op(a))F^,F^) 
= F{m,k). 

Secondly, as {mj,kj) oo with mj/kj c, we have F{mj,kj) F(c, 1). But also, the 
limit is the integral of Oq against an invariant measure. The principal symbol of F is (ao), 
which is a function on the image of the moment map. Its value at (c, 1) is by definition 
//i-i{c 1) '^o'^^;, concluding the proof. 

12.3.2. ]R"-integrable systems. It is a surprisingly large step from toric integrable systems 
such as the torus, sphere or simple surface of revolution to a general M" integrable system, 
even such as a surface of revolution with a hyperbolic orbit to a 'peanut' shaped surface in 
or an undulating surface of revolution with several 'waists'. From the dynamical point 
of view, there now exists a hyperbolic closed orbit (the waist of the peanut), and a cylinder 
of geodesies which asymptotically spiral toward the hyperbolic closed geodesic. Together 
they form a singular level set with several components. Eigenfunctions on this singular level 
Proposition (112.11) localize on level sets of the moment map rather than on individual tori or 
components. We prove this and further analyze the degree of concentration by putting the 
Laplacian into a quantum normal form around each component. We follow |TZ2j and use 
its semi-classical notation where h = \~^. 

Let 6 be a regular value of the moment map V, let 

V-\b) = AW(6)U---UA('"^')(6), 

where the A^^\b);l = 1, ...,m are ra- dimensional Lagrangian tori, and djjLi^a)^^^-^ denote the 

normalized Lebesgue measure on the torus A.^^\b). Define the semiclassical ladders 

(97) 

U,s{h) := {b,ih) := ibf\h),bfih), bfih)) e Spec(Pi, ...,P„); \b,ih)-b\ < Ch'-'}. 

Taking a sequence h 0, the joint eigenvalues in hb-silT') form a sequence tending to b which 
is the analogue of a homogeneous ladder. 
Define 



(98) ci{h;bj{h)) := {Opti{xi)Vh,{h),Vb,{h))] I = 1, mc/(6). 
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Here, xi is a cutoff function wliicli is equal to 1 in a neiglibourliood fi^'^ (6) of the torus A^'^ (6) 
and vanisfies on Uk=/=i^^''\b). 

Proposition 12.4. Let b G Breg, and let {^bj{h)} be a sequence of Li^-normalizeed joint 
eigenfunctions of Pi, Pn with joint eigenvalues in the ladder Lb^siK) ■ Then, for any a G 
^0,-00^ we have that as h 0: 

\l=l J JA^'Hb) 

Here, (i;UA(3)(fe) denotes Lebesgue measure on A^^\b). 

Tliis result says that a sequence of eigenfunctions of a quantum integrable system localizes 
on the level set of the moment map corresponding to the limit point in the image of the 
clsssical moment map, but it is not very precise about the weight attached to each component. 

12.4. Conjugation to normal form around torus orbits. We sketch the proof of 
Proposition 112.41 as an example of the use of normal forms. 

Let be a small neighbourhood of A*^'^ on which there exist action-angle variables 
(^^(0 /(0) Then there exists a 'microlocal' conjugation to quantum Birkhoff normal form. 
That is, for £ = l,...,k and j = l,...,n, there exist microlocally unitary /i- Fourier inte- 
gral operators, t/^^J : C~(M) ^ (:7°°(T"; together with symbols, ff\x;h) ~ 

Ek=ofjki^)^\ with /,o(0) = such that: 



n d 



(99) f/ir/f (A - b^'\ Pn - M"); h)Ui^l =,(0 

The simple operators on the right side are the normal forms. 

Conjugation to normal form shows that the space of (microlocal) solutions of equations 

(100) Pk^b,{h) =nm{b) bf\K)^b,(n) 

is one-dimensional. Indeed, such solutions are the same as solutions of 

/«(Pi - 6«, P„ - h)^, =^(0(,) /f (6? ^ - b^'\ &f ^ - h)^,. 
After conjugation to Birkhoff normal form (l99ll the equation becomes 

h d 

I 89 j 

and the solutions are just multiples of exp[i{n + 71^/4)6], where 7 is the Maslov index and 
n G Z. Thus, the joint eigenfunctions (pbj{h) are given microlocally by 

(101) =ni'Hb) \lci{h-Mm f/g(e^("^+-^/^)^). 

The right sides of fllOll) are the usual quasimodes or semiclassical Lagrangian distributions 
[CV2] 

We further use the normal form to prove the localization statement. Let Xii^iC) ^ 
C^{T*M);l = 1, ....,mci{b) be a cutoff function which is identically equal to one on the 
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neighborhood and vanishes on ^l^^\b) for k ^ I. For h sufficiently small, we then have 



In each term, we now conjugate to normal form. Recalling the definition fl98|) of the weights 
Cj{h, bj (h)), we have 



In the last step we used the Egorov theorem to recognize Ull*Opfi{a) o Opfi{xi)Ull as a 
pseudo-differential operator on the standard torus and then used the calculation in (IMIl to 
complete the proof. For the details of the remainder estimate we refer to |TZ2j . 

13. Concentration and non-concentration for general {M,g) 

In this section we go over results about concentration on, and norms along, submani- 
folds of general Riemannian manifolds. Then we again compare the results to the quantum 
integrable case. 

13.1. LP norms for restrictions to submanifolds. A measure of concentration near a 
submanifold is given by the norms of the restrictions of eigenfunctions to the submanifold. 
Early results were obtained by Reznikov [R] for curves on hyperbolic surfaces. The most 
general results are due to Burq-Gerard-Tzvetkov [BGTj and pertain to submanifolds of 
general Riemannian manifolds. On curves, the estimates in |BGT] are as follows: 

Theorem 13.1. jBGT] Let Wifx^Wi^ = l- Then 
• (i) If 'J is a unit-length geodesic, then 



1=1 



(102) 



ci{h- b,{h)){Ug*Opn{a) o Op;,(x/)f/ge*("^+"^/')^ e^K+-7/4)e^ 




7 



where 



(103) 




• (a) If 'J is a curve with strictly-positive geodesic curvature, 




7 



LOCAL AND GLOBAL ANALYSIS OF EIGENFUNCTIONS ON RIEMANNIAN MANIFOLDS 



63 



13.2. Non-concentration in tubes around hyperbolic closed geodesies. We now con- 
sider results which give upper bounds on the concentration of eigenfunctions in tubes around 
hyperbolic closed orbits and other invariant sets. The article |CVP] gave such an bound and 
a general principle was then proved in [BZJ (see also |Chrj ) . It is stated in terms of lower 
bounds of eigenfunction mass outside of a tube around a hyperbolic closed geodesic. 

Theorem 13.2. \BZ\ IChr] Let {M,g) be a compact Riemannian manifold, and let 'y be a 
hyperbolic closed geodesic. Let U be any tubular neighborhood of 7 in M . Then for any 
eigenfunction (px, there exists a constant C depending only on U such that 

I IV'A|W,>-^||v.A||i.. 

Jm\u iog ^ 

More generally, let A G \E'°(M) be a pseudo-differential orbit whose symbol equals one in 
a neighborhood of 7 in S*M and equals zero outside another neighborhood. Then for any 
eigenfunction (p\ 

C 
vlogA 

These theorems are proved in |BZ] using a rather abstract observability estimate: 

Theorem 13.3. Let P{h) be a family of self-adjoint operators on a Hilbert space Ti, with 
a fixed domain V. Let Hi be a second Hilbert space and suppose there exists a family of 
bounded operators A{h) : T) 7ii satisfying 



\u\\n<^\mh) + T)u\\n + gm\Mh)u\ 



Hi, 



for all T e {-b, -a) and with 1 < G{h) < Ch^^" for some No. Let x e C^{-b, -a). Then 

G{h) 
T(h) 



there exist constants cq, Co o.nd Hq > so that for any function T{h) satifying Ettt < cq, 



MPWn < Co'^£^'^ \\A{h)e-^'^('y\{Pih))u\\l,^dt. 

This formulation has applications to non-concentration of eigenfunctions in the Buni- 
movich stadium as well. 

13.3. Non-concentration around closed geodesies on compact hyperbolic surfaces. 

We sketch the proof in [CVP] in the case of hyperbolic closed geodesies on a compact 
hyperbolic surface X = H\r. 

Suppose 7 is a simple closed geodesic whose normal bundle is trivial. Then it possesses a 
tubular (collar) neighborhood Xq which is isometric to a collar around a closed geodesic 7 
in a hyperbolic surface of revolution H\r^ ~ 5*^ x My. We noted above that such a collar 
is quantum integrable and that the Birkhoff normal form of \/A around 7 is the hyperbolic 
operator yDy + Dyy. Using the explicit formulae for the eigenfunctions in the model domain, 
one can calculate the amount of mass of the model eigenfunction near 7, and finds that the 
mass increases logarithmically in A (cf. jC VP] . Proposition 16). The mass for the singular 
energy level E is explicitly given by 

2 



(104) /o*? r^-^^^-^^^^^Mf)^ ~|r(i + zf)|^e^|ln/^|. 



-1 , .Ex|2 ^421, 
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The mass of the sequence of eigenmodes of the hyperbohc surface of revolution with boundary 
grows at the same logarithmic rate, since it is microlocally unitarily equivalent to the model 
eigenfunctions. 

Now consider the eigenfunctions of the compact hyperbolic surface without boundary. 
Unlike the hyperbolic cylinder with boundary, its geodesic flow is ergodic and far from 
completely integrable. Let ip\^ denote a sequence of eigenfunctions, and let 

k 

be the Fourier series expansion of ip\^ for each y. Let | |afe;Aj I |[-i,i] denote the norm of the 
k component in [—1, 1] and let Hv^jHxo denote the norm in the collar. Then we have 



Suppose now that the sequence {v^Aj } concentrates in Xq in the sense that 

'aJIxo 



||v^A,||xo>^o>0, ||<^A,|||,->0, VirccX\7. 



Then, 

(105) J]||aA,A, ||[_i,i] >Co>0 

k 

uniformly in j. Then the only possible quantum limit measures of the sequence must have 
the form C/x^ + z/ where z/ = in the tube around 7 and C > 0. For a time-reversal symmetric 
symbol a with essential support in the collar, we must have 

{pp{a)ipxp^x,) ^ C / a^ds. 

Now let P = A~^-^ + 1, a pseudo-differential operator of order whose principal symbol 
equals 1 — The symbol vanishes along 7. Then 

(106) / y.A,P¥^A,rfV = 5^(l-^)||a,;A,||[-i,i]^0. 

Let Zj = {k: |§ - 1| < i}. Combining (fTO and ffTU^ . we have 

(107) 5^||afc;A,|P>C>0. 

Thus, the Fourier coefficients of a sequence of eigenfunctions concentrating near 7 concentrate 
around the joint eigenvalues of the modes on the hyperbolic surface of revolution which 
concentrate on the central geodesic. The mass around 7 of the associated modes can only 
grow logarithmically must charge the complement of any collar around 7 with mass at least 
of order j^^. It follows from (11071) that any sequence of eigenfunctions on the compact 

surface without boundary must also charge the complement of the collar by at least 
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14. U' NORMS AND CONCENTRATION IN THE QUANTUM INTEGRABLE CASE 

1/ norms of eigenfunctions in the quantum integrable case have a much more geometric 
theory than in general. Exphcit examples show that the they are often extremals for 
norm and concentration inequalities. 

One extremal problem raised in |YH IY2] is to determine the Riemannian manifolds which 
possess orthonormal bases of eigenfunctions with uniformly bounded L°° norms. An obvious 
example is a flat torus. The question arises whether any others exist. The following result 
shows that flat tori are the unique minimizers in the class of quantum integrable Laplacians. 

Theorem 14.1. |TZ] Suppose that A is a quantum completely integrable Laplacian on a 
compact Riemannian manifold {M,g). Then 

(a) IfL'^{X,g) =0(1) then {M,g) is flat. 

(b) Ifr{X,g) = 0(1), then {M,g) is flat. 

It is possible that fiat tori are the only compact Riemannian manifolds with a uniformly 
bounded orthonormal basis of eignfunctions. But at this time, it is not even known if the 
standard sphere possesses such an orthonormal basis. The idea of the proof of Theorem 114. II 
is that the joint eigenfunctions concentrate on level sets of the moment map, and therefore 
develop singularities at points where the projection of the level set to the base manifold 
is singular. The only case where no singularities occur is when (M, g) is a torus without 
conjugate points, and in this case Burago-Ivanov proved (the Hopf conjecture) that (M, g) 
must be fiat. 

There also exists a quantitative improvement of Theorem 114.11 which gives blow-up rates 
for norms for quantum integrable eigenfunctions concentrating on singular level sets, i.e. 
level sets which are not regular in the sense of ([7]). These eigenfunctions are the extremals 
for LP blow-up and mass concentration. In the following there is an additional technical 
assumption (Eliasson non- degeneracy) which we omit for simplicity. 

Theorem 14.2. jTZ2] Suppose that {M,g) is a compact Riemannian manifold whose Lapla- 
cian A is quantum completely integrable as in H^) . Then, unless {M,g) is a flat torus, this 
action must have a singular orbit of dimension < n. If the minimal dimension of the singular 
orbits is i, then for every e > 0, there exists a sequence of eigenfunctions satisfying: 



• A point {x,C,) is called a singular point of the moment map V if dpi A ■ ■ • A 

dPnix,0 = 0. 

• A level set V~^{c) of the moment map is called a singular level if it contains a singular 
point (x,0 e V-^ic). 

• An orbit M" ■ (x,^) of $f ([8]) is singular if it is non-Lagrangean, i.e. has dimension 

< n; 

The idea in the proof is to consider singular orbits and to conjugate to a quantum Birk;hoff 
normal form around the orbit. Hence, one calculates the mass in the normal form space. 




Here, 
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Conjugation to normal form does not preserve norms, and what is really calculated are 
not norms but rather norms in shrinking tubes. Interestingly, this method produces 
optimal results. The proof does not determine the minimal dimension i. By taking products 
of lower dimensional manifolds, it is easy to construct examples with any value of £ = 
l,...,n-l. 

14.1. Mass concentration on small length scales. We sketch the proof of Theorem 114.21 
as an illustration of mass estimates in shrinking tubes. We follow the semi-classical notation 
h= of [TZ2 J. 

Let A := M" ■ f be a compact, k < m-dimensional singular orbit of the Hamiltonian Re- 
action generated by In this section, we study mass concentration of modes in 
shrinking tubes of radius ~ for < 5 < 1/2 around vr(A) in M, where vr : T*M — ^ M 
denotes the canonical projection map. 

We denote by T^{tt{A)) the set of points of distance < e from vr(A). For < 5 < 1/2, we 
introduce a cutoff xi{x; h) G C^(M) with < Xi < 1, satisfying 

. (i) supp XiCT^.(7r(A)) 
. (ii) x? = lonT3/4?,.(7r(A)). 
• (iii) \d'^x{{x]h)\ < c„/i-^i"i. 
Under the assumption that A is an embedded submanifold of M, the functions 

(108) xiix;h) = Ciih-''d\x,n{A))) 

are smooth on T^{7r{A)) and satisfy the conditions. Here, d{., .) is the Riemannian distance 
function. Also, Ci e C^(M) with < Ci < l,Ci(a;) = 1 for |x| < 3/4 and supp Ci C (-1, 1). 

Theorem 14.3. Let G Vdh) satisfy the hounds in LemmaW^ Then for any < 5 < 1/2, 
(OMx?)^M.<^M)>|log^r"^- 

We briefly sketch the proof. Let xH^^^] K) G C^(T*M] [0, 1]) be a second cutoff supported 
in a radius tube, Q{h), around A with Q{h) C suppxi and such that = 1 on suppxl,- 
Then, clearly 

(109) xi{x,o>xi{x,o. 

for any (x, ^) G T*M. By Garding's inequality, fll09p implies 

(110) ( Opn{xi)¥>^., ) > ( Opn{xi)'P^., <^/. )• 

We now conjugate the right side to the model by the h- Fourier integral operator F of 
Lemma (1991) . Since F is a microlocally elliptic ^-Fourier integral operator associated to a 
canonical transformation k, it follows by Egorov's theorem 

(111) (Opnixi)^,,^,) = \c{h)\\Opn{xi ° ~ C.h^''' 

where c{K)u^{ii, 9; h) is the microlocal normal form for the eigenfunction ip^. Since ip^ G Vdh) 
satisfies the bounds in Lemma fl93|) it follows that |c(/i)p ^ |log^|~™ and from (II 111) we 
are left with estimating the matrix elements {Opnixl, ° 1^)^(1, u^) from below. As in the 
localization proof, the matrix elements are now in terms of elementary model eigenfunctions 
and the calculation has become easy. The normal form eigenfunctions separate into a product 
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of factors and one only has to calculate one (or two) dimensional integrals. As an example, 
in the hyperbolic case the integral has the form 

|2 



(112) > _^ (log/i)-i I «e--x-V2+*V^dx ' + C»(|log^|-i) 

> \T{l/2 + i\/h)mi - 2S) + 0{\\ogh\-^) >C(e) >0 
uniformly for h G (0, ho{e)]. 

14.1.1. Completion of the proof of Theorem \1 4 -2^ The small scale mass estimates immedi- 
ately imply lower bounds on L°° norms and norms due to the shrinking volumes of the 
tubes. For instance, 

Im \'^f^i^)?xi{x; h) dvol{x) < sup^gr^^^(^(A)) ImXH^; K) dvol{x) 

(113) 

< WvAl^- jMxi{x]h) dvol{x) 

and it follows from Lemma [14.31 that 



(114) \WA\loo-\^j^x{{x;h) dvol{x)j >C{e)\\ogh\- 
uniformly for h G (0, ho{e)]. Since 

(115) / xi{x]h) dvoiix) = Oih^^''~^'y), 



J M 

f[TT4D implies 

> C(e)n-5("-^)+^|log^|-'". 
Recalling that hr^ G {A^; Aj G Spec — V^}, this gives: 

\Wx,\\loo > Cie)X~\ 

14.1.2. Concentration of quantum integrable eigenf unctions on submanifolds . Similar meth- 
ods were used in |To2] to obtain sharp bounds on norms for restrictions to submanifolds 
in the quantum integrable case, making more precise the results of |BGT] in this special 
case. For simplicity, let us consider curves on surfaces. First is the generic upper bound: 

Theorem 14.4. |To2j Let ^x^j = 1,2,3,... he the L?' -normalized joint Laplace eigenfunc- 
tions of the commuting operators Pi = — A and P2 on a Riemannian surface {M'^,g). Then 
for a generic curve 7 such that L*P2\s*M is Morse, we have 



/ l^Xjl^ds = (logAj) 



When the curve is a geodesic, the bounds depend on the type of level set the geodesic lies 
on: 

Theorem 14.5. |To2j Let Pj{h)]j = 1,2 be a nan- degenerate quantum integrable system 
system on a surface, {M,g). Then, 
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• (i) When 7 is the projection of a geodesic segment contained in V ^{Breg), 

I |y.,^.(s)prfs = 0|,|(l), 

• (a) When 7 is the projection of a singular joint orbit in V~^{Bsing), 

Moreover, there exists a constant c^, > depending only on the curve 7, and a sub- 
sequence of joint eigenfunctions, ^\j^ = ] k = 1,2, ... such that 

(fx^^{s)\'^ds > c^A]^^ when'y is stable, 

f\j^,{s)\'^ds > c^A^^^I log Aj,, when is unstable. 

Thus the exact bound depends on the nature of the geodesic. In the general quantum in- 
tegrable case, most geodesies lie on regular Lagrangian tori in V~^{Breg) and these geodesies 
do not support large L^-bounds. But as in Theorems 114. II and I14.2[ there always exists a sub- 
sequence of joint eigenfunctions of Pi and P2 with mass concentrated along (singular) orbits 
contained in V~^{Bsing), and the associated eigenfunctions saturate the upper bounds. For 
instance in the case of a simple surface of revolution, the equator is the projection of a singu- 
lar orbit of the action generated by geodesic flow and rotation. The corresponding joint 
eigenfunctions (the analogs of highest weight spherical harmonics) satisfy \(p\^\'^ds ~ X^"^ 
along the equator, 7. The equatorial geodesic is singular and the norms along it had 
singular blowup. In the case of the meridian great circles, the closed geodesic lies in the base 
space projection of a maximal Lagrangian torus. The zonal harmonics have ^-microsupport 
on this torus and have L^-restriction bound ~ log A along any meridian great circle. 

15. Delocalization in quantum ergodic systems, I 

In this section, we discuss general results on eigenfunctions when the geodesic flow of 
{M,g) is assumed to be ergodic (see and §6.41 for definition and notation). The study 
of eigenfunctions and eigenvalues of Laplacians on manifolds with ergodic (or more highly 
mixing) geodesic flows is generally known as 'quantum chaos'. The basic question is, what 
impact do dynamical properties of the geodesic flow g'^ have on eigenvalues and eigenfunctions 
of its quantization f/^? This question has been studied over the last three decades by a 
large collection of mathematicians and physicists, using both theoretical and computational 
methods. In this section, we largely follow our recent survey |Z3j . Another exposition with 
emphasis on arithmetic hyperbolic quotients is [Sar2j . For recent computational results, we 
refer to [Bar] . 

One of the basic and most studied problem is Problem 16.61 for quantizations of classically 
ergodic systems. The main result is that there exists a subsequence {^pj^.} of eigenfunctions 
whose indices jk have counting density one for which pj^ (A) := {A{pj^, ipj^^) — > uj{A) (where as 
above uj{A) = .^^^^ Is*m ^Adfi is the normalized Liouville average of a a)- Such a sequence 
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of eigenfunctions is called a sequence of 'ergodic eigenf unctions'. The key quantities to study 
are the quantum variances 

(116) ^^^^)^=iv7Ay 5Z i(M,^.)-^(^)r- 

Theorem 15.1. jSETl ISEl [Z2l lUVl [Sul IZZ^ ICn IZ^ Let{M,g) be a compact Riemanman 
manifold (possibly with boundary), and let {Xj,(pj} be the spectral data of its Laplacian A. 
Then the geodesic flow G* is ergodic on {S*M, d/i) if and only if for every A G 'i/°{M), we 
have: 



(1) limA_ooV:4(A) = 0. 

(2) (Ve)(35)limsup;,^^^^,/fe:A^,A,<A \{A^j,^k)\'^ <e 

Since all the terms in (1) are positive, no cancellation is possible, hence (1) is equivalent 
to the existence of a subset 5 C N of density one such that Qs := {d^k : k E S} has only u 
as a weak* limit point. 

As explained in |Z9j . this ergodicity of eigenfunctions may be viewed as a convexity the- 
orem: By the Banach-Alaoglu theorem, the set of invariant probability measures for the 
geodesic flow is a compact convex set Aii. An invariant measure is ergodic if it is an ex- 
treme point of the compact convex set. The same is true on the quantum level: The set of 
£^]R of invariant states for at is a convex set. A classical invariant measure is an invariant 
state, and if it is ergodic classically it is also ergodic quantum mechanically, i.e. it is an 
extreme point of S^. Hence Liouville measure u is an extreme point of this convex set. But 
the local Weyl law says that u is the limit of the convex combination j:^je)'^\j<e Pj- 
extreme point cannot be written as a convex combination of other states unless all the states 
in the combination are equal to it. In our case, cu is only a limit of an infinite sequence of 
convex combinations, and the result is that almost all terms in the sequence tend to u, and 
that is equivalent to (1). 

Sketch of Proof of (1) 

Let 

(117) ^A)t:=^ j\;AUtdt. 
Then, 

(118) J]|(Ay,,^.,v.,^.)-.;(A)|2= J]|((AK-.;(A)y,,^,,y,,^.)p. 

Aj<A Aj<A 

Apply the Schwartz inequality for states. 



|(5^A,,¥^A,))r<Trn[o,A]5*5, 

j:Aj<A 
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where T^[o,x] is the spectral projection for vA corresponding to the interval [0, A], to the 

operator B = Hp a] [(^)t — ^^(^)]n[o a]- We then have 

(119) 

Ea,<a I((^)t - u;{A)^^^,^,X < Tr {U[o,x][{A)t - a;(A)n[o,A])in[o,A] (A)^ - ^(A)]n[o,A]) 

< Tr (n[o,A] [{A)t - oj{A)y[{A)T - ^(A)]n[o,A]) 

= ujm)r-uj{A)r{{A)r-uj{A))). 
Here, we used the Jensen inequality, 

^ -Try;(n[o,A]KA)T-^(^)]n[o,A]) < -^Trn[o,A]</^([(A)T-cc'(^)])n[o,A], 



valid for any convex function Lp in the case ^p{x) = . By the local Weyl law, we get 
li^^El((^^A,,¥^A,)-^(A))^< / 

A-»oo N[X) Js*M 

As T ^ oo the right side approaches by the dominated convergence theorem and by 
BirkhofF's ergodic theorem. Since the left hand side is independent of T, this implies the 
stated theorem. □ 

15.1. Quantum ergodicity in terms of operator time and space averages. To explain 
the term 'quantum ergodicity', we reformulate the result in terms of space and time averages. 
We assume for simplicity the generic condition that all eigenvalues are of multiplicity one 
(cf. [U]). The space and time averages are defined as follows: 
Definition Let A ^ be an observable and define its time average to be: 

kdt 



{A) ■.= w- hm / U:AUt> 



and its space average to be scalar operator 

u{A) ■ I 

The limit is taken in the weak operator topology, i.e. in the sense of matrix elements. To 
verify it we observe that 

iy:AU,dt^,,^,) = '^^^^^^{A^,,^,), 

hence the matrix element tends to zero as T — > oo unless Aj = Xj. The limit only occurs 
in the weak sense since the rate is clearly not uniform, as the spacings Aj — Xj could be 
arbitrarily small. 

Quantum ergodicity can thus be reformulated as the condition, 

(120) {A) = uj{A)I + K, where lim uJx{K*K) ^ 0, 

where oJ\{A) = TrE{X)A. Thus, the time average equals the space average moduli a term K 
whose Hilbert-Schmidt norm in the range of IIa is o{N{X)). Note that (A) commutes with 
-\/A, hence is diagonal in the basis {(pj} of joint eigenfunctions of (A) and of Ut- Also, K is 
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the diagonal matrix with entries {A(pk, fk) — uj{A). The condition is therefore equivalent to 
Va{\) 



15.2. Quantum unique ergodicity and converse quantum ergodicity. A Laplacian 
is said to be QUE (quantum uniquely ergodic) if Liouville measure is the only weak * limit 
point of the semi-classical Wigner measures. The terminology was introduced in |RSj . 

The condition may be reformulated in terms of K: Namely, QUE is equivalent to the 
compactness of K. Indeed, this would imply that {K(pk,ipk) —>■ 0, hence {Aipk,(pk) ^{A) 
along the entire sequence. A key difficulty of settling the question whether K is compact 
is that the time averaged operator (A) no longer belongs to the class of pseudo-differential 
operators, due to the very weak nature of the weak operator limit. 

It is widely conjectured there exist non QUE (M, g) with ergodic geodesic flow. We refer 
to |Bar[ IBZ2[ [Z9j for some recent discussions and references to the literature. The simplest 
example is that of the Bunimovich stadium, which possesses quasi-modes concentrating on 
the invariant Lagrangian cylinder (with boundary) formed by bouncing ball orbits in the 
middle rectangle. An analogue among Riemannian manifolds without boundary is that of a 
non-positively curved surface with a flat cylindrical part and ergodic geodesic flow [Donlj . 
which also carries product quasi-modes of the same kind. The existence of such quasi-mode 
suggests that there are nearby modes which are not quantum ergodic (see |Z9] for the precise 
suggestion). Faure-Nonnenmacher-de Bievre |FNBj have shown that QUE does not hold for 
the hyperbolic system defined by a quantum cat map on the torus, and since the methods 
available for studying eigenf unctions of quantum maps and of Laplacians are very similar, 
this negative result shows that there cannot exist a universal structural proof of QUE. A QUE 
result has been proved by E. Lindenstrauss, namely the QUE property for the orthonormal 
basis of Laplace-Hecke eigenf unctions eigenf unctions on arithmetic hyperbolic surfaces (see 
[LINDj l 

To the author, an interesting and almost completely open problem is the converse of 
quantum ergodicity: does quantum ergodicity imply classical ergodicity? More precisely, 
are there natural sufficient conditions for this? For instance, if (A) = uj{A) + K where K 
is compact, is the geodesic flow ergodic? Very little is known on this converse problem at 
present. One may imagine a system which is non-ergodic but in which 'tunnelling' causes 
eigenfunctions to become uniformly distributed in phase space. A simple model is the one- 
dimensional Schrodinger operator with the even double well potential V{x) = (1 — x^)^. 
Level sets -]-V{x) = E with E near zero have two components and hence the Hamiltonian 
flow is not ergodic on this level; yet the eigenfunctions are either even or odd and hence they 
are quantum ergodic. No such example is known in the case of Laplacians A of compact 
Riemannian manifolds. However, the example shows that no general 'abstract' proof of 
QE =^ CE is possible. 

In |Z5t IMUZj it is shown that if there exists an open set in S*M flUed by periodic orbits, 
then the Laplacian cannot be quantum ergodic. But it has not even been proved that KAM 
systems, which have Cantor-like positive measure invariant sets, are not quantum ergodic 



although there exists a positive density of quasi-modes concentrating on invariant tori [Pop 
The problem is to prove that eigenfunctions are linear combinations of not too many quasi- 
modes. 
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15.3. Quantum weak mixing. Quantum weak mixing concerns tlie off-diagonal matrix 
elements. 

Theorem 15.2. (see [Z3] for references) The geodesic flow $* of {M,g) is weak mixing if 
and only if the conditions (l)-(2) of Theorem \15. 1\ hold and additionally, for any A G \E'°(M), 

(Ve)(35)limsup^ U^,.^k)?<e (Vr G M) 

|Aj-Afe-T|<« 

The restriction j ^ k is of course redundant unless r = 0, in which case the statement 
coincides with quantum ergodicity. This result follows from the general asymptotic formula, 
valid for any compact Riemannian manifold {M,g), that 

(121) 

In the case of weak-mixing geodesic flows, the right hand side — > as T — * oo. 

15.4. Spectral measures and matrix elements. Theorem 115.21 is based on expressing 
the spectral measures of the geodesic flow in terms of matrix elements. The main limit 
formula is: 



smT{Xi~Xj~T) 
T{X,-X,-t) 



(122) / df,^,:= \im—- ^ \{Aip,,ip,)W 

Xi,Xj<X, |Ai-Aj-T|<e 

where dfi„^ is the spectral measure for the geodesic flow corresponding to the principal 
symbol of A, a a G C°°{S*M,dfi). Recall that the spectral measure of Vt corresponding to 
/ G is the measure dfif defined by 

{Vtf,f)LHS'M)= [ e^'^df^jir). 
The limit formula fll22p is equivalent to the dual formula (under the Fourier transform) 

(123) i^WX) ^ e''^'^-'^^\{A^^,^J)f = {VtaA,aA)LHs*M). 

^ ' i,j:Xj<X 

The proof of (fT23D is to consider, for A G the operator A^A G ^° with At = UlAUf By 
the local Weyl law, 

1 

a";^ N{X) 

The right side of fll22p defines a measure dmA on M and fll23p says 



lim -^;j^TiE{\)AlA= {Vt(JA,OA)L^(s*M) ■ 



e'^^dmAir) = {VtaA,(TA)LHS'M) = / e^*^rf/x^^(r). 

Jr 

Since weak mixing systems are ergodic, it is not necessary to average in both indices along 
an ergodic subsequence: 
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(124) lim {A'lAip^^ip,) = Ve**(^-^^)|(A¥;„¥;,)|2 = {VtaA,aA)L^is*M). 

j 

Dually, one has 

(125) lim J2 = / ^'df^'^A- 

Xj^OO ^ ' Jt-£ 

i : \Xi — Xj—T\<e 

For QUE systems, these limit formulae are valid for the full sequence of eigenfunctions. 

15.5. Rate of quantum ergodicity and mixing. A quantitative refinement of quantum 
ergodicity is to ask at what rate the sums in Theorem ll5.1( l) tend to zero, i.e. to establish a 
rate of quantum ergodicity. In the off-diagonal case one may view \{A(pi,(pj)\'^ as analogous 
to \ {A(pj, (fj) —uj{A)\'^. However, the sums in (11221) are double sums while those of (11161) are 
single. One may also average over the shorter intervals [A, A + 1]. 

The only rigorous result valid on general Riemannian manifolds with hyperbolic geodesic 
flow is the logarithmic decay: 

Theorem 15.3. |Z5] (see also |Schu2] for p = 2) For any {M,g) with hyperbolic geodesic 
flow, 

The proof uses the central limit theorem for geodesic flows of M. Ratner. The logarithm 
reflects the exponential blow up in time of remainder estimates for traces involving the wave 
group associated to hyperbolic flows. It would be surprising if the logarithmic decay is sharp 
for Laplacians. It was shown by R. Schubert |Schu] that the estimate is sharp in the case of 
two-dimensional hyperbolic quantum cat maps. Hence the estimate cannot be improved by 
semi-classical arguments that hold in both settings. 

A stunning asymptotic formula for Va(A) was proved by Luo-Sarnak for holomorphic forms 
of arithemetic hyperbolic quotients. Before stating the result, we review some conjectures 
in the physics literature. 



15.6. Quantum chaos conjectures. First, consider off-diagonal matrix elements. One 
conjecture is that it is not necessary to sum in j in (I125p : each individual term has the 
asymptotics consistent with (I125p . This is implicitly conjectured by Feingold-Peres in jFPj 
(11) in the form 

(126) K^,..,,)P^^L_|_i. 

where Ca(t) = f'^^e~'^'^^{VtaA,crA)dt. In our notation, Xj = h~^Ej and p{E)dE ~ dN{X). 
There are ~ CA""^ eigenvalues Aj in the interval [Xj — r — e, Xj — r + e], so (11261) says that 
individual terms have the asymptotics of (11251) . 
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On the basis of the analogy between \ {Aipi, v?j)P and \ {A(pj, ipj) —uj{A)\'^, it is conjectured 
in [FP] that 

The idea is that !f± = :^(v^t =t (fj) have the same matrix element asymptotics as eigenfunc- 
tions when Aj — Xj is sufficiently small. But then 2{Aip^, ip_) = {Acpi, ipi) — {Aipj, ipj) when 
A* = A. Since we are taking a difference, we may replace each matrix element by {Aipi, ipi) 
by {Aipi, ifi) —uj{A) (and also for ipj). The conjecture then assumes that {Aipi, ipi) —u^A) has 
the same order of magnitude as {Aipi.ipi) — {A{pj,ipj). The order of magnitude is predicted 
by some natural random wave models, as discussed below in ^ IT8l 



15.7. Rigorous results. At this time, the strongest variance result is an asymptotic formula 
for the diagonal variance proved by Luo-Sarnak for special Hecke eigenfunctions on the 
quotient H^/ SL{2, Z) of the upper half plane by the modular group |L.S.2j . What they prove 
is an asymptotic variance formula for holomorphic Hecke eigenforms. One expects that their 
proof, suitably modified, extends to smooth Maass-Hecke eigenfunctions. Therefore, as in 
|Z3] . we describe the statement for smooth eigenfunctions as a Theorem/Conjecture - i.e. it 
is a Theorem for holomorphic forms, but still a conjecture for non-holomorphic forms. Note 
that H^/S'L(2, Z) is a non-compact finite area surface whose Laplacian A has both a discrete 
and a continuous spectrum. The discrete Hecke eigenfunctions are joint eigenfunctions of A 
and the Hecke operators Tp. 

Theorem/Conjecture 1. |L.S.2j Let {ipk} denote the orthonormal basis of Hecke eigen- 
functions for H^/S'L(2, Z). Then there exists a quadratic form B{f) on C^(H^/S'L(2, Z)) 
such that 

When the multiplier / = is itself an eigenfunction, Luo-Sarnak have shown that 

5(¥.A,¥^A)=C^,(0)L(i,v9A) 

where L{^, ipx) is a certain L-function. Thus, the conjectured classical variance is multiplied 
by an arithmetic factor depending on the multiplier which has no dynamical significance. At 
this time, it is unknown whether variance asymptotics exist in the non-arithmetic case. From 
numerical experiments, it is believed that arithmetic settings behave somewhat differently 
from non-arithmetic ones (see |Sar2j ). and this could be another example of the non-generic 
behavior of arithmetic quantum chaos. 

15.8. Quantum limits on a hyperbolic surface and Patterson-Sullivan distri- 
butions. In this section, we mention a curious link between quantum limits and classical 
dynamics on a hyperbolic quotient that was observed in [AZj . It is related to the invariant 
triple products studied in [Sari IBRt IMSj . 

We write G = PSU{1, 1) := SU{1, 1)/ ±1 = PSL{2,R),K = PS0{2) and identify the 
quotient G/ K with the hyperbolic disc D. We let F C G denote a co-compact discrete group 
and let Xr = F\D denote the associated hyperbolic surface. In this context it is standard 
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to denote eigenvalues by A| = Sj{l — Sj) = \ + i"j {sj = ^ + irj) and eigenf unctions by 

{<Pirj}j=0,l,2,...- 

We recall that on a hyperbolic quotients there exists a quantization a Op{a) of symbols 
which is adapted to the non- Euclidean Fourier transform of Helgason (see f|T71) and [AZJ for 
references and details). We then define the Wigner distributions by 

(127) {a,W,r,)= I a{g)W^rM9)■■={0p{a)^^r,,^^r,)L^^,), aGC°°(5*Xr) 

On the other hand, one can define a second sequence of phase space distributions, the 
Patterson-Sulhvan distributions {PSir-} associated to the eigenf unctions {v^jr^}, by the ex- 
pression 

Tir\dh)Tir\dh') 

(128) PS,rM9) = PS,rMb',db,dt) := ' ® \dt\. 

In this definition, Tiy. is the boundary values of ipi^- in the sense of Helgason (see fl2T|) ). The 
parameters (6', h) {h ^ b') vary in B x B, where B = (9D is the boundary of the hyperbolic 
disc, and t varies in M; (6', b) parametrize the space of oriented geodesies, t is the time 
parameter along geodesies, and the three parameters {b',b,t) are used to parametrize the 
unit tangent bundle 5'D. The Patterson-Sullivan distributions PSir^ are by construction 
invariant under the geodesic flow {g^) on ST), i.e. 

(129) {g%PS„,=PS„,, 

and by using ( !22l) they can be shown to be F-invariant. Hence to each eigenfunction 
one obtains a geodesic-fiow invariant distributions on 5'Xr. We also introduce normalized 
Patterson-Sullivan distributions 

(130) FS,r, := , PS.r,, 

which satisfy the same normalization condition {l,PSir.) = 1 as Wir- on the quotient 5'Xr. 
In |AZ] the following is proved: 

Theorem 15.4. For any a e C°°(F\G), 

a{g)W,r^{dg) = ! a{g)Fs,r,{dg) + 0{rj'). 

It follows that the Wigner distributions are equivalent to the Patterson-Sullivan distribu- 
tions in the study of quantum ergodicity. Yet, the Patterson-Sullivan distributions have a 
purely classical dynamical definition: Define the classical dynamical zeta functions. 



(131) 



1 10 



(zz) C{s-a):=Y.,-^{j\,o), (Ree.>l) 



where the sum runs over all closed orbits, and 70 is the primitive closed orbit traced out by 
7. The sum converges absolutely for Ree s > 1. 
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Theorem 15.5. Let a be a real analytic function on the unit tangent bundle. Then C{s; a) 
and £2(5; a) admit meromorphic extensions to C. The poles in the critical strip < Ree s < 
1, appear at s = 1/2 + ir, where as above 1/4 + is an eigenvalue of A. For each zeta 
function, the residue is 

{a, PSirj)sXr, 

where {PSir } are the normalized Patters on- Sullivan distributions associated to an orthonor- 
mal eigenbasis {ipir^}. 

Thus, the quantum hmit problem is the same (for compact hyperbohc surfaces) as the 
problem of finding the limiting behavior of the residues of classical dynamical zeta functions 
as the pole moves up the critical line. 

16. Delocalization of eigenfunctions: II: Entropy of quantum limits on 

MANIFOLDS WITH AnOSOV GEODESIC FLOW 

We now describe the recent results of Anantharaman |A] , Anantharaman- Nonnenmacher 
|AN] (see also [ANK] ) on entropy of quantum limits for {M,g) with Anosov geodesic flow. 
These articles give lower bounds on entropies of the quantum limit measures that arise from 
sequences of eigenfunctions. We closely follow the presentation in |ANl lANKj . and refer in 
particular to the partly expository article [ANK] and to the recent Bourbaki seminar report 
|CV4j for an exposition of the results. 

As discussed above, it is known that quantum limits are invariant probability measures for 
(?*, but many such measures exist for any geodesic flow. To pin down the possible quantum 
limits, one needs to add constraints on the possible limit measures. The entropy bounds of 
Anantharaman et al provide almost the only additional constraints known at this time. The 
lower bound on entropy rules out such possible limit measures as periodic orbit measures 
(which have entropy zero) or finite sums of such measures. However, the entropy results 
leave open the possibility that a sequence of eigenfunctions could tend to a limit of the form 
a/i-y + (1 — a)dfi if a is small enough (here, n is Liouville measure). 

We recall (see ^Dthat a geodesic flow g^ is Anosov on on S*M if the tangent bundle TS*M 
splits into g^ invariant sub-bundles E^{p) © E^{p) ® R-^i/(p) where is the unstable 
subspace and E^ the stable subspace. The unstable Jacobian J"(p) at p is defined by 
J''(p) = det ^El(gi^)). 

The goal is to give a lower bound for entropy of quantum limits. Entropy is complicated to 
define, and we only provide a brief sketch here. Classically, entropies are defined for an invari- 
ant probability measure p for the geodesic flow and measures the average complexity of p- 
typical orbits. In the Kolmogorov- Sinai entropy, one starts with a partition V = {Ei, . . . , E^) 
of S*M and deflnes the Shannon entropy of the partition by hp{p) = J2j=i f^i^j) ^og f^i^j)- 
Under iterates of the time one map g of the geodesic flow, one reflnes the partition to 

P^" = {E^, n g-'E^, n ■ ■ ■ n g''^+'E^^_,}. 

One deflnes hn(V, p) to be the Shannon entropy of this partition and then deflnes h^sifJ', 'P) = 
hm„^oo ^KifJ', V). T hen hKsjlA = supp hKsi^-, 'P)- 
The main result of |AN^ lANKj is the following 
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Theorem 16.1. Let jj, he a semiclassical measure associated to the eigenfunctions of the 
Laplacian on M . Then its metric entropy satisfies 

(132) hKs{^)> ^ (^-1) 



- Ar 



2 - ^max 5 



I S*M 

where d = dimM and Amax = lini|f|_^oo ^logsup^g^- {dg^l is the maximal expansion rate of the 
geodesic flow on E. 

In particular, if M has constant sectional curvature —1, this means that 

(133) hKs{li) > 

The proof is based on a quantum analogue of the metric entropy, and in particular on a 
development of the 'entropic uncertainty principle' of Maassen- Uffink. There are several 
notions of quantum or non-commutative entropy, but for applications to eigenfunctions it is 
important to find one with good semi-classical properties. 

Let {H, (., .)) be a complex Hilbert space, and let Wi/jW = ^J {ip, ip) denote the associated 
norm. The quantum notion of partition is a family vr = {T^k)k=i,...,N of operators on Ti such 
that Yl!k=i ^fc^fe = Id. If ll?/'!! = 1, the entropy of ip with respect to the partition vr is define 

by 

k=l 

We note that the quantum analogue of an invariant probability measure /i is an invariant state 
p, and the direct analogue of the entropy of the partition would be ^ pijikT^l) logp(7rfc7r^). If 
the state is p{A) = {Aip.ip) then pinki^l) = W'n'li'W'^- 

The dynamics is generated by a unitary operator U on Ti. We now state a simple version 
of the entropy uncertainty inequality of Maasen- Uffink. A more elaborate version in [XI lANt 
lANKj gives a lower bound for a certain 'pressure'. 

Theorem 16.2. For any e > 0, for any normalized (zTi, 

K{Ui)) + K{^1}) > -21ogc(W) , 

where 

c{U) = sup |(efc, Uej)\ 

is the supremum of all matrix elements in the orthonormal basis {ej}. In particular, h.^{ilj) > 
— \ogc(U) if ip is an eigenfunction ofU. 

In the application to eigenfunctions, one fixes a partition {M^} of M and a corresponding 
partition T*Mk of T*M. One then defines a smooth quantum partition of unity Pk by 
smoothing out the characteristic functions of M^. The partition is refined by 

(134) P« = P„„_, {n - l)Pa,Mn - 2) ■ ■ ■ P„„, 
where a = {ao, . . . , ttn-i) and where P{k) = W^PU'^. One then specifies: 

(1) 14 = e*"^^^ is the wave operator at the 'Ehrenfest time' Tr = ^^sA, Qr from a 
semi-classical (where h = where the Hamiltonian is H = h'^A and where the time 
evolution is e'i^), U = e^"^('^)'^^ with nE(h) = 
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(2) iph is an eigenfunction of lA ; 

(3) Km =T.\a\=nJPM?^Og\\P:,IJn\?. 

With these specifications, 

(135) c{U)= max | |Pa't/"^PaOp(x^"^^)| I 

\a\=\a'\=nE 

where x^"^'' is a very sharp energy cutoff supported in a tubular neighborhood := H~^{1 — 
e, 1 + e) of £ = S*M of width for a given 5 > 0. 

We now give a very sketchy outhne of how the entropy uncertainty inequahty (Theorem 
116. 2p is used to prove the lower bound on the entropy of the limit measure (Theorem 116. ip . 
The argument is technical and difficult and the outline only gives the fiavor of the estimates; 
the interested reader should consult [ANK] for a complete exposition. 

The next step is to link c{U) to the classical dynamics. The authors introduce a discrete 
'coarse-grained' unstable Jacobian 

jr(«o,«i) := sup{ J^(p) : p e T*Q^,nS' : g'p G T*fi,J, 

for ao, ai = 1, . . . , K. Here, Qj are small open neighborhoods of the partition sets Mj. For 
a sequence a = {ao, . . . , of symbols of length n, one defines 

J^{a) := Ji{ao, ai) ■ ■ ■ J{'(a„_2, ftn-i)- 

Theorem 16.3. Given a partition V^^^ and 5, 5' > small enough, there exists hp(o) such 
that, for any h < ^(o) 55/, for any positive integer n < riEih), and any pair of sequences a, 
a' of length n, 

(136) \\P:, f/"P,Op(x("))|| < Cr('^-i+^^) Vj«(a)J„"(a') . 
Here, d = dimM and the constants c, C only depend on {M,g). 

To prove this, one shows that any state of the form Op(x^*'')^ can be decomposed as a 
superposition of essentially h 2— normalized Lagrangian states, supported on Lagrangian 
manifolds transverse to the stable leaves of the fiow. The action of the operator on such 
Lagrangian states is intuitively as follows: each application of U stretches the Lagrangian 
in the unstable direction (the rate of elongation being described by the unstable Jacobian) 
whereas each multiplication by Pa^ projects onto a small piece of the Lagrangian. This 
iteration of stretching and cutting accounts for the exponential decay. 

Combined with the entropy uncertainty inequality, one obtains 

Proposition 16.4. Let d = dimM and let {ipf^jfi^Q he a sequence of eigenf unctions. Then 
there exist 5, 5' so that, at time n = nE{K), 

(137) h^{i,n) > 2(d - 1 + c5) log h + 0{l) > _2 ^d - 1 + c5)X^^, ^ ^ ^^^^ 

(1 - d') 

Now suppose that the Wigner measures W^,^ of a subsequence {'ipfi)h~>o of eigenfunctions 
converges to the semiclassical measure /i on S. Consider the limit h ^ of hTj{iph){so that 
nE^h) — > 00). For any sequence a of length ue, each ||P* converges to /i({a}), where 
{a} is the function P^^ (P^^ o g^) . . . (P^^ o 51"^) on T*M. Then for any Hq < Ue, hno{iph) 
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semiclassically converges to the classical entropy 

|a|=no 

where 

Here, is a smoothing of the characteristic function of the indicated set. 
Using Proposition 116.41 one obtains the lower bound 



138) 



Un no (1 - d') Uo 

ao,ai 



R 



Here, 5 and 5' could be taken arbitrarily small, and at this stage they can be set equal to 
zero. 

The Kolmogorov-Sinai entropy of /i is by definition the limit of the left side of (11381) when 
Uo —>■ oo. Then let no —>■ oo, and let the diameter £:/2 of the partition tend to 0. Then the 
first term in the right hand side of (11381) converges to the integral — log J'^{p)dfi{p) as 
£ — i> 0, proving (I132p . 

□ 



17. Real analytic manifolds and their complexifications 

In this section, we consider eigenfunctions on a real analytic compact Riemannian manifold 
{M,g). The advantage of real analyticity is that one can complexify the manifold and 
analytically continue the eigenfunctions to the complexification of M. This allows one to use 
methods of holomorphic and pluri-subharmonic function theory to obtain sharper results on 
volumes and distribution of nodal hypersurfaces than are possible for C°^{M,g). The gain 
in simplicity is two fold, reflecting the relative simplicity of real polynomials over smooth 
functions, and of complex zeros of polynomials over real zeros. This point of view has been 
taken in |DF[ ILinl IZ6j among other articles. 

A real analytic manifold M always possesses a unique complexification generalizing the 
complexification of as C™. The complexification is an open complex manifold in which 
M embeds l : M ^ as a totally real submanifold (Bruhat- Whitney). As examples, we 
have: 

• M = WjlT" is Mc = 'C^jlT'. 

• The unit sphere defined by + ■ ■ ■ + x'^_^^i = 1 in W^'^^ is complexified as the 
complex quadric Sq = {{zi, . . . , Zn) G C""*"^ : zl + ■ ■ ■ + z^_^i = 1}. 

• The hyperboloid model of hyperbolic space is the hypersurface in M"+^ defined by 

JT = {xl + ---xl- x^+i = -1, x„, > 0}. 

Then, 

= {(^1; • • • ; ^n+l) ^ C""^^ : zf + ■ ■ ■ Z^ — Z^j^^ = —1}. 

• Any real algebraic subvariety of M"* has a similar complexification. 

• Any Lie group G (or symmetric space) admits a complexification Gc- 
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The Riemannian metric determines a special kind of distance function on Mc \GS1\ IGS21 
ILSll ILS2| IGLSj . The metric g determines a plurisubharmonic function ^/p = ^/p^ on Mq as 
the unique solution of the Monge- Ampere equation 

Here, SM^,dVg is the delta-function on the real M with respect to the volume form dVg, i.e. 

/ — i> j^.jfdVg. In fact, yp(C) = "^V^^^ICTO where r'^{x,y) is the squared distance function 
in a neighborhood of the diagonal in M x M. 

One defines the Grauert tubes = {C G : y/p{C) ^ t}- There exists a maximal tq 
for which ^/p is well defined, known as the Grauert tube radius. For t < tq, is a strictly 
pseudo-convex domain in Mq. 

The complexified exponential map (x, ^) exp^iC, defines a diffeomorphism from B*M to 
Mj. and pulls back ^ to \^\g. The one-complex dimensional null foliation of dd^/p, known 
as the 'Monge- Ampere' or Riemann foliation, are the complex curves t + ir T7(i), where 
7 is a geodesic, where r > and where T'j{t) denotes multiplication of the tangent vector to 
7 by r. We refer to [LSI] for further discussion. 

17.1. Analytic Continuation of eigenfunctions. Let A{t) denote the operator of an- 
alytic continuation of a function on M to the Grauert tube Mr. It is simple to see that 
A{t) = UciiT)e'^^ where U{iT,x,y) = e~'^^{x,y) is the Poisson operator of {M,g), i.e. 
the wave operator at positive imaginary time, and Uc{iT,(,y) is its analytic continuation in 
X to Mt-. In terms of the eigenfunction expansion, one has 

oo 

(139) U{tT,C,y) = J2e-^'^iff{0vj{y), iC,y)eM,xM. 

j=0 

To understand the analytic continuability of the wave kernel, we first consider Euclidean M" 
and its wave kernel U{t, x, y) = Q^Ai\ ^^{i^^-y) which analytically continues to t + ir, = 
X + ip & C+ X C" as the integral 

Uc{t + iT,x + ip, y)= [ e'(*+*")l«le'<«'"+'P-^>c/e 

The integral clearly converges absolutely for \p\ < r. 

Exact formulae of this kind exist for S"™ and H™. For a general real analytic Riemannian 
manifold, there exists an oscillatry integral expression for the wave kernel of the form, 

(140) U{t,x,y)= [ e**l«l««e'<«'^"P«"'(")U(t,a;,i/,0(ie 

Jt*m 

where A(t,x,y,^) is a polyhomogeneous amplitude of order 0. For background, we refer to 
|Bet ID.Gt [T] . The holomorphic extension of fll40p to the Grauert tube |C| < t in x at time 
t = ir then has the form 

(141) Uc{iTX,y)= I e--l«l«^e^«''=-P«'^(^))A(t,C,l/,0^e {C = x + ip). 

Jt* 

y 

Since 



(142) 
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the analytic continuability of the Poisson operator to imphes that every eigenfunction 
analytically continues to the same Grauert tube. It follows that the analytic continuation 
operator to is given by 

(143) Ac{T) = Uc{ir)oe^. 

Thus, a function / G C°°(M) has a holomorphic extension to the closed tube ^/p{C) ^ 
if and only if / G Dom{e^^), where is the backwards 'heat operator' generated by 

\/A (rather than A). That is, / = Yl'^=o^n'^Xri admits an analytic continuation to the open 
Grauert tube M^. if and only if / is in the domain of e"^^, i.e. if |a„pe^'^'^" < oo. Indeed, 
the analytic continuation is Uic{iT)e^^ f . The subtlety is in the nature of the restriction to 
the boundary of the maximal Grauert tube. 

This result generalizes one of the classical Paley- Wiener theorems to real analytic Rie- 
mannian manifolds jBoutlGS2] . In the simplest case of M = S*^, / ~ 'Ylm&z '^n^^^^ G C"^(S'^) is 
the restriction of a holomorphic function F ~ Xlnez '^n^" the annulus S"^ = {| log | < r} 
and with F G L'^{dS].) if and only if |/('^)P e^'"''^ < oo. The case of is more compli- 
cated since it is non-compact. We are mainly concerned with compact manifolds and so the 
complications are not very relevant here. But we recall that one of the classical Paley- Wiener 
theorems states that a real analytic function / on M" is the restriction of a holomorphic func- 
tion on the closed tube |ImC| < r which satisfies \F{x + i^)\'^dx < C ioi ^ < r ii and 

only if />|i"^«l G L2(M"). 

Let us consider examples of holomorphic continuations of eigenfunctions: 

• On the flat torus R'^/Z™', the real eigenfunctions are cos(/c, x), sin(/c, x) with k G 
27rZ™'. The complexified torus is C^/Z™ and the complexified eigenfunctions are 
cos(fc, C), sin(fc, Q with = x + i^. 

• On the unit sphere S**", eigenfunctions are restrictions of homogeneous harmonic 
functions on M"*"*"^. The latter extend holomorphically to holomorphic harmonic 
polynomials on and restrict to holomorphic function on S*™. 

• On H"*, one may use the hyperbolic plane waves e'^*'^"'"^)^^'''^, where (2, h) is the (signed) 
hyperbolic distance of the horocycle passing through z and h to 0. They may be 
holomorphically extended to the maximal tube of radius 7r/4. 

• On compact hyperbolic quotients H^/F, eigenfunctions can be then represented by 
Helgason's generalized Poisson integral formula [Hj . 

Jb 

Here, z E D (the unit disc), B = dD, and dTx G T>\B) is the boundary value of 99^, 
taken in a weak sense along circles centered at the origin 0. To analytically continue 
^px it suffices to analytically continue (2, h). Writing the latter as (C, 6), we have: 

(144) ^^(C) = / e(*^+^)«'''>rfTA(6). 

Jb 

The results of Sarnak [Sarj on the exponential decay of integrals /^^^ ipxip^^dVg and subse- 
quent results of Miller-Schmidt |MS] and Bernstein-Reznikov |BR] are closely related to the 
Paley- Wiener theory and this analytic continuation formula for eigenfunctions on hyperbolic 
quotients. 
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17.2. Maximal plurisubharmonic functions and growth of ip"^. There are natural 
analogues in the setting of Gruaert tubes for the basic notions of pluripotential theory on 
domains in C"^. One may view the Grauert tube function y^p as the analogue of the pluri- 
complex Green's function or Siciak maximal PSH (pluri-subharmonic) function. 

In the case of domains Q C C"^, we recall that the maximal PSH function (or pluri-complex 
Green's function) relative to a subset E G Q is defined by 

Ve{0 = snp{u{z) : u G PSH{Q),u\e < 0,u\an < 1}. 

This maximal function controls the Bernstein constant B{f,E,il) = j j.j for any holo- 
morphic /. 

An alternative construction of the maximal PSH function due to Siciak is defined by taking 
the supremum only with respect to polynomials p. We denote by the space of all complex 
analytic polynomials of degree and put = {p E : \\p\\k < 1, WvWn < e}. Then 
define 

log$^(C) =sup{— log|pjv(C)l -.peV^}, log^s = limsuplog$^. 

Here, = sup^g^ |/(z)|. Siciak proved that log'l'^; = Ve (see [K], Theorem 5.1.7). 

Intuitively, there are enough polynomials that one can obtain the sup by restricting to 
polynomials. 

On a real analytic Riemannian manifold, the natural analogue of is the space 

= {p= ^ O'j'fXj, ai, • • • , ttNix) e M} 

j:\j<X 

spanned by eigenfunctions with frequencies < A. Rather than using the sup norm, it is 
convenient to work with based norms than sup norms, and so we define 

7V(A) 

= {p= "j^^^' wpWl^m) = Yl i^jI^ = 

i:Aj<A j=l 

We define the A-Siciak extremal function by 

<^Uz) = sup{|V^(z)|V^: i; e Hx; Mm ^ 1}, 
and the extremal function by 

<I>m{z) = snp<^lj{z). 

X 

We may also define a natural analogue of the pluri-complex Green's function by putting 
E = M and = and defining 

VgiCr) = snp{u{z) : u G PSH{M^),u\m < 0,u\9m^ < r}. 

Although it does not seem to have been proved at this time, it is easy to guess that Vg = 
since the latter solves the homogeneous Monge-Ampere equation (dd^)"^ = on Mt-\M. 
Moreover, it is not hard to prove that 

(145) $M = Vg, 

generalizing the so-called Siciak-Zaharjuta theorem in the special case where the boundary 
conditions are placed on all of M. To see this, we consider the analytic continuation of 
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the spectral projections kernels Il[o^x]{x,y) = J2j-Xje[o \]Vj{^)fjiy)- Its complexification 
evaluated on the anti-diagonal equals, 

(146) n[o,A](C,C)= E l^i(Or 

j:XjeO,X] 

By using a Bernstein- Walsh inequality 



it is not hard to show that 

(147) $m(^)= lim |logn[o,A(C,C)- 

To evaluate the logarithm, one can show that the kernel is essentially times the tem- 
perate projection defined by the Poisson operator, 

(148) ^'[o,A](C,C)= E e''-^^^^^'^\^fiC)\'. 

j:A,e[0,A] 

The equality (11451) follows from the fact that limA^oo j log-P[o,A](C) C) = 0. 

17.3. Analytic continuation and nodal hypersurfaces. We now survey some methods 
and results on nodal hypersurfaces of eigenf unctions on real analytic compact Riemannian 
manifolds (M, g) of dimension m. The principal result on volumes is due to Donnelly- 
Fefferman [DFj : 

Theorem 17.1. |DFj (see also |Linj ) Let {M,g) be a compact real analytic Riemannian 
manifold, with or without boundary. Then there exist Ci,C2 depending only on {M,g) such 
that 

(149) ciA < n^-\Z^J < C2A, (Aipx = AVa;ci,C2 > 0). 

A very readable exposition of the proof is contained in [HJ. The upper and lower bounds 
require rather different arguments. The upper bound is simpler and may be sketched as 
follows: By a local Crofton's formula, the real volume of the nodal hypersurface equals 
the mean number of intersections it has with a random line in a coordinate chart, i.e. by 
the number of zeros of ifx along each line. This number is obviously bounded above by the 
number of complex zeros of on the complexification of the line. Hence, the upper bound is 
reduced to bounding the number of complex zeros of a family of complex analytic functions 
of one variable in a disc. For each disc, the number of zeros is in principle estimated by 
Jensen's formula, which bounds the number of zeros in a disc of a holomorphic function by 
the logarithm of the modulus of the holomorphic function. One can then use the frequency 
function estimates, Carleman estimates or Bernstein- Walsh type inequalities to obtain the 
doubling estimates in Theorem 13.21 and they bound the growth of the log modulus of <y9^ 
by ex. Jensen's formula does not directly apply, since it concerns the growth of a single 
holomorphic function in an expanding family of domains, while we are interested in a family 
of holomorphic functions y?^ in a single domain M^. Donnelly- Fefferman find a suitable 
replacement using Blashcke product factorizations. Lin gives an alternative bound on zeros 
in terms of growth using the frequency function for functions of one complex variable. Thus, 
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the upper bound reflects the growth rate of Lp\ combined with the upper bound on the 
number of complex zeros by the growth rate. 

The lower bound is of a somewhat different nature. A general analytic function of expo- 
nential growth (e.g. e^) need not have any zeros, i.e. zero might be an exceptional value. 
Jensen's formula equates the growth of the log modulus to the sum of the number of zeros 
and the 'proximity to zero', i.e. in our setting to values where |v5a(C)I ^ e~^^ It is necessary 
to rule out the possibility that zero is an exceptional value of complexified eigenfunctions. 
In effect, this is possible because real eigenfunction must have a zero in any ball of radius y 
by Theorem 14. 1[ But then one needs a lower bound on the hypersurface volume of the nodal 
set in such a small ball. More precisely, in [DF] a lower bound for the hypersurface volume is 
proved for a certain A-independent proportion of a covering by small balls of radius y". The 

global result then follows by summing the volume in the these small balls. Besides |DF] . we 
refer to |HLj for a detailed discussion. 

17.4. Nodal hypersurfaces in the case of ergodic geodesic flow. We now consider 
global results when hypotheses are made on the dynamics of the geodesic flow. Use of 
the global wave operator brings into play the relation between the geodesic flow and the 
complexified eigenfunctions, and this allows one to prove gobal results on nodal hypersurfaces 
that reflect the dynamics of the geodesic flow. In some cases, one can determine not just the 
volume, but the limit distribution of complex nodal hypersurfaces. 
The complex nodal hypersurface of an eigenfunction is defined by 

(150) Z^c = {C G 5;M : y,^(C) = 0}. 

There exists a natural current of integration over the nodal hypersurface in any ball bundle 
B^M with e < eo , given by 

(151) ([Z^^l^) = ±[ 5aiog|^^rA^= / V, ^eV(^-'^--'\B:M). 

In the second equality we used the Poincare-Lelong formula. The notation p(™~^'™~^)(i?*M) 
stands for smooth test (m — 1, m — l)-forms with support in B*M. 

The nodal hypersurface Z^c also carries a natural volume form \Z^c\ as a complex hyper- 

surface in a Kahler manifold. By Wirtinger's formula, it equals the restriction of ^Jl^_iy to 
Z^c. Hence, one can regard Z^c as defining the measure 

^ Jzc "^-1! 

We prefer to state results in terms of the current [Z^c] since it carries more information. 
We will say that a sequence {fji^} of L^-normalized eigenfunctions is quantum ergodic if 

(153) ^,,) - -^^^^ 1^^^^ aAdfi, VA G v&°(M). 

Here, \E''*(M) denotes the space of pseudodifferential operators of order s, and dfi denotes 
Liouville measure on the unit cosphere bundle S*M of (M, g). More generally, we denote by 
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d/ir the (surface) Liouville measure on dB*M, defined by 
(154) di^r = -TTTT on dB;M. 

We also denote by a the canonical action 1-form of T*M. 

Theorem 17.2. Let (M, g) be real analytic, and let {v^jj,} denote a quantum ergodic sequence 
of eigenf unctions of its Laplacian A. Let {B*^M, J) be the maximal Grauert tube around M 
with complex structure Jg adapted to g. Let e < eq. Then: 

-^[Z^c ] ^ -dd^ weakly m V'^^'^\B*M), 
in the sense that, for any continuous test form ip e 'D^"^~^'"^~^^{B*M), we have 



^Jk JZ r JB^M 



Equivalently, for any ip G C{B*M), 



\k Jz c i^- 1)' ^ JBtM [m- 1)! 



Corollary 17.3. Let {M,g) be a real analytic with ergodic geodesic flow. Let {fj^} denote 
a full density ergodic sequence. Then for all e < Cq, 

-^[Z r ] ^ -dd^, weakly in V'^^^^\B*M). 

^jk 

The proof consists of three ingredients: 

(1) By the Poincare-Lelong formula, [Z^c] — idd log \ if'^\. This reduces the theorem to 
determining the limit of ^log |v^^|. 

(2) i log 1 I is a sequence of PSH functions which are uniformly bounded above by 
^/p. By a standard compactness theorem, the sequence is pre-compact in L^: every 
sequence from the family has an L^ convergent subsequence. 

(3) when properly L^ normalized on each dMr is a quantum ergodic sequence on 
OMt-. This property implies that the L^ norm of |v^^P on dQ is asymtotically ^/p. 

(4) Ergodicity and the calculation of the L^ norm imply that the only possible L^ limit 
of J log \ (p'x\- This concludes the proof. 

We note that the first two steps are valid on any real analytic (M, ^f) . The difference is 
that the norms of ip'^ may depend on the subsequence and can often not equal ^/p. That 
is, j\p'x\ behaves like the maximal PSH function in the ergodic case, but not in general. 
For instance, on a flat torus, the complex zero sets of ladders of eigenfunctions concentrate 
on a real hypersurface in M£. This may be seen from the complexified real eigenfunctions 
sm{k,x + i^), which vanish if and only if {k,x) e 27rZ and (k,^) — 0. Here, k e N"* is a 
lattice point. The exact limit distribution depends on which ray or ladder of lattice points 
one takes in the limit. The result reflects the quantum integrability of the flat torus, and a 
similar (but more complicated) description of the zeros exists in all quantum integrable cases. 
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The fact that j log \ (p'^\ is pre-compact on a Grauert tube of any real analytic Riemannian 
manifold confirms the upper bound on complex nodal hypersurface volumes. 

17.5. Analytic domains with boundary. Many of the basic estimates on nodal hyper- 
surfaces, such as Theorem I149[ apply to real analytic manifolds with boundary as well as 
boundaryless manifolds. We are concentrating on manifolds without boundary, but mention 
one result on nodal lines where the boundary effects are central. Namely, in the case of real 
analytic plane domains, it is visible from computer graphics (see |FGS] for references) that 
only a small proportion of the nodal components touch the boundary. Most are small nodal 
loops in the interior. 

We thus consider Neumann (resp. Dirichlet) eigenfunctions (fx. on (piecewise) real analytic 
plane domains C M^, i.e. solutions of ([T]) with boundary conditions 

(155) d^(px, = (resp. (px, = 0) on dQ, 

Theorem 17.4. |TZ3] Let Q be a piecewise analytic domain and let nan(Aj) be the number 
of components of the nodal set of the jth Neumann or Dirichlet eigenfunction which intersect 
dfl. Then there exists Ca such that naai^j) < CfiXj. 

In the Dirichlet case, we remove the boundary before counting components. For q generic 
piecewise analytic plane domains, zero is a regular value of all eigenfunctions (pxj, i.e. Vv^a^ 7^ 
on Z<^^ [U], and the nodal set is a disjoint union of connected components which are 
homeomorphic either to circles contained in the interior Q° of Q or to intervals intersecting 
the boundary in two points. We call the former 'closed nodal loops' and the latter 'open nodal 
lines'. Thus, the theorem states that the number of open nodal lines is 0{Xj). As mentioned 
above, the Courant nodal domain theorem implies that the number of nodal componants 
is of order O(A^). When the upper bound is achieved, the number of open nodal lines in 
dimension 2 is of one lower order in Xj than the number of closed nodal loops. This effect is 
known from numerical experiments of eigenfunctions and random waves [FGSj . 

The proof of Theorem 117.41 is to complexify the Cauchy data (i.e. (pxlan in the Neumann 
case or d^ipxldn in the Dirichlet case) as holomorphic functions on the complexification of 
the boundary, and to count the number of its complex zeros. This number is bounded by the 
growth rate of the log modulus, which can be directly estimated in terms of the eigenvalue. 
Thus, the proof is of a similar form to the upper bound half of |DF]. On the other hand, no 
comparable lower bound exists: in the case of the disc (see §2.61) . the boundary values are 
of the form sinm^ or cos m6, where m is the angular momentum rather than the frequency 
A. So the number of boundary zeros remains bounded for some sequences of eigenfunctions 
as the eigenvalue tends to infinity. 

As previously mentioned, I. Polterovich [PoJ used Theorem ll7.4l to resolve an old conjecture 
of Pleijel, to the effect that Pleijel's bound on the number of nodal domains (see Theorem 
14. 6p is valid for Neumann as well as Dirichlet boundary conditions. Polterovich observes 
that Pleijel's argument applies to any nodal domain that does not touch the boundary. By 
Theorem 117. 4[ the nodal domains which do touch the boundary are of lower order than the 
bound, concluding the proof. 
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18. RiEMANNIAN RANDOM WAVES 

We have mentioned that the random wave model provides a kind of guidehne for what to 
conjecture about eigenfunctions of quantum chaotic system. In this final section, we briefly 
discuss random wave models and what they predict. 

To define Riemannian random waves, we first consider the standard sphere. We choose an 
orthonormal basis {^PNj}'j=i for Hn- We endow the real vector space Hn with the Gaussian 
probability measure •jn defined by 

\ djv/2 dpf 

^ ' i=i 
Here, dc is ciTv-dimensional real Lebesgue measure. The normalization is chosen so that 
E^^ (?/', ?/^) = 1, where E^^ is the expected value with respect to 7Ar. Equivalently, the 
real variables Cj (j = \,...,d^^ are independent identically distributed (i.i.d.) random 
variables with mean and variance i.e., 

"AT 

We note that the Gaussian ensemble is equivalent to picking ■j/'^r G T^at at random from 
the unit sphere in ?-^7v with respect to the l? inner product. The latter description is more 
intuitive but it is technically more convenient to work with Gaussian measures. In the cutoff 
ensemble, we put the product Gaussian measure n^^j^7„ on 0^=^ 'Hn- 

On a general compact Riemannian manifold (M, g) of dimension m, the analogue of the 
space Tijq oi spherical harmonics of degree N is played by the space 7i/^ of linear com- 
binations of eigenfunctions with frequencies in an interval In := [iV, iV + 1]. The precise 
decomposition of M into intervals is not canonical on a generic Riemannian manifold and 
the results do not depend on the choice. We choose In = [N, + 1] only for notational 
simplicity. In the special case of Zoll manifolds (all of whose geodesies are closed), there is 
a canonical choice where the intervals are centered in the middle of the eigenvalue clusters, 
i.e. the points '^N + ^ where L (resp. j3) is the common length (resp. Morse index) of 
the geodesies. Henceforth we abbreviate Ti.N = Hij^ on general Riemannian manifolds. We 
continue to denote by {LpNjYj=i an orthonormal basis of T-Ln where dN = dimHAr- We equip 
it with the Gaussian measure (11561) and again denote the expected value with respect to 
(^iv,77v) by E^^. 

18.1. Levy concentration of measure. Levy concentration of measure occurs when Lip- 
schitz continuous functions / on a metric probability space (X, d, /x) of large dimension d 
are highly concentrated around their median values M.f- In the fundamental case where X 
is the unit X-sphere 5*^ with the usual distance function, and /i is the SO(A^ + l)-invariant 
probability measure, the concentration of measure inequality says that 



(157) Proh{xeS^ ■.\f{x)- Mf\>r] <ei^^^ 



where ^ /( )| 

^^^^^ ^ d{x,y)>Q d{x,y)\ 
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is the Lipschitz norm. (See, e.g. [L].) 

A key point in the proof is the following well-known and simple fact that the mass the 
sphere concentrates in a small tube around any 'equatorial' sphere. 

Lemma 18.1. Let A G S'^'^~^ C C^, and give S'^'^~^ Haar probability measure. Then 

Proh{P G S^''-^ : |(P,A)| > A} ^ e~^d-i)x^ _ 

Another basic fact about concentration of Gaussian measures is the following Levy con- 
centration of Gaussian measure result due to Sudakov-Tsirelson and Borell: 

Lemma 18.2. |STj Let F C I-Ln and let dj^j = dimT^Tv Let F^p be the p-tube around F. 
Then 

18.2. Concentration of measure and norms. The purpose of this section is to il- 
lustrate the use of concentration of measure inequalities in one of its simplest applications 
to random waves: to determine the asymptotic behavior of norms of random spherical 
harmonics of degree N. The main functionals we consider are the norms on ST-Cn'- 

We only consider the case p = oo in detail; for p < oo see |SZ] . 

Theorem 18.3. For each of the above complex ensembles, there exist constants C > such 
that: 



sn e SHn : sup ^jnI > C^/hgN | < O (^-^^ 



In fact, for any k > 0, we can bound the probabilities by 0{N ^) by choosing C to be 
sufficiently large. 

As a corollary we obtain almost sure bounds on the growth of norms for independent 
random sequences of £^-normalized spherical harmonics. To state the result, we introduce 
the probability sequence space S = Y[n=i'^'^}^ with the measure z/ = nw'=i ^a^- The 
estimate of Theorem 118.31 immediately implies that 

lim sup — I^^n I ^ ^ almost surely . 

N-^oo V log 

Hence we have: 

Corollary 18.4. Sequences of sections ipN € STIn satisfy: 

ll^'A^lloo = 0{^J\ogN) almost surely. 

Results of the latter type were first proved by Salem-Zygmund in the case of random 
trigonometric polynomials on the circle , and by Kahane for random trigonometric poly- 
nomials on tori. Vanderkam |Vanj generalized the results to the case of random spherical 
harmonics by a geometric method that seems special to the sphere. Neuheisel |Neu] adapted 
a method of Nonnenmacher-Voros on holomorphic sections to simplify the proof of the sup 
norm estimates of | Vanj . Here we follow an approach of [SZ] to derive them again from Levy 
concentration of measure. It also gives concentration results on norms: 
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Theorem 18.5. Let 2 ^ p < oo. Then the median values of the C norm on the unit spheres 
SHn are bounded by a constant a = a{p,m), and 

unH^n e SHn ■■ C^ii'N) >r + a}< expi-Cr^N^'^/P) , 

for some constant C > 0. Hence, sequences of sections ipN ^ SHn satisfy W'iPnWp = 0(1) 
almost surely. 

18.3. C°° norms: Proof of Theorem [IQ 

Proof. We refer to §2.31 for notation, but drop the dimensional subscript. Throughout this 
section we assume that H'^^Arll^a = 1. 

We define $jv(x) = {ipN,i, ■ ■ ■ , '^n^m) '■ ^"^ ~^ M.'^^ . It is well known (and easy to prove) 
that this eigenmap is an isometric minimal embedding of (S"™, Qq) as a subsphere of a sphere 
of (M'^^,(is^) where ds'^ is the standard Euclidean metric. In fact, 

(158) nds' = d^fix) ® d^fix) = ^yl^^^^^ go - C^N^+'go. 

Indeed, by SO{m + 1) invariance, the metric (11581) is a constant multiple of the standard 
metric qq. The constant can be calculated using the fact that IIn{x, x) is a constant function. 
By integrating over S"™ one sees that nAr(x, x) = yj^^gm-^ ~ v^I(s™0^™~^' ^-^ taking A of 
this formula one can determine the constants in (I158P (see e.g. ^Chj for background). Then 

(159) IiN{x, y) = Y. <(^)^f(2/) = {^n{x), ^N{y)) . 

Let ipN = (X^kiP = 1) denote a random element of SHn, and write c = 

(ci, . . . , Crfj^). Recall that 



(160) -ipNix) = / IIn{x, y)sN{y)dy = Cjipf (x) = c ■ ^n{x) . 



Thus 

|c ■ ^ n{x) I 

(161) \i'N{x)\ = ||<I>iv(x)|| cos6'a: , where cos6'a; = . 

\\^n{x)\\ 

The angle 9^ can be interpreted as the distance in S'^'^~^ between [c] and $Ar(x). Also, 
||$^(a;)|| = ^Un{x,x) ~ C^M— D/^. 
Now fix a point x G 5™. By Lemma [18. H 

UN{^N:cose,^CN-^"^-'y^y^h^} ^ exp(^-{dN-l)^j^^^ 

(162) = js^-C'N-i^^-^HdN-i) _ 
We can cover S*™" by a collection of kN balls B{z^) of radius 

(163) Rn := ^ 
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centered at points z^, . . . , z^'^ , where 

kN ^ OiR-"") ^ OiN"^"). 

By ffTB2]) . we have 

(164) z/jv jv-JV : max cos ^ CAr-('"-i)/2^/k^| ^ A;;vA^-^'^"*'"~''('^^-i) , 

where denotes a point in X lying above zK 

Equation (11641) together with (I16ip imphes that the desired sup-norm estimate holds at 
the centers of the small balls with high probability. To extend (I164p to points within the 
balls, we consider an arbitrary point G B{z\ A^"™). To estimate the distance 5]^, between 
$Ar(z-') and (^n{w^) in 5''^'^""^ we let 7 denote the geodesic in S*™ from z^ to and use (11581) 
to get 

K < [ ds = I v/$^ = CmN"^ I ds 

(165) ^ C^AT^AT-'" = C'„A^-('™-i)/2 _ 

By the triangle inequality in 5"^^~^, we have 16*2^3 ~ (^yA ^ ^n- Therefore by (I165p . 

(166) cos e^, ^ cos Oy, - 5^ ^ cos Oyj - CmN-^'""^^/^ . 
By ([166]), 

(£+lVbg^ ^ iC + l)Vh^-C^ ^ CVEiN 

and thus 

{^/^iv e 57i:jv : sup cos ^ ^ (C + l)Ar-(™-i)/2 7b^} 

C {sn G 57^iv : maxj cos^^. ^ CA^-("-i)/2yi5^j _ 

Hence by (fTM]) . 

(167) z/^ 1^/;^ G ^Tl^iv : sup cos ^ (C + l)Ar-('"-^)/2 0^| ^ fc^vA^-^'^"''""'^^^^'-^) . 
There exists Am > so that rf^v ~ A^N"^'^, so one has 

^^iv jv^jv G 57i:jv : sup iV^jvl ^ (C + 2)v/logAr| 

Choosing C so that C^Am — rn"^ > is sufficiently large concludes the proof. 

□ 
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18.4. Sup norms on small balls. We use the same method to prove a claim in [NSj on 
the sup norm in small balls on S"^. This corresponds to m = 1 

Corollary 18.6 (Estimate of the maximum). Given p > 0, there exists A such that, for 
any xq E §, P{ max 1/1 ^ A) ^ i 

Again by Lemma [18. H for any x G B{xo, p/N), 

unI^Pn -cose,^ N-^A} ^ exp(^-{2N-2)^'^ . 

We can cover B{xo, p/N) by a collection of kQ{e) of balls B{x^ , j^). By fll68p . we have 



(168) UN 



ipN '■ max COS 



A2 



Again, we need that (11681) implies that the desired sup-norm estimate holds on all of 
B{xq^ jf), so we consider an arbitrary point e B{z^) and estimate the distance 6j^, between 
$iv(x-') and $Ar(?/-^) in S'^'^~^. As before, we have 



and thus 



or 



A A 

cos Oyj ^ — ^ cos e^j ^ — 
^ iV N 



jV'Tv e SHn ■■ supyg(5(^, ^) COS 9y ^ A/N^ 

C {ipN G SHn) '■ niaxj cos^^-j ^ A/N} 



UN {-^N & SHn : sup By ^ A/N } ^ ko{e)e- 



It follows that 



un {ipN ^ SHn '■ sup IV'AfI ^ ^ > ^ ko{e)e 

[ B{xo,l/N) J 

Choosing A sufficiently large will make the right side < |. 

18.5. Relation to Levy concentration. The estimate in Theorem 118.31 is very closely 
related to Levy's estimate. The proof shows that 

(i) £^ is Lipschitz continuous with norm ^^^^^j^^ < ll-C^^IUip < A^(™~i)/2. ; 

(ii) The median of satisfies: A^£^ < CmV^ogN for sufficiently large A^. 

Indeed, Lipschitz continuity follows from equivalence of norms on finite dimensional vector 
spaces. To estimate the Lipschitz norm, we recall the well-known fact that the £^-normalized 
'coherent states' $7(7(2;) = , g^^e the global maxima of on SH^iM^ L^), as follows 

from the Schwartz inequality applied to the reproducing identity s{z) = fj^^ ^n^z, w)s{w)dV{w). 
Moreover, \\^'^{z)\\c^ = ^Un{w,w) ~ A^('"-i)/2. It follows that 

111^1 + S2II00 - llsillool <3Ar(— 
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Now let si have norm < C\/\ogN and let Si = for some w. Then we see that 

jY(m-l)/2 

+ S2II00 - ||si||oo > 

V log 

It obviously follows from (i)-(ii) combined with the Levy estimate (11571) that (for any 
C > 0) 

(169) un{^ G SHn ■■ i:^(^) > Cy/\ogN} < exp{-C{dN - 1) \ogN/2N'^-^). 

Since ^ N"^~^, this is essentially the same estimate as in Theorem 1 18. 3 [ 

18.6. Nazarov-Sodin Theorem on the mean number of nodal domains of random 
spherical harmonics. In this section we review a recent theorem of Nazarov-Sodin [NSJ 
on the mean number of nodal domains of random spherical harmonics on S"^. We closely 
follow their exposition. 

We recall that a nodal domain of an eigenfunction / is a connected component of M\Zf 
where Zf is the zero set of /, i.e. the nodal hypersurface. The number of nodal domains of 
/ is denoted A^(/) The main result of |NS] is: 



Theorem 18.7. There exists a constant a > such that, for every e > 0, we have 



P 



>e} ^ C{e)e- 



-c(e)n 



where c(e) and C(e) are some positive constants depending on e only. 

Remark: As noted above, the total length of the nodal line Z{f) of any spherical harmonic 
/ G Tin does not exceed Cn for some C > [DFj . Since the typical spherical harmonic has 
~ an^ nodal domains, the perimeter of most of its nodal domains must be of order ^ and 
the diameters of the nodal domains must be comparable to ^. 

The proof of Theorem 118.71 consists of the following steps: 

(1) Proof that the lower bound EA^(/) ^ const n^. 

(2) Proof of exponential concentration of the random variable N{f)/n^ around its me- 
dian. The proof uses the uniform lower continuity of the functional / ^ N{f) with 
respect to the L^-norm outside of an exceptional set E G Ti. of exponentially small 
measure and Levy's concentration of measure principle. 

(3) Proof of the existence of the limit lim KN{f) /n^ . In this part, we use existence of 

n— >oo 

the scaling limit for the covariance function E|/(a:;)/(?/)}. 

18.6.1. Lower bound for 'EN(f). The first step is to show that EA^(/) > n'^. By Corollary 
118. 6[ for any p > 0, there exists A > such that, for all x G S, 

P| max I/I ^ A) ^ i. 

Furthermore, the zonal spherical harmonic b^ with pole at x, i.e. b^ = ^^^-''^^ , has the 
following property: there exists p and ci such that, for each sufficiently large 



%\\ = 1, b^{x) ^ CiViV, and &a;|g2?(^,p/iv) ^ -Ci'/N ■ 
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Proof of the lower bound for N{f): Fix x G S. We have / = C,obx + fx where is a 
Gaussian random variable with E^q = jw+i^ ^ Gaussian spherical harmonic built 

over the orthogonal complement to in ?ijv and normalized by E||/^. |p = ^n+i ■ We choose 
a Gaussian random variable ,^0 independent of ,^0 and of f^ with E,^q = E,^q = ^nj^-i ? and set 

f± = ±^obx + fx- These are Gaussian spherical harmonics having the same distribution as 
/. Note that 

/ = eo&x + ^ (/+ + /-) , 

and that by Corollary 118.61 

P({max ^A}n {max |/_ | ^ A}) ^ 1 - + ^) = ^ . 

Now, consider the event Vt.j, that f{x) ^ A and /|gx,(^,.p/jv) ^ Clearly, Vl^ must contain 
a nodal component of /. The event Vt^ happens provided that 



^oViV ^ 2c^^A and max |/±| ^ A . 

V{p,x) 



Since the variance of the Gaussian random variable ^o\N is of constant size, there exists 
K > such that 



P(fi^) ^ P(CoviV ^ 2c^M) -PffmaxIZ+l ^ A} n {max |/_| ^ A}) ^ k > Q . 

'D{p,x) V(p,x) 

Now pack S"^ with ~ A^^ disjoint disks Dj of radius 2p/N centered at points Xj. Each 
of them contains a component of Z{f) with probability at least k, and since the discs are 
disjoint, A^(/) > Ylj -^if'i -^j) where N{f,Dj) is the number of nodal components in the 
disc Dj. Hence, 

EiV(/) > Y,^N{f,D,) > 5^P(1^.,) > a,p N'. 

j j 

□ 

18.6.2. Exponential concentration near the median. This exponential concentration follows 
from Levy (-Sudakov-Tsirelson) concnentration of measure plus the delicate fact that A^(/) 
is stable under slight perturbations of / for all but an exponentially rare set of /, called 
'unstable' spherical harmonics. 

Lemma 18.8. For every e > 0, there exists p > and an exceptional set E C TIn of 
probability F(E) ^ C{e)e~^^'^^^ such that for all f E Hn \ E and for all g G T-Cn satisfying 
\\g\\ ^ p, we have N{f + g) ^ N{f) - eA^^^ 

We refer to [NS\ for the proof. Lemma 118.81 in conjunction with the Sudakov-Tsirelson 
concentration of measure theorem (Lemma 118.21) implies an exponential concentration of 

near its median a^. Recall that A^p is the p-tube around A. 
Corollary 18.9. P{/ g Hn- \N{f) - a^l > eN^} ^ 2e-^^'^. 
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Proof. Let F = {f e Hn: N{f) > {a^ + e)N^}. Then for / G (F \ E)+p, we have 
A^(/) > ajy N"^, and therefore, P((F \ -E)+p) ^ |. By the Sudakov-Tsirelson inequahty 
(Lemma [TE2D, P(i^ \ ^) ^ 26"'='''^, hence 

P(F) ^ 2e'^'''^ + C(e)e-'=(")^ ^ C(e)e-'^(^)^ . 

On the other hand, let G = {f eHn- N{f) < (a^v - e)N'^}. Then 

G+p C {/ G : iV(/) <a%}UE 

and so 

P(G+,) ^ i + C(e)e-(^)^ < ^ 
for large A^. It follows again that P(G') ^ 2e~'^^^^ for large A^. 

□ 

Further, it is proved in [ NSj that 

Proposition 18.10. The sequence of medians {ajy^ converges. 

Since the random variable N{f)/N'^ exponentially concentrates near its median and is 
uniformly bounded, it suffices to show that the sequence of means KN{fN)/N'^ converges; 
the sequence of medians {ajy^ converges to the same limit. The authors of [NSj show that 
{EA^(/„)/n^} is a Cauchy sequence. 

19. Appendix on Tauberian Theorems 

We record here the statements of the Tauberian theorems that we use in the article. Our 
main reference is |SVj . Appendix B and we follow their notation. 

We denote by the class of real- valued, monotone nondecreasing functions A^(A) of poly- 
nomial growth supported on M+. The following Tauberian theorem uses only the singularity 
at t = of dN to obtain a one term asymptotic of A^(A) as A ^ oo: 

Theorem 19.1. Let N e and let ip g S{R) satisfy the conditions: ijj is even, iplX) > 
for a// A G M, G C^, and V^(0) = 1. Then, 

ilj*dN{\) KAY =^ \N{\)-N*^{\)\<CAX', 

where C is independent of A., A. 

To obtain a two-term asymptotic formula, one needs to take into account the other singu- 
larities of dN . We let ip be as above, and also introduce a second test function 7 G 5 with 
7 G and with the supp 7 C (0, 00). 

Theorem 19.2. Let A^i, A^2 G and assume: 

(1) iV,*^(A) = 0(A-),(j = l,2); 

(2) N2*i^{\)=Ni*ij{\) + o{r)- 

(3) 7 * dN2{\) = 7 * dNi{\) + o{\r 

Then, 

iVi(A - 0(1)) - o{r) < N^iX) < N,{X + 0(1)) + oiX"). 
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This Tauberian theorem is useful when the non-zero singularities of dN2 are as strong as 
the singularity at t = and N2 does not have two term polynomial asymptotics. 
Alternatively, we may use the Tauberian lemma of |HoI-IV] : 

Lemma 19.3. Suppose that ^ is a non- decreasing temperate function satisfying /i(0) = 
and that u is a function of locally hounded variation such that z/(0) = 0. Suppose also that 
m ^ 1 and that ip G i5(M) is a fixed positive function satisfying J {p{X)dX = 1 and 0{t) = 0, 
t ^ [—1, 1]. If (p„[\) = a~^ip{X/ a) , < (J ^ (Jo; assume that for A G M 

(170) \dv{X)\ ^ (Ao(l + |A|r + ^1(1 + |A|)"^-^) dX, 
and that 

(171) \{{d^i - du) * <^^)(A)| ^ 5(1 + |A|)-^ 
Then 

(172) IMA) - u{X)\ ^ C^{Aoa{l + \X\r' + A.ail + \X\r'-' + B), 
where Cm is a uniform constant depending only on ctq and our m ^ 1. 
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